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Buamsanne (WHaMM4eCKHX CBOMCTB B3AaMMOJeHCTBYIONINX NOACHCTeM Ha 9BOJUOWHIO POPMHPOBAHHA 


RR 
H30upaTeJIbHOrO TepeHoca B y31ax TpeHuA 


B. JI. 3akoBoporupiii', B. E. TBnnpKnans’, IL. C. Konoaknn? 


1,2,3 


Introduction. Selective transfer is a typical example of the 
self-organization processes in tribosystems. In this case, joint 
surfaces of the servovite film are formed in the contact area, 
which changes fundamentally the friction and wear conditions. 
To form selective transfer in the area of mating surfaces, some 
power of the irreversible transformations of the input energy is 
needed, which depends on the elastic-dissipative properties of 
the contact surfaces. 

Materials and Methods. The mathematical model of the dy- 
namic system considering the evolutionarily changing servo- 
vite film is given. Its formation depends on the phase path of 
the irreversible transformations power in the area of surface 
matching, and it is represented by the Volterra integral opera- 
tor of the second kind. 

Research Results. The outcome analysis including dependenc- 
es of the servovite film formation on the dynamic parameters 
of interacting subsystems is provided. The mathematical simu- 
lation of the evolution of the friction unit properties with the 
formation or destruction of the servovite film is first consid- 
ered. 

Discussion and Conclusions. During the evolutionary process 
of the servovite film formation, the dynamic coupling parame- 
ters generated in the friction unit, change. Consequently, the 
dynamic properties of the system also change. The tribosystem 
dynamics is first considered under the process of forming the 
selective transfer. 


Keywords: dynamic friction system, selective transfer, servo- 
vite film, evolution. 
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Jlouckoii rocyyapcTBeHHbIit TeXHH4YeCKHM yHuBepcurTert, Poctos-Ha-Jlony, Poccniickaa Degepauna 


Beedenue. Tunw4HbIM MpHMepoM upoleccos camMoopraHH3a- 
UHH B TpHOocHcTeMax ABIIAeTCA U30upaTesIbHbIii MepeHoc. B 
3TOM Cily4yae B OOACTH KOHTaKTa (OpMUpyroTca compsarato- 
UHe MOBepXHOCTH CepBOBUTHOM MJICHKH, YTO IIPHHUMMMAIbHO 
MeHAeT YCIOBHA TpeHHaA HM W3HalMBaHHa. Ja OOpa3z0BaHHa 
w30upaTeIbHOrO HepeHoca B 30HE COMpmKeHHA KOHTAaKTUpy- 
FOWIMX MOBepXHOCcTeii HeOOXOAMMa HeKOTOpad MOLIHOCTb 
HeoOpaTuMBbIx MpeoOpa30BaHHit MOABOAMMON 3HEprun, KOTO- 
pad 3aBHCHT OT yipyro-4MCCHMaTHBHbIX CBOMCTB KOHTAKTH- 
PY!OWHx MOBepxHoctTei. 

Mamepuaavi u memodoi. Tipusoqutca MaTeMaTHueckad MO- 
Wesb WHHaMM4eCKO CHCTeMBI C YYCTOM 3BOJIFOIMOHHO H3Me- 
HArOUlelicaA CepBOBUTHOM WieHKH. Ee PopMupoBaHHe 3aBHCHT 
oT (pa30B0ii TpaeKTOpHH MOUHOCTH HeoOpaTuMBIx Mpeodpa- 
30BaHHii B 30He COMps»KeHHA MOBepXHOcTeli HM MpescTaBsAeT- 
cA B BHe MHTerpalbHoro onepatopa Bosbreppbl Broporo 
poyla. 

Pezynomamet uccaedoeanus. Tipwpoyutca aHalM3 pe3yibTa- 
TOB HCCJIeOBaHHA, B TOM YMCIIe 3aBMCHMOCTH (OpMUpoBaHHA 
CepBOBHUTHOM MIeCHKH OT JHHAaMM4eCKHX MapaMeTpOB B3aHMO- 
WelicTByroulMx MoycHctem. Bueppble paccmatpHBaetca mpo- 
OsleMa MaTeMaTH4YeCKOrO MOeNMpOBaHHA IBOIFIOUMM CBOLCTB 
y3Jla TpeHHa c OOpa30BaHHeM WIM pa3spylleHHeM CepBOBUTHON 
TWICHKH. 

O6cyarcdenue u 3axmouenue. B xoqe 3BOMIOWMOHHOTO Tpo- 
yecca o0pa30BaHHaA CepBOBHTHOM MJIeHKH MeHAFOTCA Mapa- 
MeTpbI JMHaMH4ecKoli cBA3H, POPMUpyeMoii B y3JIe TpeHHA. 
CreqoBaTesbHO, MeHAFOTCA HM MHaMMuecKHe CBOliCTBa CH- 
cTembl. BnepBble npoaHasM3upoBaHa JHHaMuKa TpHOocucTe- 


MBI B Iipolecce oOpa30BaHHa vw30upaTesIbHOrO TlepeHoca. 


KurroueBbie c1I0Ba: JMHamMMyeckas CHcTema TpeHHA, w30upa- 
TeJIBHBIM TlepeHoc, CeCpBOBHTHad MWICHKa, SBOJIIOUHA. 


O6opa3zeu Ona YUMupoeanua: 3aKoBopoTHBIi, B. JI. Bausanue 
TWHaMHYeCKHX CBOMCTB B3aMMOJeHCTBYIOMMX MOACHCTeM Ha 
SBOMOUMIO POpMUpOBaHHA W3ONpaTebHOTO WepeHoca B y3- 
max tpenua/  B. JI. 3axopopotusml, B. E. [puayxuma, 
II. C. Konogxun // Bectuuk JJon. roc. TexH. yH-Ta. — 2019. 
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Introduction. After publication of the papers by I. Prigogine [1-3] and G. Haken [4, 5], lots of issues on the 
operation of technical systems interacting with different media are considered from the view point of their self- 
organization [6-12]. If you follow the synergetic analysis paradigm, then when studying such systems, you should first 
perform the procedure for extending the state space dimension [7, 8]. That means that it is necessary to additionally 
consider a model of environment in the coordinates of the system state. There is an exchange of information, material, 
energy, etc., with the environment. Here, coherent interactions of various physical nature are possible, that is, synergetic 
phenomena. A typical example of self-organization is the effect of selective transfer discovered by I.V. Kragelsky and 
D.N. Garkunov in 1956 [13]. 

The ideas of self-organization were formulated by B.I. Kostetsky [14] in the 60s of the XX century, and then 
developed in papers on structural adaptability together with L.I. Bershadsky [15] and N.A. Bushe [16]. In a number of 
papers, the formation of a servovite film in the tribosystem when it enters the selective transfer mode is studied [17-22]. 
It is shown that definite tribochemical reactions are necessary for the formation of selective transfer, as well as some 
power trajectory of irreversible transformations in the system-environment interface [17—24]. When entering the selec- 
tive transfer mode, self-oscillations are observed. In some cases, chaotic attracting sets of deformation displacements of 
the contacting pairs are formed [17, 22—24]. 

Materials and Methods. Thus, when studying the transition to the selective transfer mode, it is necessary to 
consider the tribosystem dynamics in the unity of the subsystems interacting through friction and the dynamic coupling 
formed by the friction node. The parameters of such a dynamic link and the formation of a servovite film depend on the 
power path of irreversible transformations on the work done. The mathematical modeling of such an evolutionary sys- 
tem is considered, and the dependence of evolution on the dynamic parameters of interacting subsystems is analyzed. 

Research Results 

Mathematical system modeling. In mathematical modeling, we use the previously obtained results. In [17, 
23, 24-26] it was shown that the basic dynamic properties of the friction system can be disclosed on the basis of the 
following assumptions: 

- the sample is absolutely tough; 

- indenter deformations are considered in the plane normal to the contacting surface and passing through the di- 
rection of the relative slip velocity. 

When studying the dynamics, we can limit ourselves to the first vibration modes. Then, the system model is 
represented by the equation 

2 
ma Kh 40x = FeO, (1) 
dt dt 





m 0 my Ay, Ci C24 ae , ; en . 
where m = , A= CS are positively defined, symmetric matrices of inertial, velocity, 
0 m hn hop C2 €22 


and elastic coefficients; X ae 8 we Cr is vector of the deformation displacements of the indenter tip; 


BO=fhrO.A> (17 is vector of forces affecting the indenter tip (Fig. 1). 





Fig. 1. Dynamic friction system schematic 


We will follow the synergistic concept of the analysis, that is, we will represent Fy (ft) forces in the coordi- 


nates of the state. To do this, we introduce the concept of tribo-medium - this is the third body that is formed in the tran- 
sition zone between the contacting surfaces. When approaching, the conditions of their interactions change, which phys- 
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ically manifests itself in changing parameters such as the actual contact area, heat production, friction coefficient, mo- 
lecular interaction conditions, diffusion processes, etc. [27-30]. However, the primary reasons for all changes are inter- 
actions due to the contact mechanics. 

In this case, due to changing the actual contact area, the normal pressure forces grow disproportionately fast, 
which prevent the surfaces from approaching. It is convenient to consider modeling these forces [23, 24] as the follow- 
ing function of approach 

F(X) = Fo exp[-a(4y)]-U,  X1 € (0, +00), (2) 


where © is coefficient of the contact force rate in | mm | ; U isan external force. 


The above papers show that when a servovite film is being formed, a potential barrier appears in the function of 
approach, which actually determines the carrying capacity of the servovite film of the friction unit. The formation of a 
potential barrier can be conveniently represented as 


2 
D(X) = D9 exp[-0 (X1 — X19)", (3) 
where a, is the parameter characterizing the potential barrier slope in | mm]; X19 18 the potential barrier coordi- 


nate; @, 9 is the evolutionary parameter in [kg]. Thus, the function of convergence is Fi y(X 1) =®(X1) + (41) sum. 
On modeling the tangential component of the force, let us take into account that when forming a servovite film, the fric- 
tion force is by an order of magnitude less. Therefore 
1 
KOR y(t), at: X1(X1 9; 


(2) (4) 
kp Fis (t), at: X1)X\ 95 


Fy (@)= 
1 2 
where &0)y)kO) | 
For example, in the steel — glycerin — brass friction system, after the formation of a servovite film, the friction 
coefficient decreases by more than an order of magnitude. Previously it was shown that the forces of contact interaction 
are characterized by a lag in the tangential component of the force (friction force) with respect to the normal pressure 
forces. In this case, instead of (4), the following relation occurs 
1 1 
kK Fy @-T), at: XX 93 


(5) 
KF 5 (¢-T), at :X1)X1 0, 


FyxO= 
where 7) ow () are time-lag constants depending on the properties of the tribo-medium formed in the contact area, 


and on the relative slip velocity. In (5), T Oyyr 2) condition is always satisfied. In addition, T and 7) decrease 


with increasing the relative slip velocity. Later, the hypothesis was adopted, according to which the potential barrier is 
formed depending on the work paths and the power of irreversible transformations as the indenter moves relative to the 
sample in the direction of the relative slip velocity. Moreover, ® 9 is affected not only by the current value of power, 


but also by the preceding values of power. Therefore the following is true: 


t 
P(N) =P] we-E)N(EdE, (6) 
0 


. : : : 1 : 
where N = F5(V +dX,/dt) is the power of irreversible transformations; w(t—&) = {B ovl-70- Or is the core of 
the integral operator considering the effect of the preceding power values; T is the parameter in [s] characterizing a 


long-time impact of power on the potential barrier; { is the parameter in [m"| ? 


The equations (1) - (6) characterize the mathematical model of a dynamic friction system considering the evolution of 
its properties during the formation of a servovite film. 

Dynamic properties change depending on the potential barrier parameters and the dynamic properties 
of the indenter. First, consider the properties of the frozen system on the assumption that the parameters of the servo- 
vite film are preset and constant. Suppose also that 7 ® =0 and T° =0. Consider an example of changing system 
properties depending on the stages of the evolutionary transformation (Fig. 2). 
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dX,/dt, mm/sec 


120 
100 
80 
60 
40 


Wi = 0.22 um 





—0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 Xi, um 
Fig. 2. Transformation of phase path projections onto X; — dX;/dt plane 


At the initial stage of transformations (Fig. 2-1), a unique equilibrium point is formed in the system, which cor- 
responds to the friction with a traditional coefficient in the range of 0.2—0.3. Moreover, this point has the property of 
global attraction. Then (Fig. 2-2), two equilibrium points (Xj;,X1.2) are formed during the formation of a servovite film 
at the initial stage. The domain of attraction of Xj, point is bounded by a saddle-shaped separatrix. The rest of the do- 
main is characterized by the attraction to X/ point. It is important that at X,, point of equilibrium, a servovite film is 
already formed, but its value and properties (determined by ® ¢ ) do not allow for the attraction of the trajectories of the 


whole phase space. Therefore, depending on the initial conditions or disturbances (for example, fluctuations or formable 
attracting sets of deformation displacements in variations relative to the equilibrium point), it is possible to generate 
properties of a friction unit that differ fundamentally from each other. Subsequently (Fig. 2-3), the domain bounded by 
the saddle separatrix decreases and disappears (Fig. 2-4). In the latter case, the servovite film thickness and its carrying 
capacity increase. In the model, this is shown as an increase in the potential barrier, and the second equilibrium point 
acquires the properties of global attraction. In this case, a stable selective transfer is formed in the system. 

The illustrations characterize the system dynamics considering that the servovite film properties are frozen, 
that is, they do not evolve. However, even in this case, the domains of attraction of equilibrium points corresponding to 
the friction with a servovite film and without it depend fundamentally on the dynamic parameters of the indenter, pri- 
marily, on the matrices of its stiffness (c) and dissipation (/) (1). In actual practice, there is a time evolution of the 
system. It depends on the relative slip velocity 

(Vo ), the force (U ), the initial state of the surface of the contacting bodies, etc. The evolution simulated by an 
integral operator (6) is affected by many physical factors the integral count of which is determined by the kernel of the 
operator. Therefore, when varying the external conditions (for example, relative slip velocity), it is necessary to consid- 
er the velocity — integral operator parameters relationship, as well as the velocity — the power trajectory of irreversible 
transformations relationship. It is experimentally shown that with increasing speed, 7;,i=1,2 parameter in (6) decreas- 
es. In addition, frictional forces (and, consequently, the power path of irreversible transformations) depend on the rela- 
tive slip velocity. Note that it is the time required to establish a stationary servovite film that determines one of the indi- 
cators which allow identifying the parameters of the integral operator kernel. 

Cursorily, we give an example of changing the phase trajectories of deformation displacements in Xj -direction 


as the indenter dynamic parameters and the tribological medium properties vary (Fig. 3). 
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Fig. 3. Transformation of phase path projections onto X, — dX;/dt plane depending on system parameters and friction conditions 


8 
Mechanics 


First, consider friction with a sufficiently high velocity of relative slip Vg = 2.0 u/c. In this case, when form- 


ing a servovite film, the power of irreversible transformations decreases and the destruction is not observed, since the 
relative slip velocity is rather high (Fig. 3, a). At the initial stage, under various conditions, the trajectories are asymp- 109 
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totically attracted to the equilibrium point without a servovite film (point A in Fig. 3, a). After the potential barrier has 
been established, all trajectories evolve to AY equilibrium point on the servovite film. Fig. 3, b, c, and d show evolu- 
tionary curves with one initial value corresponding to B point. At this, the deformation displacements and the corre- 
sponding properties of the friction system deform to one of the two equilibrium points. However, evolutionary curves 
differ depending on the properties of the tribo-media and the Q-factor of the oscillatory circuits of the indenter subsys- 
tem without friction. Fig. 3, b shows an example of the evolutionary trajectory under a conditions, but at the relative 
slip velocity of Vy =1.0 m/s. In this case, the power released in the friction zone is not enough to maintain the servo- 


vite film, and friction is periodically observed in the system with a servovite film and without it. In the deformation 
displacements, this causes the effect of low-frequency oscillations. Specifically, the formation time and the destruction 
time of the servovite film differ essentially. In Fig. 3, this is indicated through a significant decrease in the rate at which 
the indenter tip returns from A" point to A point. 

The Q-factor of the indenter subsystem can be increased (Fig. 3, c). In this case, complex oscillatory displace- 
ments are formed in the system. Some of their components characterize the low-frequency material exchange between 
the lubricant and the friction surface, as well as vibration displacements at the frequencies of the indenter oscillators. 

In the direction orthogonal to the relative slip velocity, cyclic forces are generated. Due to this, and also due to 
the dynamic interaction, a short-term formation of forces exceeding the potential barrier is possible. The situation is 
even more complicated if we additionally consider the lag in the variations of the tangential force components to their 
normal components (Fig. 3, d). Depending on the lag, various attracting sets of deformation displacements and, conse- 
quently, the forces of contact interaction are formed. The generation of chaotic dynamics in the system, which introduc- 
es great uncertainty in the regularities of servovite film formation, is also observed. 

Discussion of the results. It is known that moving of the tribosystem to the selective transfer mode is deter- 
mined by tribochemical reactions and the material exchange between the contacting surfaces and the lubricant. In order 
for the tribosystem to enter the selective transfer mode with the formation of a servovite film and to maintain it in the 
process of friction, some particular power of irreversible transformations of the input energy of the mechanical system 
is required [17, 21-24]. During the evolution of the dynamic friction system, the carrying capacity of the servovite film 
varies, which is simulated by the size of ®, g potential barrier preventing direct contact of the friction surfaces. If the 


forces of contact interaction in the direction normal to this surface exceed ®, 9, then a dynamic restructuring of the fric- 


tion system with partial destruction of the servovite film is observed. Therefore, the friction conditions in the selective 
transfer mode depend not only on tribochemical reactions and the power of irreversible transformations, but also on the 
dynamic mode. Here, the stability of the equilibrium point and various attracting sets of deformation displacements 
(limit cycles, invariant tori, chaotic attractors) formed in its neighborhood are of fundamental importance. In addition, 
the dynamic system is perturbed, for example, by the out-of-true-running. Moreover, due to a fundamental change in the 
friction conditions during the direct contact and at the contact through a servovite film, the power of irreversible trans- 
formations in the surface matching changes, which changes the evolution and the maintenance conditions of the formed 
servovite film. 

Conclusion. So, in conclusion: the stability of selective transfer is affected by both the dynamic parameters of 
the subsystems interacting through the friction node, and by disturbances that always occur in a real machine. There- 
fore, when studying selective transfer to ensure the friction unit wearlessness, it is required to obtain identity of the dy- 
namic systems on a friction machine and in a real machine. For this purpose, the well-known techniques [27, 31] can be 
used. 

The data show that the formation of a servovite film in a tribocontact node is affected not only by the tri- 
bochemical characteristics of the mating surfaces, but also by the parameters of the indenter and sample subsystems 
interacting through the tribology. Therefore, in each specific case, there is a bounded domain of parametric space of the 
dynamical systems, as well as disturbances in which the formation of selective transfer is a rule. 
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Introduction. Functional properties of diamond powders are 
determined by a large-scale structural factor since it affects 
the formation of structurally sensitive mechanical properties 
— stress limit and yield value. Considering the qualitative 
correlation between yield value and hardness, it is possible to 
predict an increase in hardness including highly rigid materi- 
als. 

Materials and Methods. Physical characteristics of the basic 
types of fillers that make up the reinforcers are considered, 
systematized and tabulated. M2-01 tin bronze (20 wt. % tin, 
80% copper) was used as a bond. Ultradisperse natural dia- 
mond (UDND, 0.5—4 wt. %) was added to it, as well as pow- 
ders of natural diamond (3/2 ym fraction, 7/5 um, —-40 pm) 
obtained through processing diamonds at the enterprise of 
“Sakha Daimond” JSC. The above materials were made on 
the crushing and screening equipment and shaking tables. 
The stages of obtaining powders were recorded using the 
raster electron microscopy. Vibroscreens were applied for 
the grain-size classification of diamond powders. Physical 
and mechanical characteristics of the produced samples were 
tested by standard methods. VLTE-500 electronic fourth- 
class laboratory balance was used for weighing. Density was 
determined by MK 0—25 mm micrometer according to GOST 
6507-78. 

Research Results. Porosity was calculated through the actual 
and theoretical density. It was found that with a decrease in 
the filler size, an improvement in the physicomechanical 
properties of the binder modified with a diamond powder is 
observed. The best performance was observed in the samples 
with the UDND filler. 

Discussion and Conclusions. As a result of the study, it was 
recorded that the calculated data differ from the experimental 
data since they show an increase in the material hardening 
pro rata to the amount of the diamond particles introduced 
into the volume. An assumption has been made that the con- 


sidered hardening model (Orowan model) does not take into 





* The research is done within the frame of the independent R&D. 


Beedenue. DyHKYHOHAbHbIC CBOMCTBAa aJIMa3HbIX MOPOWIKOB 
OOYCIOBIICHbI MaCIUTaOHBIM CTpyKTypHbIM (bakTOpoM, T0- 
CKOJIbBKy OH BJIMAeT Ha (opMMpoBaHHe cCTpyKTypHO- 
YYBCTBHTCJIBHBIX MCXaHHYeCCKHX CBOMCTB — lmpeeoB 
IIpO4HOCTH HW TekyyecTH. YYMTbIBAad KAYCCTBeHHY!O KOppella- 
WHIO MeXK ly WpesesioM TeKyyecTH H TBepJOCTbIO, MO2XKHO 
IIpOrHO3SHPOBaTb MOBBILIeCHHe TBEPAOCTH, B TOM 4HCIIC BbICO- 
KOTBepJIBIX MaTepHasIOB. 

Mamepuaavi u memoovi. Paccmotpeubl, CHcTeMaTH3MpOBaHbl 
MW TipeycTaBeHbI B Bue TaOMMIbI (u3sH4eckHe xapakTepH- 
CTHKM OCHOBHBIX THIOB HanojIHHTeseH, BXOJAMJMX B COCTAaB 
yupounutesieH. B KayecTBe CBA3KH HMCHOJIb30BaHa OJIOBAHH- 
ctat Opou3a M2-01 (20 mac. % omosa, 80% mMegn). B Hee 
WoOaBIAIH yiIbTpayqucnepcubli mpupogzHEi anma3 (YJIIA, 
0,54 mac. %), a TakoKe MOpPOMIKH MpHpoyvoro asima3a 
(ppaxuuu 3/2 mKM, 7/5 MKM, —40 MkM), HosyyeHHbIe mpu 
TepepaOotke ayimMa30B Ha UpequpuatuHu OAO «Caxazaii- 
MOH)». Ha3BaHHble MaTepHaJibl W3rOTaBIMBaJIMCh Ha JIpo- 
OMWJIbHO-KaCCHUKalMOHHOM OOopyoBaHHuH Hu BHOpOCTO- 
max. CTaqwM MouyYeHHA NOpowKoB @uKcHpoOBaIIMch c M0- 
MOLIbIO pacTpoBOH 3IeKTpOHHOM MuKpockonMu. Jia KyIac- 
CHuKalHM aIMa3HbIX MOPOWKOB 10 3epHHCTOCTH MpHMeHsA- 
JM BHOpocuta. Du3sH4eckHe M MexaHHueckve XapakTepHcTu- 
KH H3FOTOBJICHHBIX OOpa3i{OB HCIIBITbIBaIM 110 CTaHapTHbIM 
MeTOquKaM. J[y1f B3BCLUIMBAHHA HCHONb30BaNIM JabopaTopHBEle 
3JICEKTPOHHEIe BeCbI 4YeTBeproro Kiacca BJITS-500. Inor- 
HOCTB onpeyeuasH MK 0-25 mm 
m0 TOCT 6507-78. 

PesynbTaTbl UccneqoBaHua. Uepe3 dakTHyeckylo HU TeopeTH- 


MHKpOMeTpOM 


YeCKYIO IIIOTHOCTH pacc4untTaHa MopucTocTs. BsiacHHsocb, 
4TO C YMeHbINeHHeM pa3Mepa (ppakUHH HamoJHHTesIA 
HaOOfaeTcA yilyalleHHe (H3HKO-MexaHW4eCcKHX CBOMCTB 
CBA3KH, MOAM@HUMpOBaHHO! aJIMa3HbIM MOpOMIKOM. 
Hannyuutue moKa3aTesIM OTMeYyeHbI y OOpa3L[OB C HallosHH- 
Temem u3 YJITIA. 

O6cyarcdenue u 3akxmouenua. B pe3ymbTaTe MpoBeyeHHoro 
HCCIeHOBaHHA OTMC4CHO, YTO pacueTHble aHHbIe OTIM4AIOT- 
CA OT 9KCIePHMeCHTAJIBHBIX, T. K. HOKa3bIBaIOT yBelIM4eHHe 
yIpOuHeHHA MaTepHasla MpONOpUMOHaIbHO KOJIM4eCTBY BBO- 
MMBIX B OObeM YacTHI aiMa3a. BelqBHHyTO MmpesMos0x«Ke- 
HHe, UTO paccCMOTpeHHad MOJeNIb yIpodHeHHA (Moses Opo- 


* E-mail: marisafon_2006@mail.ru, arkhangelskaya@yandex.ru, Fedot_andrey@mail.ru 
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account the formation of carbon and the agglomeration of BaHa) He YUHTBIBaeT OOpaz0BaHHe yryiepoya HW aroMepal{Hio0 


diamonds into larger objects in the matrix volume under an ayIMa30B B Ooylee KpylHble OObeKTHI B OObeMe MaTPHUbI Tp 
TIOBIMICHHH KOJIM4CCTBAa BBOJMMBIX ajIMa30B. Ecru oObeM 
yacruy YA yocturaet 3 %, B MaTepHasie pacTeT cojlep- 
2%KaHHe yruepoga. B pe3syibTaTe YacTHIbI HaNOHUTeIA MOsI- 


increase in the number of input diamonds. If the UDND par- 
ticle volume reaches 3%, the carbon con 


tent in the material increases. As a result, the filler particles 
HOCTBIO HE OKHCJIAFOTCA, TEM CaMBbIM YBCJIM4MBad KOJIM4eCTBO 


are not fully oxidized, thus increasing the number of pores in 
Top B MaTepualie. 


the material. 


Keywords: bond, metal matrix, composite, hardener, ul- K.novesBble C10Ba: CBA3Ka, MeTaJWIM4ecKad MaTpHia, KOMIIO- 

trafine particles, hardening mechanisms. 3HT, YNPOUHHTelIb, yIbTpayMchepcHble 4acTHIbl, MCXaHH3MbI 
ynpouHenna. 
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Introduction. It is known that the physicomechanical properties of diamond powders are determined by a 
larger specific surface area and grain dispersion. In particular, this is shown by the Hall — Petch equation which is per- 
formed over a wide range of grain sizes (up to 1 um). Functional properties of these materials are determined by the 
large-scale structure factor since it affects the formation of structure-sensitive mechanical properties — stress limit and 
yield value [1, 2]. Given the qualitative correlation between yield strength and hardness, it is possible to predict an in- 
crease in hardness including highly rigid materials: the finer the filler, the fewer defects it has and, accordingly, the 
higher the strength. 

Materials and Methods. Strengthening mechanisms in highly rigid materials depend on the interaction pat- 
tern of the introduced particles or hardener fibers and the matrix material. A successful application of the dispersion 
strengthening effect is shown in [3—6, 7]. Under such strengthening, a structure that makes the dislocation motion hard- 
er is created in the materials. Discrete particles of the second phase characterized by high strength and melting point 
inhibit the dislocation motion especially strongly. Considering the two-phase structure and high hardness of the materi- 
als obtained, it should be expected that their wear resistance will also be higher than that of the nonhardenable matrix 
materials. Besides, the strength of the components interface is an important factor for the wear resistance of composites. 
A strong adhesive bond at the interface ensures a composite with high rigidity and higher static strength [3, 6-8]. 

According to the mechanism described by Orowan, ultrafine diamond particles distributed in the binder vol- 
ume affect the strain hardening of the composite material. Particles of the particulate filler introduced into the matrix 
inhibit the dislocation motion in a metal increasing its strength at the standard and elevated temperatures. In addition, 
they represent a mechanical obstacle in the path of crack propagation that may appear in the matrix, and increase the 
fracture resistance of the composite material. 

Another characteristic that determines the reinforcing filler — matrix relationship is the thermal linear expan- 
sion coefficient. For solids at constant pressure and temperature, the thermodynamic equilibrium criterion is Gibbs min- 
imum potential (or energy). This value shows an energy change during a chemical reaction and demonstrates the possi- 
bility of chemical reactions between the material components [5]. Thus, a minimum change in Gibbs energy corre- 


sponds to a stable equilibrium between the components of the system (Table 1). 
Table | 


Physical characteristics of the main types of fillers [9, 10, 11] 





























pubstanes AG" sp jor Stas imel Microhardness, x10? Mpa __| Temperature stability, °C 
and state Gibbs energy change 
C (diamond) 2.377 1000 650-700 
W 0 258 3300-3400 
Al,O3 —1582.3 180-220 1500-1700 
BN —226.8 800-900 1100-1300 
SiC —60 300-320 1200-1300 
BeO —579.9 152 2500 
Be,C —948 780 2150 
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The table shows that the optimal filler is diamond. It has a fairly low value of Gibbs energy change, the highest 
microhardness value, but the lowest temperature stability. Diamond has a high adsorption capacity [12] and is the least 
chemically active compared to other forms of carbon. These properties are important advantages when using diamond 
as a hardener. 

The work objective is to study features of the mechanisms of forming the metal-matrix composition structure. 

Objects of Study. We used M2-01 tin bronze (20 wt. % tin, 80% copper) as a basic binder. Ultradisperse natu- 
ral diamond (UDND, 0.5—4 wt. %) was added to it, as well as powders of natural diamond (3/2 um fraction, 7/5 um, — 
40 um). 

Research Methodology. Natural diamond powders (NDP, 3/2 um fractions, 7/5 um, —40 um) and the UNDN 
submicropowder were obtained during the processing of diamonds at the enterprise of “Sakha Diamond” JSC. The 
powders were made on the crushing and screening equipment and shaking tables under the optimum conditions. Fig. 1 
shows the stages of obtaining the powders. The images were taken using the scanning electron microscopy (SEM). 


i ” / f 
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Fig. 1. SEM photographs of: fine-grained powder of 7 um or less (a); submicropowder (b) 


Vibroscreens were applied for the grain-size classification of diamond powders. . Physical and mechanical 
characteristics of the produced samples were tested through standard methods. VLTE-500 electronic fourth-class labor- 
atory balance was used for weighing. Density (p) was determined by MK 0-25 mm micrometer according to GOST 
6507-78. 

Research Results. Hardness was measured on Equotip 3 device, Proceq, according to the corresponding meth- 
od of the GOST. The proportionality limit in compression and the elasticity modulus were determined according to 
GOST 25.503-97. The test results are shown in Table 2. 




















Table 2 
NDP impact on physicomechanical properties of M2-0lalloy 
Proportionality limit ; Compressive 
Actual densit 
ND, % Porosity, % Hardness, HB in Compression Opp, one oo modulus Ec, 
MPa Pee MPa 
weg 2 Z 2 Z Z 
QOze-ia lex /l/F<lalisi/e/F#l/alale!/Fsl|Blale|e]é6 a 
8 He a é | Z ate | Z a é | Z a é | Z Z 
0 41.0} 41.0) 41.0) 41.0} 41 | 41 | 41 | 41 | 9.7 | 9.7 | 9.7 | 9.7 | 7560 | 7560 | 7560 | 7560 5735.94 
1 29.0 | 30.0 | 30.0) 28.0} 43 | 47 | 48 | 48 | 10.5] 11.5} 11.0) 11.5 | 7630 | 7620 | 7610 | 7670 6771.03 
2 27.0 | 28.0 | 28.0) 26.0} 43 | 43 | 51 | 53 | 10.5] 11.5] 12.0) 12.5 | 7710 | 7700 | 7760 | 7750 6953.50 
3 25.0 | 28.0 | 27.0 | 26.0} 47 | 43 | 51 | 53 | 10.2} 10.5} 12.1) 12.2 | 7730 | 7720 | 7770 | 7750 6580.27 






























































To determine the elasticity modulus and proportionality limit, the samples were compressed on presses to a 
relative deformation of 15-16%, at a loading rate of 0.2 kN/s. 
The porosity (P) was calculated through the actual and theoretical density using the formula: 
Il = (1 — p/p,) x 100%, (1) 


where p, is the theoretical (calculated) density of nonporous material; p is the actual density of the sample. 
The theoretical density was obtained by the formula: 
Pr = 100/(Cy/pi + Co/p2 + C3/p3 + Ca/pa), (2) 
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where C}, Cy, C; and Cy are concentrations of copper, tin, NDP and UNDN in the powder mixtures according to their 
density P), P2, P3 and py. 

The densities taken into account are: 

- copper: 8.96 x 103 kg/m’, 

- tin: 7.28 x 103 kg/m’, 

- diamond: 3,5 x 103 kg/m’, 

- UDND: 3.1 x 103 kg/m’. 

The volume of samples was calculated by the formula: 

V=(m,—m)/ Px, (3) 
where V is pressing volume; 7m is mass in air; my is mass in water; p,, is liquid density. 
In the course of the study, it turned out that with a decrease in the filler size, there was an improvement in the 
physicomechanical properties of the bond modified with diamond powder. 

The elasticity modulus values were determined for the samples with the addition of the UDND and a pure 
binder. The best indicators are shown by the samples with the filler from UDND. At the same time, the physicomechan- 
ical properties under consideration deteriorate if the volume of the added UDND particles exceeds 2%. 

The metallographic studies of the samples allowed us to establish how diamond particles affect the matrix 
structure. Fig. 2 shows the images of microstructures of the deformed samples. 





Fig. 2. Images of ground surface of deformed sample with addition of 2% of diamond powder particles under x 1000 magnification. 


The pictures show clearly visible narrow and branched microstructure objects. Hypothetically, these are 
boundaries between the grains or microcracks formed under the deformation. 

Inside the grains, there are also point microscopic objects that form a disperse substructure. Compared to the 
boundary-distributed point objects, their density is much less, but significantly more than in the original matrix which 
has no diamond fillers. 

Based on the results of metallographic studies, it can be argued that the strengthening of the matrix material 
has two mechanisms — dispersion and grain-boundary ones. 

If we are talking about a dispersion mechanism, the volume of NDP introduced into the matrix material can be 
calculated using Orowan equation [13]: 

Gy = ko Ins, (4) 
where A is the nearest distance between the particles; G is matrix shear modulus; b is Burgers vector; ko is coefficient 
characterizing the pattern of interacting atoms with dislocation. 

The following values are chosen: G = 0.367 x 10° MPa for bronze; b = 2.564 A for copper; ky coefficient is equal 
to 0.85. 

The nearest average distance between the particles, depending on the content and dispersion, is calculated by the 
formula proposed in [14]: 





4 cose) _ | d, (5) 


1,91Ly 


where L;; is weight fraction of the filler; d is the diameter or thickness of the filler particles. 
Table 3 shows the calculations of an average distance between filler particles depending on their volume and size. 
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Table 3 
Nearest average distance between filler particles depending on their volume and size, and hardening according to 
Orowan equation under introduction of diamond particles 

















Particle size 7/5 3/2 —40 UDND 
Particle content, % 1 2 3 1 2 3 1 2 3 1 2 3 
A, um 125.3 | 98.19 | 85.02 | 52.19 | 40.91 | 35.42 | 417.56 | 327.30 | 283.40 | 6.26 | 4.91 | 4.25 
Cv, MPa 0.13 | 0.16 | 0.18 | 0.28 | 0.35 | 0.40 | 0.04 0.05 0.06 | 1.91 | 2.38 | 2.70 












































The calculated data were substituted into Orowan equation, and thus the strengthening was determined through 
introducing the diamond powder particles into the matrix material. 

According to the calculations, the greatest strengthening is provided through introducing the UDND into the 
matrix, which is generally validated by the experimental data. 

When the grain geometry changes due to the agglomeration of filler particles at the interfaces in the material, it 
is advisable to calculate the material properties change according to the theory of grain-boundary strengthening [15, 
16]. 

To determine the quantitative increase in the strength of the material through adding particles of diamond 
powders due to the grain-boundary strengthening, calculations were made using Hall — Petch empirical relationship 
[17]: 

Ao; = kd, *””, (6) 
where k is Hall — Petch coefficient for this material; d, is grain size. 

For calculations, we used the samples showing the greatest increase in strength according to Orowan theory. 
The calculations were performed according to the data obtained from the processing of the surface microstructure imag- 
es through the technique proposed in [18]. 

The Hall — Petch coefficient is applied to copper. According to [19], it is a variable value; it depends on the av- 
erage grain size and varies in the range of 0.01-0.24 MPa x m’”. The calculations show that the greatest strengthening 
is provided by the introduction of ultrafine NDP into the matrix material. In general, this is confirmed by the experi- 
mental data. 

The average grain size is calculated according to the metallographic studies of the sample surface: 


d; = (res) 
where S,o,, is total area of objects; Ny5,, is total number of objects. 
With an average grain size of about 10”' um, the Hall — Petch coefficient is about 0.01 MPa x m!”. 
The calculations based on Hall — Petch ratio indicate an increase in the yield strength of the material with the 
addition of particles of NDP. The yield strength reaches the maximum design value when the content of fillers is 1% 


(Fig. 3). 


; (7) 


15 


Aor, Mila 
n 


0 1 2 3 4 
UDND content, % 


Fig. 3. Yield strength - UDND content dependence 


Discussion and Conclusions. The calculations show that the samples with the addition of NDP have a smaller 
grain size compared to the initial ones. This fact can be explained as follows. The diamond particles, settling at the grain 
boundaries of the material, help to reduce their average size. As a result, the geometry of the boundaries between the 
grains changes; barriers to dislocations are formed; and thus, the potential capability of the material to resist plastic de- 
formation is enhanced. 
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The yield strength increases by about 12—13 MPa which correlates with the calculated data obtained using Or- 
owan theory for dispersion strengthening. If we are talking about polycrystalline material, then, in principle, the yield 
strength increases with decreasing the grain size. Diamond particles added to the matrix volume increase the yield 
strength since they change the geometry of grains reducing their average area and size. 

The calculated data differ from the experimental results since they show an increase in the material strengthen- 
ing pro rata to the number of diamond particles introduced into the volume. It can be assumed that the Orowan strength- 
ening model does not consider the formation of carbon and the agglomeration of diamonds into larger objects in the 
volume of the matrix when increasing the number of diamonds introduced. 


The decrease in the number of pores when adding the UDND particles in the amount of 1—2% can be explained 
by the high sorption properties of the filler. During sintering of compacts obtained through the powder metallurgy, the 
UDND particles absorb oxygen contained in the powder mixture with the formation of CO and CO, reducing gases. 
These gases destroy the oxide film covering the powder mixture particles and prevent oxidation during sintering, there- 
by reducing the total volume of gases in the powder mixture. At the same time, reducing gases accelerate the sintering 
process of the material. The combination of these factors ultimately reduces the residual porosity in the material which 
is validated by the calculated data. If the UDND particle volume reaches 3%, an increase in the carbon content in the 
material occurs. As a result, the filler particles are not fully oxidized, thereby increasing the number of pores in the ma- 
terial. 
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Introduction. The development of computational-and- 
experimental methods for evaluating the force distribution pat- 
tern across the width of the toothed rim and in the circumferen- 
tial direction of the toothing of a large wave gear is considered. 
The study is based on the results of the test tooth tensometry 
using scale modeling of prototype units. The work objective is 
to create a reliable experimental-theoretical model of the teeth 
force interaction in a large wave gearing. Such a solution in- 
volves the transformation of model sample deformations into a 
distributed load between teeth which will eliminate the basic 
uncontrollable nonlinear errors and improve the accuracy of 
estimation of force factors in the toothing area. 

Materials and Methods. An improved power analysis procedure 
of a large wave gearing, optimized by accuracy criteria, is de- 
veloped. The accuracy of the research results is enhanced 
through improving physical and computational models. This 
approach enables to obtain reasonable dependences of the pow- 
er factors distribution in a large wave gearing. 

Research Results. The design shape of the control tooth is sim- 
plified; an invariant profile is introduced over the full width of 
the ring gear. Thus, non-linear distortions of the experimental 
results introduced by a variable tooth shape across the width of 
the ring gear are excluded. In this case, the installation of ten- 
soresistors across the full width of the test tooth is possible. In 
addition, the proposed solution can establish the dependence of 
the teeth deformation across the full width of the ring gear, and 
not only in the extreme areas as suggested by the well-known 
techniques. The development of perfect physical and mathemat- 
ical models enables to increase accuracy of the results of theo- 
retical and experimental studies on power processes in the large 
wave gearing. The scientific-based two-parameter dependences 


of the force distribution in gearing are obtained. 





Beedenue. CtTaTba MocBalleHa pa3padoTKe pacueTHO- 
SKCIHePHMCHTAJIbHBIX MCTOJOB OL[CHKH XapakTepa pacripeye- 
JI@HHA CHJI 10 WHpHHe 3yOuaToro BeHWa HM B OKpy2KHOM 
HallpaBlleHHH 3yOuaTOrO 3al[eMIeHHA KPYyNHO BOJHOBOI 
nepenaun. UccneqoBaHve OCHOBaHO Ha pe3yJIbTaTaX TeH30- 
Me€TPHM KOHTPOJIbHOrO 3yOa C HCIOb30BaHHeEM MaclTaoHo- 
TO MOJeJIMpOBaHHA ONBITHBIX OOpa3yoB. Lem padoTsr — 
co3qaHHe OObCKTHBHOM 9KCIepHMeHTasIbHO-TeopeTHyecKolt 
MOJeIM CHIOBOTO B3aHMOelicTBHA 3yObeB B 3alleIvIeHHu 
KpynHOol BOMHOBOM Mepeqaun. Takoe pelieHue mpeznosara- 
eT MpeoOpa30BaHHe edopMallMii MOJebHEIX OOpa3l[oB B 
paciipeyjeneHHy!lo Harpy3ky MexTy 3yObAMH, 4TO MO3BONT 
YCKIIOYUMTh OCHOBHBI€ HeKOHTpOJIMpyeMble MOorpeliHocTu 
HeJIMHeHHOTO XapakTepa, NOBbICHTb TOYHOCTb OLCHKH CHJIO- 
BBIX (PAKTOPOB B Hose 3al|eWIeHHA. 

Mamepuaaei u memoodvi. PaspadotaHa yrouHeHHad MeTOAH- 
Ka CHJIOBOTO aHaliM3a 3yOuaTOrO 3aleIWIeHHA KPyNHOM BOsI- 
HOBOM Mepeqauu, ONTHMU3HpOBaHHad 0 KpHTepHAM TOUHO- 
cTH. TOYHOCTR pe3yIbTaTOB VCCeqOBaHHA NOBbIMIeHa 3a 
C¥eT COBEPUICHCTBOBaHHA (PH3M4eCKHX MW PaCueTHBIX MOJIe- 
nei. Tako MOzXo MO3BONH NONyYNTh OOOCHOBaHHBIC 
3aBHCHMOCTH paciipelesleHHA CHJIOBLIX PaKTOPOB B 3yOua- 
TOM 3alleIvIeHHH KpyMHOl BOUHOBOM Mepeqayn. 
Pezyiemamel uccaedoeanua. YipowjeHa KOHCTPpyKTHBHaA 
(opMa KOHTpombHOro 3y6a: NO BCeli WHpHHe 3yOuaToro 
BeHI[a BBeJeH HeH3MeHHbIM mpodusb. Takum oOpa30m Hc- 
KJIOUCHbI HeNMHeMHbIC UCKaKeCHHA Pe3yIbTATOB IKCHepH- 
MeCHTOB, BHOCHMBIe TepeMeHHOM (opMbI 3yOa No WupHHe 
3yOuaToro BeHya. B Takom cylyuae BO3MO%KCH MOHTAK TeH- 
30pe3HCTOpoB M0 BCei WIMpHHe KOHTpoubHOrO 3yOa. Kpome 
TOPO, WpeyaraemMoe pellieHve MO3BOJIAeT YCTAHOBHTb 3aBH- 
CHMOCTb JeopMaluu 3yObeB 10 BCei WIMpHHe 3yOuaToro 
BeHIa, a He TOUBKO Ha KpaliHHXx ydacTKax, Kak MpeyiararoT 
W3BeCTHbIe MeTOQMKH. Pa3spaboTka COBepIIeCHHbIX Pu3H4e- 
CKHX HM MaTeMaTHY¥eCKHX MOJesIei MO3BOIMIa MOBbICHTb 
TOYHOCTh pe3yIbTAaTOB TeOpeTHYeCKUX HW IKCIePHMeHTAIIb- 
HBIX HCCI€QOBaHHM CHJIOBbIX IIpOWeccosB B 3yOuaTOM 3al{er- 
JI¢HHH KpynHOl BouHOBOM nepeyauu. TomyyeHbt HayyHo 
oOocHOBaHHbIe JByXMapaMeTpHyeckHe 3aBHCHMOCTH pac- 
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Discussion and Conclusions. Approximation of the involute- 
tooth profile in the trapezoidal profile has simplified evidence 
of identity of the elasticity equations and the boundary condi- 
tions of mathematical models. The results obtained are applica- 
ble in the mathematical simulation of the planar stress state of 
teeth with nonlinear profiles. Comparative evaluation of errors 
introduced by deviations of geometry and dimensions of physi- 
cal models and mathematical analogues supports the experiment 
correctness and the validity of the quantitative data obtained. 
The research results can be used in the improved calculation of 
the design parameters of the gear components in the engineering 


process of large heavily loaded wave reducers. 
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TipeqeueHHA CHI B 3yOUaTOM 3allersIeHHH. 

O6cyorcdenue u 3akiiouenua. ATITIPOKCHMalMA 3BOJIbBCHT- 
HOrO Mpodusa 3yOa B TpallelleHAabHBI Mpopusb yupoctTu- 
Jla OKa3aTesIbCTBO TOKJCCTBCHHOCTH ypaBHeHH yupyro- 
CTH HW rpaHW4HBIx yCNOBui MaTeMaTH4eckux Mozesen. Ilo- 
JIYYeHHBIe pe3yIbTATbI MpHMeCHHMbI TIpH MaTeMaTH4ecKOM 
MOJeIMpOBaHHH TIOCKOHAaIIpsxKeHHOTO COCTOAHHA 3yObeB C 
HeJIMHeHHBIMH TMpodusiamu. CpaBHHTeIbHat OeHKa T0- 
TpelHOcTeH, BHOCHMBIX OTKJIOHCHHAMM TeOMeTpHYeCKHX 
(popM H pa3sMepoB Pu3H4eCKHX MOJeseH HM MaTeMaTH4eCKHX 
aHaJIOrOB, HOATBep2KaeT KOPPCKTHOCTb MOCTAHOBKH 9KCIIC- 
PHMeHTa HW OOOCHOBaHHOCT MOJYYCHHBIX KOJIMYCCTBCHHBIX 
WaHHBIX. Pe3synbTaTbl paOoTb! MOryT ObITb HCTIOb30BaHbI 
IIpH yTOYHCHHOM pacueTe KOHCTPYKTHBHBIX TWapaMeTpoB 
QIEMCHTOB 3YO4aTOTO 3alleIWIeHHA B Mpolecce mpoeKTHpo- 
BaHHA KpYIHBbIX TAKEO Harpy2KCHHBIX BOJIHOBBIX PeYKTO- 
pos. 


Ksrouesble c10Ba: BONHOBad Mepeqaya, 3yO4uaToe 3arlerse- 
Hue, CHIOBOM aHas3, pactipeyesienve Harpy3ku. 
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BOHOBOH Mepegauu / B. H. Crpenpuuxos, A. VW. Bonomuny, 
M. I. Cyxos // Bectuuk JJon. roc. TexH. yH-Ta. — 2019. — 


T. 19, Ne2. — C. 120-129. https://doi.org/10.23947/1992- 
5980-2019-19-2-120-129 


Introduction. The intermediate transformation of the rotational motion into a continuous wave deformation of 
the flexible wheel has changed the established principles and forms of conjugation of the engagement elements of high- 
er kinematic pairs. At this, a small tooth difference in the internal gearing of wave transmission (two teeth) differs from 
the traditional ideas on the conditions of geometric synthesis and functioning of a gear pair. Therefore, from the view- 
point of interference of the teeth, the engagement of the wave transmission is beyond the permissible limits of the exist- 
ence of an involute internal gear with rigid tooth-wheels. The interference of second-type teeth arising in clamping is 
enhanced by the deformation of the flexible wheel from the transmitted load, which limits the carrying capacity of wave 
gear. The negative impact of the scale factor contributes to the occurrence of tooth interference in large wave transmis- 
sions, which causes jamming and slipping of teeth in the gearing. 

The basic concept of wave gear (strain wave gearing, SWG) was patented by K. W. Musser in 1957. The solu- 
tion remains unchanged for the stock-produced designs of limited capacity despite many design improvements. This is 
confirmed by wave gears manufactured by the leading companies in the USA (Harmonic Drive Technologies Inc., 
United Shoe Machinery Corp.), Japan (Harmonic Drive Systems Inc.), Germany (Harmonic Drive AG), and others in- 
cluding licensed products of Harmonic Drive Systems Inc., China [1, 2]. In large wave gearheads, high capacities are 
realized, and torques exceed (0.3 ... 1.5) x10° Nm. Mining, metallurgical equipment and large machines produced by 
Novokramatorsk Machine Building Plant (NKMZ) are completed with such parts. The standard design of a cam wave 
generator with a flexible ball bearing is inoperative and inapplicable in such products of heavy engineering [3]. Large 
wave gears use a disk wave generator. In contrast to the cam generator which simultaneously contacts along the entire 
perimeter of the flexible wheel, the disk wave generator interacts with the flexible wheel in diametrically opposite areas. 
This creates a higher level of freedom of deformation of the flexible wheel and increases possible deviations from the 
specified position. 

Under the action of the wave generator discs and torque, the flexible wheel becomes curved in the form of a 
slightly twisted taper. This causes tooth misalignment over the width of the ring gear relative to the rigid wheel. The 
uniformity of the gap setting in the gearing across the width of the gear rims is disturbed. 

Under heavy-duty service, the deformations of the flexible wheel exceed the gaps in the gearing and create in- 
terference conditions for the teeth of the second kind. With wide toothed crowns, a small module of the teeth under 
heavy loads, the gaps in the local zones of the deformed gearing take on negative values, activate the interference, cause 
jamming and breakthrough. Constructive features of large wave gears are mainly due to the scale. This factor exacer- 
bates the negative processes observed in the gearing of higher kinematic pairs, and reduces the specifications and per- 
formance of heavily loaded wave gearboxes. 

The work objective is to increase the load capacity, the specifications and performance of large wave gear re- 
ducers that meet the technical requirements of heavy engineering. 
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Theoretical and practical developments in the study of large wave transmissions are limited, the results ob- 
tained are contradictory, and they cannot always be repeated. Besides, we note significant construction features, lack of 
geometric similarity in the size series of kinematic and large wave transmissions among others. All this does not allow 
using well-known computational methods in the design of wave transmissions for the purposes of heavy engineering 
[4-6]. Synthesis of gearing is subject to the condition of constancy of a given gear ratio, governed by the Euler — Savary 
formula. Permanent deformation of the flexible wheel does not agree with the fundamental gearing theorem, which con- 
tributed to the development of conflicting methods for the synthesis of wave gearing [7—11]. In some papers, a general 
case of gearing is taken as the basis [12], while others consider a wedge mechanism with a complex relative movement 
of links [13-15]. In some techniques, the geometry and kinematics of wave gearing are simplified [16, 17]. The com- 
bined approach is synthesis of the general theory of gearing with elements of the wedge mechanism [1 ]. 

To satisfy the conditions of the fundamental theorem of gearing in relation to the synthesis of wave transmis- 
sion, E.G. Ginzburg used the kinematic approach [18]. The concept of “angular velocity of a point” introduced by him 
is wrong, the arguments based on it are incorrect. 

In the monograph by N.A. Kovalev [19], complexity of the force and kinematic processes occurring in wave 
transmissions, absence of reliable dependencies that can adequately reflect the impact of external factors on internal 
processes in kinematic pairs are specified. It prevents from developing an effective method for synthesis of wave gear 
that optimizes the basic parameters of the bearing links. According to N.A. Kovalev, with an inextensible middle sur- 
face, only the point of intersection of the tooth axis and the middle surface moves uniformly, which contradicts the con- 
clusions of E.G. Ginzburg [18]. 

Maximum radial elastic displacement (w,) on the middle surface of the flexible wheel is set by the wave gener- 
ator in the absence of torque. The w, offset must exceed half the working height of the tooth. The toothing dimensions 
should be adjusted by w, value, not vice versa [19]. 

There is no consensus on the value of w, parameter. Some authors associate maximum radial deformation 
(Wmax) With the transmission ratio of wave gear [20] and a disk wave generator. To increase efficiency, it is recommend- 
ed to assume W,,,,, = m. E.G. Ginzburg recommends to determine Wma, through m, conditional module, different from m 
module [18]. 

Materials and Methods. To prevent the interference of teeth, the following well-known techniques are used in 
the stock-produced wave gearboxes: 

- tooth correction, 

- widening of the tooth space, 

- angle increase in the original tooth contour. 

This, along with a narrow width of the gear rims and relatively small transmitted torques, in many ways settles 
practical issues of the approximate gearing synthesis in terms of eliminating the tooth interference for the stock- 
produced low-loaded wave gears. However, theoretical studies on these issues are far from complete. So far, there is no 
satisfactory solution regarding the engagement of higher kinematic pairs of large wave gears with the gear rim width of 
100-200 mm, a tooth module of 1.5—3 mm and a disk wave generator. Large wave gearboxes, designed by well-known 
methods and manufactured at NK MZ, turned out to be inoperable due to the tooth interference. Solving such problems 
requires new technical solutions based on reliable results of the theoretical and experimental research. 

A scientific research has been carried out to assess the boundary conditions for the interference phenomena and 
tooth breakthrough in the gearing, as well as to prevent these negative developments. Their results enable to create a 
generalized model of the tooth interaction in the engagement of a large wave transmission. A refined power analysis 
procedure of such toothing is developed. The research results accuracy is enhanced through improving the physical and 
computational models. This allows obtaining reasonable dependences of the distribution of power factors in the gearing 
of a large wave transmission. 

Deformation of stressed structural elements, including teeth, is usually measured using strain gauge methods. 
However, it is impossible to fix strain gages on the wave transmission teeth with a small module. According to the 
known methods [20, 21], two parallel slots are made in the rim of the hard wheel removing two adjacent teeth (Fig. 1). 
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Fig. 1. Localized tooth elements of rigid wheel tilt gear of MP-600AC mobile mixer 


Across the wheel width, the slots can be end-to-end or blind, made at both ends of the rim so that the middle 
part remains intact. According to the authors of common methods, through deep slots can reduce significantly the wheel 
stiffness, which is unacceptable. In the known papers, an experimental model is used, where the middle part of the teeth 
is not removed and it works alongside with all others. 

The ideas of an unacceptable reduction in wheel stiffness through end-to-end slots [22] are inconclusive, since 
the thickness of the rim of a rigid wheel is not regulated; and it can be taken sufficient under the experimental condi- 
tions that the slots have no real impact on rigidity. In addition, the hard wheel is pressed into the hull design and, when 
evaluating stiffness, is considered together with the gearbox housing. . 

According to the methods [20, 21], the slots of the rigid wheel increase the deformation of the control tooth 
relative to its middle part. The change in the rigidity of the control tooth along the length increases the deformation of 
the selected elements at the ends, distorts the deformed state of the control tooth and does prevents with sufficient accu- 
racy the estimation of the force characteristics of the gearing along the crown width and in the circumferential direction. 
A multiple decrease in the rigidity of the selected tooth elements is not considered by the well-known methods [22], 
which violates the objectivity of the experimental results and complicates their analysis. 

Standard techniques for measuring the tooth deformation of the wave transmission [8] differ little from the 
method of conventional gears [23]. In [22], it is shown that the difference in rigidity of the model, due to blind slots on 
the hard wheel, distorts the deformation of the control tooth under load and is considered when conducting the experi- 
ment and processing of the data obtained. This reduces the accuracy and reliability of the results. 

The strain gauge is not mounted in the middle of the control tooth, so it is impossible to get a full picture of its 
deformation, qualitative and quantitative characteristics of the force distribution in the toothing. Moreover, when con- 
sidering the stepwise change in the rigidity of a physical model by the known methods, the boundary conditions are not 
taken into account, the impact of which is enhanced by the scale factor [18]. 


Consider VZ-1120 wave reducer of the tilt drive of MP-600AC mobile mixer with the capacity of 600 tons of 
molten metal. Assume that load in the gearing of the wave transmission is applied in the middle of the tooth height, 
because there is ~ 20+25% of teeth in the internal gearing of one wave. 

Let us determine a possible error of allowance in the experiment. Consider a rigid annulus internal gear: 


— modulus of teeth m = 1.5 mm; 
— number of teeth of the rigid wheel Z, = 762; 
— number of cutter teeth Z, = 68; 
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— wheel pitch circle diameter d= 1143 mm; 
— wheel tooth depth H = 2.985 mm; 
— wheel outside diameter D = 1155.12 mm; 
— root diameter of wheel teeth D, = 1161.09 mm; 
— pressure angle at the tooth point on the pitch circle d, a = 20 °; 
— tool addendum modification coefficient x = +4.953 mm; 
— face width b = 100 mm; 
— normal rack tooth profile as per GOST 13755-81 standard; 
— accuracy degree 7D (GOST 1643-81). 
The relative error of R assumption concerning the conditional application of normal force in the middle of the 
tooth depth: 
i 5000 — ng aay 
D+D, 1155.12 +1161.09 
The accepted assumption has no real impact on the transmission of forces in gearing, since it changes slightly 
the diameter of the loading appliance. For approximation of the involute tooth profile in a straight line shape, we deter- 
mine the thickness of the tooth (S,) of the rigid wheel along an arbitrary circle of radius (7,): 


R 








COS Oy 7 


. . r 
S,=m +A, +Z) (inva, —invdg ), cos a, = cul dg, 


COS 0 x 


where a, is pressure angle at the tooth point, where radius is r, vector; A, is ratio of the tooth thickness on any wheel 
circumference in the normal section under toothing by the cutter of medium wear: 


Ay =(Z, — Z4) (inv dg — inv ac). (1) 
The relative error is defined as the modulus of ratio of the absolute error to the approximation value. In this 


case, maximum possible absolute error is equal to oe . The approximate value is the radius of the middle tooth circle of 


the rigid wheel: (D+D, ) /4 . The angle cosine of the machine tool a, in the expression (1) when cutting a rigid wheel 


by the cutter of medium wear is determined considering the addendum modification coefficient of the cutter: 














Zn —Z4 
COSAC7 = COSAg aes 
The tooth thickness values (Sy) along the arc of an arbitrary radius (r,) of the rigid wheel are given in Table 1. 
Table 1 
Arc tooth thickness Sy of rigid wheel in different sections with radius ry, mm 
ry 579.239 579.426 579.612 579.799 579.985 580.172 580.358 580.545 
Sy 2.613 2.765 2.918 3.071 3.223 3.377 3.530 3.685 





























Arc tooth thickness differs from the chordal tooth thickness at the given parameters of a rigid wheel by 5-10" 
> mm, so with sufficient accuracy, we can equate the arc tooth thickness with r, radius and the corresponding chordal 
tooth thickness. 

We determine the tooth thickness S,, approximated in a straight profile at the radius of ry = 579.799 mm. It is 
equal to the half-sum of the extreme tooth thickness values along the radii 7, = 579.053 mm and 7, = 580.545 mm: S,, = 
3.074 mm. 

The tooth thickness of the involute profile of the rigid wheel on the circle of radius r, = 579.799 mm is equal to 
S, = 3.071 mm. Absolute error of the approximation of the tooth profile of a rigid wheel in a straight-line profile: 

A, = 8, — S,, =— 0.003 mm. 
Relative error of the straight-line approximation of the involute tooth profile of a rigid wheel: 


Sy = Dae 


R= -100 =0.1%. 











ax 
The approximation of the involute tooth profile of a rigid wheel with a straight-line profile does not introduce a detecta- 
ble error in the design scheme. The equation of the straight-line profile of the tooth model in the XOY coordinates is 


shown in Fig. 2. 
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I =15.947 





Fig. 2. Design model of control tooth of rigid wheels approximated by straight-lineprofile 
In Fig. 2, 1 and 2 indices indicate membership of the coordinate functions X and Y to the extreme | and 2points 
of the straight-line tooth profile. XY and Y values are inserted in the equation (2) at 1 and 2 points as shown in Fig. 2, and 
the equation of the straight tooth profile with linear dimensions in millimeters is written: 


Y =7.767—0.410 X. 


Thickness and rigidity of the approximated straight tooth of a rigid wheel in an arbitrary section, mm: 





Soe = 2Y =15.534—0.820X (3) 
3 
EI, =2.1-10° =, (4) 


Specific rigidity of the approximated straight tooth of a rigid wheel in an arbitrary cross section is determined 
by substituting the values (3) at /= 17.5 mm into the expression (4): 


EI,, =0.169-10°.(18.943—X). 
Deflection of the elastic line of the control tooth model under N; force is determined at_X = /,,= 15.947 mm: 
Y, = 0.819-10* Ny [mu]. (5) 
Tension stresses o* and compression o~ of the control tooth model in the area of installation of strain gages at 


the distance of a, = 7.2275 mm from the cantilever fastening (pinching) of the model: 


My Xs, 
gio HSU (ie 


. Me. ze H 
o =-—4+4-%-0,1983N,,—, (7) 
Wr #F mm? 


Here, M,, is bending moment in the area of installation of strain gauges; M, = M,— Y,-a,; W,is section modulus of the 


linear tooth model: 


2 
W, = ~~ = 39.61 mm’. 


Using the ratio 
o = Es, 
where ¢« is the relative deformation of the strain gauge, as well as the formulas (5—7), a relationship is established 


between the relative deformation of the strain gauges on the stretched and compressed surfaces of the control tooth 
model and the deflection of the elastic line under the force: 
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Y, =5013.33-e, Y> =472.07-e . (8) 

Research Results. Theoretical and experimental studies, practical developments are mainly associated with 
the development, production and operation of relatively small wave gearboxes [25-32]. The accumulated experience 
cannot be used without considering the scale factor for large wave gears with torques of 5-10° Nm and more. 

The flexible link enhances the negative effect of the scale factor through disturbing the regular operation of the 
gearing of higher kinematic pairs. Under heavy-duty service, interference of the second kind can cause jamming and 
breakthrough of the teeth. Interference also activates force processes in the wave generator area increasing energy losses 

The fragment of the oscillograph pattern corresponds to the turnover of 2m wave generator (Fig. 3). 


zZ 





Fig. 3. Oscillogram of rigid wheel tooth deformation under load M, = 3- 10° Nm and rotation frequency @, = 1.667 st; 1, 2,3,4,5 
are deformation oscillograms of localized control tooth elements 


The number of tooth pairs in Zs; gearing is determined by the duration of the engagement of the control tooth 
under a complete rotation of the wave generator (a'+ a"): 


a't+a" 
Zy = , Zp 100%. (9) 





The number of simultaneously engaging tooth pairs and the load shape across the width of the ring gear depend 
on the design parameters of the wave transmission, which determine the conditions of toothing, and represent some 
torque functions. 

The technique of power analysis of gearing for large wave transmission is developed. Two advantages differ it 
from the known solutions. 

1) It enables to eliminate the distortion of the deformed state of the control tooth under load due to the stiffness 
constancy lengthwise. For this purpose, through parallel slots are made on a rigid wheel symmetrically to the axis of the 
control tooth. 

2) It becomes possible to significantly improve the accuracy of the evaluation of gearing power characteristics 
across the width of the crown and in the circumferential direction due to the rigidity reduction of the control tooth. For 
this purpose, physical models of a control tooth and a scaled tooth of a rigid wheel are developed. 

Discussion and Conclusions. The control tooth shape is taken invariant across the full width of the ring gear 
(this differs the presented approach from the techniques [20—22]). Hence, the distortion of the experimental results in- 
troduced by the variable width of the tooth shape is eliminated. 

The solution provides the installation of strain gauges across the full width of the control tooth, which estab- 
lishes a complete dependence of the tooth deformation across the full width of the ring gear, and not just in the outer- 
most sections, as the well-known technique suggests. 

The updating of physical and mathematical models made it possible to improve the accuracy of the theoretical 
and experimental results of force processes in the gearing of a large wave transmission [8—10]. The scientific-based 


two-parameter dependences of the force distribution in gearing are obtained. 
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Mathematical temperature simulation in tool-to-work contact zone during metal turning” 


E. V. Bordatchev', V. P. Lapshin” 


'Don State Technical University, Rostov-on-Don, Russian Federation 
? National Research Council Canada, Ottawa, Canada 


Martemartnueckoe MOjesIMpoBaHne TeMiepatypbl B 30He KOHTAKTAa HHCTpyMeHTa HU 131es1HA pu TOKapHol 


cred 
oOpadorkKe MeTaJLI0B 


E. B. Bopaues’, B. 1. lam 


' HarmonasbHBlit Hay4HO-HccileqoBaTesbCKuli copeT Kanai, OTtapa, Kanaya 


? JloucKoii rocyzapceTBeHHEIi TexHM4ecKHit yHuBepcuter, Pocros-Ha-Jouy, Poccuiickan Penepauna 


Introduction. Two factors of metal turning are compared: the 
dissipated temperature and the power of irreversible transfor- 
mations in the material of the product and the tool. The paper is 
devoted to the issues of mathematical modeling of their link. 
Materials and Methods. The mathematical apparatus is based 
on the modification of the Volterra equation which involves the 
use of double integral. It shows how the thermal energy re- 
leased earlier during cutting affects the current state of tempera- 
ture in the tool-to-work contact zone. In addition to the pro- 
posed new basic mathematical model, the processing effect of 
the observed data on the power of irreversible transformations 
and the measured temperature in the tool-to-work contact zone 
under metal turning are used. The experiments were carried out 
on 1K625 machine and STD.201-1 stand. A specialized soft- 
ware tool for processing information arrays describing the pro- 
cesses occurring during cutting (reaction forces, tool vibrations 
and power of irreversible transformations) was created in the 
Matlab package. The same tool has performed the temperature 
calculation in the tool-to-work contact zone. 

Research Results. The procedure of parametric identification of 
the proposed basic mathematical model is carried out. The re- 
sulting model showed a high degree of proximity of the exper- 
imental data on the temperature in the cutting zone and the sim- 
ulated level of thermal energy; but in the initial section of the 
measurable temperature dependence, the results of these two 
approaches are in rather poor agreement. This can be explained 
by an error of the experimental temperature measurement based 
on the estimate of the thermoelectromotive force (thermal 
EMF) output which is generated as a result of the dynamic 


thermocouple formation in the tool-to-work contact zone. 


Beedenue. ConoctaBienbl 1Ba dakTopa ToKapHOH oOpadboTKH 
MeTaJUIOB pe3aHHeM: pacceHBaemMasd TeMIlepaTypa HM MOLIHOCTB 
HeoOpaTHMBbIxX TpeoOpa30BaHHli B MaTepHasie W3eIMa HM HH- 
cTpyMeHta. CtTaTba NocBallleHa BolpocaM MaTemMaTH4yecKoro 
MOJ{eJIMpOBaHHa HX CBA3H. 

Mamepuanei u memooovi. Matematuueckuli almliapaT OCHOBaH Ha 
MoyMpuKalluH ypaBHeHua Bompteppsi, KOTOpad MpesmouaraetT 
VCIOJb30BaHHe AByKpaTHOro HHTerpasa. OH MOKa3bIBaeT, Kak 
BbIeeHHad paHee Ip pe3saHHM TelIoBad 9Heprua BIIMAeCT Ha 
TeKylllee COCTOAHMe TeMIepaTypbl B 30H KOHTaKTa HHCTpy- 
MeHTa Cc OOpaOaTEIBaemol JeTambro. B paboTe NOMHMO Tpe- 
JIO%KCHHOM HOBO Oa30B0H MaTeMaTHYeCKOM MOE HCIOIb30- 
BaHbl pe3yIbTaTbl OOpaOoTKH IKCICPpHMCHTAJIbHbIX JJaHHBbIX O 
MOIIHOCTH HeoOpaTHMbIxX Tpeobpa30BaHuli HW OO H3MepsxeMoli 
Temilepatype B 30H€ KOHTaKTa HHCTpyMeHTa C JeTasIbIO pu 
TOKapHoH oOpaboTKe MeTasIa. OKCHepHMeHTHI MpOBOTHWIHCh 
Ha cTaHKe 1K625 u crenye STD.201-1. B naxere Matlab Opa 
co3aHa ciielMasIu3HpoBaHHasd MporpamMa JyId OOpaOoTKH Mac- 
CHBOB HHopMallMu, OMMCHIBaFONIMX MIpolecchl, mpoTeKaronHe 
IIpH pe3aHun (CHIBI peakIMM, BHOpallMM MHCTpyMeHTa HM MOII- 
HOCTb HeoOpaTHMBIX TpeoOpa30BaHuii). B stot »%e MporpamMe 
BbIMOJIHeCH pacdeT TeMIepaTypbI B 30He KOHTAaKTa MHCTPyMeHTa 
Mi JleTasmn. 

Pezyabmamoi uccaedoeanua. Upopeyena mpouenypa mapamet- 
puyeckol HAeHTHUKalMN Mpelaraemoli 6a30B0i MaTeMaTH- 
yeckol Moyemu. IlomyaeHHad MOJeJIb MOKa3asia BbICOKYyIO CTe- 
TICHb OJIM30CTH 9KCIePHMeHTaJIbHBIX JaHHbIX O COCTOAHHM 
TeMIlepaTypbl B 30He pe3aHHaA HM MOJeMpyeMoro YpoOBHA TelI- 
JIOBOM 9Heprvu, OHaKO Ha HayasIbHOM y4acTKe H3Mepsemoli 
TeMMepaTypHOM 3aBHCHMOCTH pe3yIbTaTbI ABYX ITHX MOAXOTOB 
He BIONHe cormacyroTca. ITO MOXKHO OObACHHTb MorpelHo- 
CTbIO 9KCIICPHMCHTaJIBHOTO H3MepeHHA TeMIIepaTypbI Ha OCHO- 
Be OIe€HKH BbIpabaTbIBAeMOM TePMOIICKTPpOABMKYyMIeH CHIBI 
(tepMoD JIC), koTopat dopmupyetcsa B pe3sybTaTe OOpa3z0Ba- 
HHA ECTeCTBCHHOM TepMonapsI B 30He KOHTaKTa HHCTpyMeHTa 
nerTasn. 


* The research is done with the financial support from RFFI (grant no. 19-08-00022 A). 
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Discussion and Conclusions. The proposed mathematical model 
enables to adequately describe the conversion of the mechanical 
component of the cutting energy into the thermal component 
through the indicator of the total output power of the mechani- 


cal interaction in the cutting zone for all the processing time. 


Keywords: tuming, cutting, power of irreversible transfor- 
mations, metal working, thermal conductivity of metals, Volter- 
ra integral equation, thermoelectromotive force (thermal EMF). 
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O6cyatcdenue U 3AKIIOYEHUA. IIpeanaraemaa MaTeMaTH4ecKaAv 
MOJICJIb MOSBOJIACT ATCKBATHO ONMChIBATh Mpolecc Tipeo6pa3o- 
BaHHa MexaHHyeckon COCTAaBJIAIOINIeH SHeprHu pesaHhs B TeII- 
JIOBY!O COCTABJIAFOML YO Wepe3 NOKa3zaTejib cyMMapHolt BBITCJIA- 
emoli MOIJHOCTH Me€xaHiNdecKoro B3aMMOJICHCTBHA B 30HE pes3a- 


HHA 3a BCe BpeMa OOpadoTKH. 


Ksouepbie copa: TOUeHHe, pesaHHe, MOWJHOCTh HeoOpatTu- 
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Introduction. Thermal processes occurring under metal processing on the machine tools (both in metals and in 
other environments) are hard to describe. Their analysis is based on the Fourier heat condition equation [1]. In cases 
related to metal cutting, you need to have an idea about the temperature in the tool-to-work contact zone [2-5]. At that, 
there is no need to control the heat propagation beyond this zone since it does not affect the processing conditions. 
Similar reasoning can be applied to another widely used method for describing thermal processes, which considers the 
formation of heat source-drains in a workpiece [6]. Such approaches do not allow using the cutting process dynamics in 
the description of emerging thermal phenomena in the contact zones and developing a system of equations to represent 
the conversion of mechanical energy into thermal energy. 

During operations with metals on the machine tools, irreversible transformations in workpieces and tools cause 
the heat generation, which, according to the second law of thermodynamics, is dissipated in space. Thermal processes in 
metals are more inertial than the treatment process dynamics, therefore dissipation lasts long enough. In the case in 
hand, the heat generation itself is associated with the power of irreversible transformations [11], which depends on the 
magnitude of the reaction forces to the tool forming movements [2—5] in the tool-to-work contact zone. At the same 
time, heat generation is the primary mechanism of dissipation of the released energy of material conversion [12]. 

The cutting process dynamics and the issues of its simulation are directly related to the representation of forces 
and reactions in the coordinates of the process state. In this case, as a rule, only the mechanical component is consid- 
ered, as can be seen from foreign [13-18] and national [19-21] papers. Thus, it is advisable to introduce coordinates 
describing the temperature variation in the tool-to-work contact zone into the simulation of dynamics of the metal pro- 
cessing processes through consideration of the irreversible transformations. This approach will increase adequacy of the 
used mathematical apparatus. 

Materials and Methods 

Mathematical Model Validation. The forces preventing the tool from penetrating into the workpiece material 
have a complex spatial arrangement in the coordinate system associated with the tool strain axes. In the papers on the 
analysis of the reaction forces occurring in the cutting zone, it is mostly common to resolve such a reaction into compo- 
nents [4—6, 10, 19-21]. In the case under consideration, the following relations between the components of the reaction 


forces will be valid [5]: 
R=[P2+P2+P2, (1) 


where R is general force response from the cutting process to the forming movements of the tool; ved is projection of the 
force response on the axis of the main rotational cutting movement in the feed direction; e is projection of force reac- 


tion directed along the radius of the main rotational cutting movement of at the cutting top; P. is projection of the force 
reaction on the axis coinciding with the main movement speed at the tool tip. 
The relationship between P., P, ,P. forces introduced depends on many factors: the tool geometry, the tool wear 


factor, etc. [4]. Thus, in [5], when machining with a sharp cutter with the parameters of y = 15°, 9 = 45° and 1 = 0°, the 
ratio is on average equal to 
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P,, Py, P. = (0,3 — 0,4), (0,4 — 0,5). (2) 
Similar reasoning can be provided on analyzing the speeds of the relative tool movement along the workpiece, 
and we obtain the total speed of such movement: 


V= We +V, +V,7. (3) 


With these arguments in view, we define the power of irreversible transformations: 


n=,/p? +P? +P2 V2 +V,7 $V, (4) 


Suppose that heat generation during cutting in the metal-cutting machines linearly depends on the power of ir- 
reversible transformations: 
T, =krN . (5) 


Here, T - is the value describing the temperature increments in the cutting zone (°C); k,. is the coupling coefficient be- 


tween the power of irreversible transformations in the cutting zone and the magnitude of the increment of thermal ener- 
gy in the tool-to-work contact zone with c°C/Nm dimension. However, the total temperature value at a specific instant 
in time and in a given tool-to-work contact zone will be determined not only by the current value of heat gain. Due to 
the heat dissipation of in the material of the workpiece and tool, the impact of the following should be considered: 

- power of irreversible transformations throughout the previous path, 

- product processing time. 

Imagine processing in the form of a discrete process graph where M(n) is current value of the power of irreversible 
transformations; L(n) is current value of the tool path during machining; f(n) is current value of the processing time. 
According to (5), at each point of the discrete graph of the treatment process, heat increment proportionate to the power 
of non-reversible transformations at this point in space and time will be generated. The current temperature value at n 
point will be determined by both the heat gain at that point and the effect of temperature increments that previously oc- 
curred on the tool path (Z) during processing (¢). With the increase in time and distance passed, the impact of these fac- 
tors on the process of heat generation in the contact zone under study weakens; therefore, an approach based on the use 
of the Volterra operator of the second kind is convenient for the mathematical description [10]. It cannot be applied 
directly because of the complexity of describing the heat propagation in metals; therefore, we will take the multiplica- 
tive criterion for assessing the effect of the previous heat gain on the current value in the form of a double integral as a 
basic model: 

T, =O, + ky [J wp (yD) (n-D Ny. )dyan. (6) 
D 


Here, 7’ , 1s heat value in the tool-to-work contact zone; 0, is ambient temperature; w,(y—ZL) is the kernel which char- 
acterizes the effect of the previously generated power of irreversible transformations along the processing path on the 
current temperature value; w,(1—f) is the kernel which characterizes the impact of the previously generated power of 
irreversible transformations in terms of processing time on the current temperature value; D is the domain of integration 
characterizing the space-time heat dissipation; y has the distance dimension, m; 1 has the processing time dimension, 


s. The kernels of the integral operator themselves are dimensionless gain constants. 
The integral operator represented as a double integral (6), can be reduced to a multiple integral of the following 


form: 
L(t) t 
T, =O, +kp J we(y-L)dyfwj-DN)dn. (7) 
0 0 


The tool path is a function of time because each time value can be assigned the value of the distance passed. 
Therefore, the approach (7) is valid. 
It is convenient to use exponential functions of the following form [10] as kernels of the integral operator pro- 
posed in (7): 
w,(y—-L) = eh) (8) 
w,(n-t) = eral) * 
where a, G2 are scaling parameters of the integral operator to be identified. 
Thus, the integral operator proposed in (7) takes the following form: 


L(t) t 
T, =O, +k } et Day | eM N(n)dn. (9) 
0 0 
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This operator has the solution for the stationary case when the power of irreversible transformations is con- 
stant: N, = M(t). 
kr No 
00-9 





T, =, + (lee ya). (10) 


As can be seen from (10), under a stationary machining process, temperature in the 
krNo 
O09 

We illustrate these arguments by setting @, =25 °C; ky = 0.0026 s°C / Nm; a, =0.03;a, =0.01; N, =20 
Nis (Fig. 1). 





T, =O, + (11) 


Tz(C°) 





Fig. 1. Implementation of stationary temperature increase in the cutting zone 


Fig. lshows that the stationary heat process in the cutting zone increases exponentially and — having reached 
some equilibrium state — does not change anymore. Real processes occurring in machine tools are not stationary [11— 
12]. Based on these considerations, the process shown in Fig. | can be regarded as a movement towards which the tem- 
perature of the tool-to-work contact zone tends in a steady case. At that, the equation (10) is a generator equation for the 
integral operator (9). 

Research Results 

Parametric identification of the model from experimental data and simulation results. For the parametric 
identification of the model, a full-scale experiment was conducted on 1K625 lathe with the installed STD.201-1 stand. 
The stand is designed to study dynamic and thermal processes occurring during metal cutting on lathes in various 
modes. 

The basic objective of the experiment was to study the impact of cutting dynamics on thermal processes in the 
tool-to-work contact zone. The change in cutting dynamics is understood as the cross effect of tool wear and cutting 
process associated with the force reaction occurring in the cutting zone. Measurement of this reaction and temperature 
through thermoelectromotive force (thermal EMF) generated in the tool-to-work contact zone allows evaluating the 
impact of the cutting dynamics on the thermal processes in the specified zone. 

During the experiment, continuous turning of a shaft-type part was performed. At this, all processing parame- 
ters were preserved and the change in the force response from the cutting process was monitored. The conditions and 
the basic results of the experiment are described in detail in [22]. 

The data obtained as a result of a full-scale experiment are of discrete nature; therefore, it is impossible to di- 
rectly apply the model represented by the integral operator (9) to them. For this reason, it is required to modify the basic 
mathematical model that describes the thermal process dynamics in the tool-part contact. To begin with, we leave in 
reasoning only a member of the integral operator depending on the power of irreversible transformations: 
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L(t) t 
TN = kp } FD dy[ OO NAN. (12) 
0 0 


Since the experiment is already carried out, its completion time (f, ) and Lit ,) path value corresponding to it 


are constants. Therefore, it will be valid to factor the exponents with these values outside the integrals: 


L(t,) t, 
‘iad = kpe the Ou eldy fe ™N(n)dn (13) 
0 0 


Assume that there is a certain discrete set of the calculated power values of irreversible transformations ob- 
tained as a result of analyzing discrete experimental information on the cutting process dynamics, that is 


N (t ) =N ID N. a» N. 3> N. Aeteeeeees On this basis, the integral operator (13) takes the form of the following sum: 
L, t L, t 
TN = kpe the Hh TN, i) e*'dy| e"Ndn+ No } eldy |e "dn+ S 
L 
L, = ar so : , (14) 
+N; } eldy[ edn + Ng } ody [edn + ee ] 
L t, L is 


2 2 3 


Assuming that NV, = 0), we get the final solution for the case described by the expression (14): 


—0,L,) Oh, 





eal = kre [N> (e%EG) ae h) e% — eXtah + road 
4A . (15) 


= +N; (e%EG) — p%hG) etal — eh) 4... ] 

The expression (15) is very convenient for processing large arrays of discrete information on the cutting pro- 
cess. In fact, this is a sum of powers of irreversible transformations weighted relative to the final value. 

Consider only one experiment from the series in which there was no wear on the back surface of the tool. In 
the following cases, the tool wear changed, and hence, the process dynamics (from both a mechanical and a thermody- 
namic point of view). Note that the results obtained in this experiment are generalized for all other cases, but the length 
of the paper does not allow them to be considered. 

Fig. 2 shows the process response to the forming movements of the tool along x axis. 

Fx, kg 











-5 


0 s 10 15 20 25, 


Fig. 2. Process response to forming movements of the tool along x axis in the first experiment step 
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As can be seen from Fig. 2, the reaction to the forming tool movements during turning is essentially non- 
stationary, but there is a certain attracting variety around which a perturbed toolpath is generated. Similar to that shown 
in Fig. 2, reactions to the shape-forming tool movements along the remaining axes were obtained in the experiment. In 
addition, information is available on tool accelerations relative to the workpiece. Integration of these data with consid- 
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eration of the known cutting elements allowed us to obtain the values of the machining speeds and the toolpath along 
the workpiece. 

A specialized software tool for processing information arrays to describe the processes occurring during cut- 
ting (reaction forces, tool vibrations and power of irreversible transformations) was created in the Matlab package. The 
same tool has performed the temperature calculation in the tool-to-work contact zone according to (15). The model pa- 
rameters are presented in Table 1. 








Table 1 
Parameters of the identified model 
cy Oly kr e°C/HM ©,°C 
0.00099 0.00078 0.000159 24.8 




















Fig. 3 shows the simulation results of the expression (15) considering the parameters of Table 1, as well as the 
experimentally measured characteristic. 
ies 6 


350 


300 


250 


Simulated 
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0 5 10 15 20 25 
Fig. 3. Comparison of the results obtained in the experiment and calculated from the model 


Comparison of Fig. 3 and Fig. 2 shows that the measured temperature value, when the tool enters the work- 
piece, increases from 25° C (ambient temperature) to 330-340 ° C. According to Fig. 3, first, a faster growth of the 
measured thermal EMF value is observed, and at the end of the process, the thermal EMF is stabilized. 

This is explained by the fact that at STD.201-1 stand, it is not the temperature as such is measured in the con- 
tact zone, but the thermal EMF of the natural thermocouple formed at the tool-workpiece contact. At the beginning of 
the experiment, due to an overshot of the temperature gradients, the thermal EMF shows too large and fast temperature 
rise, and then, as the temperature gradients diverge, it is stabilized. This effect is widely known. For example, in [5], it 
is said: “The method of a natural thermocouple is unreliable since the cutting temperature is stabilized within 2-3 s”. 

In view of dynamics, the temperature variation (see Fig. 3) obtained as a result of simulating the expression 
(15) is close to the processes occurring during cutting (see Fig. 2). The graph shows the continuing growth of tempera- 
ture due to the continuing increase in reaction forces (see Fig. 2). The graph also reflects all the dynamic features of the 
process including low-amplitude power jumps. 

From this we can deduce the following: the proposed model describes the cutting temperature variation more 
accurately than the data of a dynamic thermocouple. Fig. 3 also shows that the simulated characteristic describing the 
change in the cutting zone temperature, in general, coincides qualitatively with the experimental dependence. The ex- 
perimental and simulated characteristics coincide most fully after stabilization of the cutting process, within 15-24 se- 
conds. 
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Discussion and Conclusions. The study results provide an answer to the following question: how to make the 
description of the hereditary nature of the heat transfer and accumulation during metal turning be adequate to the exper- 
imental data. For this purpose, it is required to simulate the temperature in the tool-to-work contact zone on the basis of 
the proposed modified Volterra operator. 

In science terms, the proposed basic mathematical model and the results of its identification on the basis of the 
experiment on real metal processing are of the main interest. Here, there is a high degree of coincidence of the simulat- 
ed values of thermal energy in the tool-to-work contact zone in the period of steady-state processing. 
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Introduction. Microhardness and wear resistance of the alumi- 
num alloy (AK5M7) after its surface treatment through electro- 
spark doping (ESD) using A (Cu) and P (Cu+P) electrodes with 
different applied energy is considered. The work objective is to 
study the relationship of the physicomechanical and strength 
properties of the surface of AK5M7 alloy obtained after ESD 
with copper and copper-phosphorus electrodes. 

Materials and Methods. X-ray diffraction and X-ray fluores- 
cence analysis methods were used to determine the phase com- 
position, coherent-scattering regions (CSR), and the surface 
microstrain. 

Research Results. It is established that with increasing energy, 
the values of microhardness and wear resistance increase when 
using both electrodes. The sample was modified by P electrode 
with the energy of 0.79 J. In this case, the maximum value of 
microhardness increased 5.3 times, wear resistance - 1.6 times. 
It is found that with ESD, new intermetallic phases of Al2Cu 
and Cu3P are formed on the surface of the alloy under study. If 
the value of the energy used is maximum (0.79 J), the concen- 
tration of the Al2Cu phase increases 5 and 9 times with A and P 
electrodes, respectively, and the concentration of the Cu3P 
phase increases more than 4 times when using P electrode. The 
experimental data obtained suggest size reduction of the CSR 
and an increase in the microstrain values for all the Al, Al2Cu, 
and Cu3P phases on the surface. 

Discussion and Conclusions. An increase in energy with the 
ESD involves an increase in the defective factors and a change 
in the phase composition of the newly obtained surfaces. This 


may explain the increase in wear resistance of the surface. 


Keywords: electrospark doping (ESD), coherent-scattering 
region (CSR), microstrain, aluminum alloys, wear rate. 





* The research is done within the frame of the independent R&D. 


Beedenue. CraTba MocBaAlleHa MCCIeqOBaHHIO MHKPOTBepO- 
CTH MW W3HOCOCTOMKOCTH ayIIOMHHMeBoroO criiaBa (AKS5M7) 
TlocsIe OOpaSoTKH ero MOBeEPXHOCTH IIEKTPOHCKPOBbIM JIerH- 
popanuem (OVJI) snexrpogamu A (Cu) u P (Cu + P) c pasnoit 
TIpHIOKeHHOM 9Heprueli. 

Ilemb padoTbr — wccieqoBaHve B3aHMOCBA3H ()H3HKO- 
MeXaHH4eCKHX HW MpOYHOCTHBIX CBOHMCTB MOBeEpXHOCTH craBa 
AKSM7, osyyenHot mocne OWJI menHbIM uM MeqHO- 
cbocdbopHsiM 31eKTposzaMu. 

Mamepuaavi u memoovi. Ucnionb30BaHbl MeTO]bI PeHTTeHO- 
WudpakYHOHHOrO UM peHtrenodyopeceHTHOrO aHau3a IA 
ompeyenenua ca30Boro coctaBa, oOsacTeii KorepeHTHOTO 
paccesuusa (OKP) u MukpoyxedopmaljHn MoBepxHOCcTH. 
Pesyiemambel ucciedoeanua. YcTaHOBIeHO, 4TO TIP HCTIOJIb- 
30BaHHM OOOMX BHOB 3ICKTPOAOB C yBeIMYeHHeM 9Heprun 
BO3PacTaIOT MHKPOTBepAOCTb HM H3HOCOCTOHMKOCTA. Obdpa3erl 
OBL MoaMuMpOBaH 2IeKTponoOM P c dsHeprueli 0,79 Tx. B 
9TOM Cilyyae MaKCHMaJIbHOe 3HaYeHHe MUKPOTBeEpAOCTH yBe- 
M4MI0Ccb B 5,3 pa3a, H3HOCOCTOMKOcTH — B 1,6 pa3a. Ycra- 
HOBJIeHO, 4TO MpH IVJI Ha noBepxHocTH UccHexyemoro cria- 
Ba OOpa3yIoTCA HOBbIe HHTepMeTamIMyecKHe da3b1 Al,Cu u 
Cu3P. Ecim 3HaueHve HcHomb3yeMol 9HeprH MaKCHMaJsIbHO 
(0,79 Jx), To kKonyeHTpayua chaser Al,Cu yBemmunBaetca B 5S 
u 9 pa3 c olexTpogqamu A u P cooTBeTCTBeHHO, a KOHI[eHTpa- 
ua (azb1 Cu3P yBermunpaetca Oosee yem B 4 paza pH HC- 
TIONb30BaHHH 9eKTpoya P. IlomyaeHupie skcTepHMeHTasIb- 
HbI¢ JaHHbIe CBHeTeIbCTBYIOT 006 YMeHbINeCHHH pa3MepoB 
OKP u yBemmueHHu 3HaueHHi MUkKpoepopMayyu Wa BCcex 
MMeIOWHXcaA Ha HoBepxHocTH paz Al, Al,Cu u CuzP. 
O6cyorcdenue u 3akmouenua. YBenm4enve sHeprun upu OVI 
IIPHBOAMT K MOBbILMeHHIO WeseKTHOCTH MH H3MeHeHHIO ()a30- 
BOTO COCTaBa BHOBb MOJIyYCHHBIX NOBeEpXHOCTeli. ITHM MOxX- 
HO OOBACHHTb MOBBIMNeHHe VM3HOCOCTOMKOCTH WaHHOH M0- 


BepXHOCTH. 


KsouesBble copa: 9sIeKTpOHcKpoBoe MerupoBaHue (OVID), 
oOjacTb KorepeHTHOro paccesHua (OKP), mukpoyedopma- 
IMA, AIIOMMHHeBBIe CIJIaBbl, HHTCHCMBHOCTb W3HalIMBaHHs. 


“E-mail: vice-rector@adm.mrsu.tru, srorismael@gmail.com, Kyachkin@mail.ru 


wee S 
PaOota BbINOHeHa B paMKax HHHUMaTHBHOM HUMP. 


Fomin N. E. et al. Effect of electrospark doping on mechanical properties of Al-Si-Cu allos 





For citation: N.E. Fomin, et al. Effect of electrospark doping O6opaszeu Ona YUumupoeanua: Domun, H. E. Biauanne sy1ex- 

on mechanical properties of Al-Si-Cu alloys. Vestnik of DSTU, TPOHCKpOBOrO JIeTHpOBaHHA Ha MexaHHyeckue cBolicTB Al-Si- 

2019, vol. 19, no. 2; pp. 138-142. Cu cnmapop/ H.E. Momuu, UV. X. Xacan, B. M. Kaurkun // 

https://doi.org/10.23947/1992-5980-2019-19-2-138-142 Bectuuk JIon. roc. TexH. yH-Ta. — 2019. — T. 19, Ne 2. — 
C. 138-142. https://doi.org/10.23947/1992-5980-2019-19-2- 
138-142 


Introduction. Various coatings and surface modification methods are used to increase the wear resistance of 
workpieces. One of them is electrospark doping (ESD) of conductive materials (authors: B.R. Lazarenko and N.I. Laza- 
renko [1—4]). 

Al-Si-Cu alloys are widely used in the automobile industry. They are well molded, have all the necessary me- 
chanical properties, and are characterized by a low coefficient of thermal expansion and a high ratio of strength and 
weight. However, the hardness and wear resistance of aluminum alloys are relatively low. This is the reason for the ur- 
gency of the problem of improving these properties through surface hardening [5—10]. The point is that it is required to 
create hardened layers with improved mechanical-and-physical properties. In the present paper, the results of ESD by 
various electrodes at different energies to harden the surface of AKSM7 alloy are shown. 

Materials and Methods. Samples of AK5M7aluminum alloy with dimensions of 15x15x4 mm were made for 
research. Their surfaces were processed through the ESD method with copper and copper-phosphorus electrodes. The 
elemental composition of the alloy and electrodes was determined through X-ray fluorescence analysis on an ARL Per- 
form’ X 4200 instrument. The results are shown in Table 1. 














Table | 
Elemental composition of alloy and electrodes 
Al, % Cu, % Si, % P,% 
AKSM7alloy 88 7 5 
P-electrode 92.9 7.1 
A-electrode 0.1 99.9 























Electrospark doping was carried out on ALIER-31 installation (SCINTI, Moldova), which allows performing 
the ESD in a wide energy range. 

X-ray phase analysis was performed to determine the phase composition of coatings applied to aluminum alloy 
through the ESD technique. 

The PANalytical EMPYREAN X-ray diffractometer with CuKa radiation equipped with a nickel filter with an 
automatic primary beam divergence slit was used for the survey. The radiograph was decoded using the HighScorePlus 
program. We studied: 

- quantitative content of phases detected on the surface; 
- the average dimension of coherent-scattering regions (CSR); 
- values of microdeformations after ESD of samples. 

The microhardness of the modified layer was measured on the surface obtained after ESD. The research complex 
based on HVS-1000 microhardness meter and a digital video camera were used. The measurements were carried out 
according to GOST 9450-76 with the Vickers pyramid indenter (the load was P = 25 g). At least ten measurements 
were performed for each sample. The base alloy and the surfaces of its samples after the ESD treatment were subjected 
to tribological tests (a lubricant was preliminarily applied on the surfaces of the samples). 

Research Results. The results obtained are presented in Table 2 and in Fig. 1-4. 


























Table 2 
Phase concentration of base alloy and surfaces after ESD depending on energy 
Phase concentration 
Energy, J AKSM7alloy A-electrode P-electrode 
Al ALCu Si Al ALCu Si Al ALCu Si Cu3P 
0 92 2 4 
0.07 97 2 1 94 3 2 1 
0.2 96 3 1 70 22 5 3 
0.39 94 5 1 65 25 5 5 
0.79 89 10 1 59 28 5 8 









































Machine building and machine science 


— 
1eS) 
\o 


http://vestnik.donstu.ru 


140 


Vestnik of Don State Technical University. 2019. Vol. 19, no. 2, pp. 138-142. ISSN 1992-5980 eISSN 1992-6006 





Table 2 shows phase concentrations found on the surfaces of the base alloy and samples after the ESD at vari- 
ous energies. 

The base alloy has three phases: Al, Al,Cu and Si. With increasing the ESD energy, an increase in the Al,Cu 
concentration is observed when using both electrodes, and a decrease in the Al phase, while the Si phase remains con- 
stant for the A- electrode at all energies. In this case, Si concentration increases 2.5 times with the energy rise for the P- 
electrode. In the samples treated with the P-electrode, the Cu;P phase appears. Its concentration increases 8 times with 
the pulse energy rise to 0.79 J. 


——P (Al) = =A (Al) 
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Fig. 1. Dependence of CSR Al-phase size on energy for A and P electrodes 
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Fig. 2. Dependence of Al-phase microstrain value on energy for A and P electrodes 


Fig. 1 and 2 show that when using both electrodes, the CSR decreases with the ESD energy rise, and the mi- 
crostrain grows. The coherent-scattering regions are practically faultless surface areas. Therefore, an increase in CSR 
indicates improvement in the microstructure of the modified layer, and reduction of its defect level. However, it is 
known that a change in the structure imperfection of metallic materials causes a change in their mechanical-and- 
physical properties [11, 12]. The microstrain value indicates the degree of the crystal lattice distortion, and thus, the 
structure imperfection ratio can be defined by its value [13, 14]. Similar CSR dependences and microstrain were ob- 
tained for the Al,Cu and Cu3P phases after the ESD of AK5M7 alloy surface with A and P electrodes. It is found that 
the patterns of changing the CSR and microstrain are the same as for the Al-based phase. Therefore, the data in Fig. 1, 2 
suggest that with increasing energy, the structure of the resulting surface is distorted more heavily. 

Fig. 3 shows the measurement results of the surface microhardness after the ESD by the copper and copper- 
phosphorus electrodes. 
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Fig. 3. dependence of microhardness on energy for A and P electrodes 
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It is seen that the microhardness increases with the energy rise, and for A-electrode, at the peak value of ener- 
gy, nearly 3-fold increase is observed, for P-electrode — 5-fold. 

It is known that during the ESD, melting of materials of the anode and the sample surface occurs, and then rap- 
id cooling (quenching). As a result, new phases (in our case, Cu;3P), dislocations, point defects, high degree of surface 
stress, fine-grained structure, etc., can form on the surface. Accordingly, it can be assumed that with increasing energy, 
the microhardness grows due to the above-mentioned distortions arisen during the ESD process (see Fig. 3). 

Fig. 4 shows the measurement results of the wear rate of AK5M7 alloy and the sample surfaces after the ESD 
with copper and copper-phosphorus electrodes at the energy of 0.79 J. 
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Fig. 4. Wear resistance rate of base AK5M7 alloy and its surface after coating at energy of 0.79 J 


Fig. 4 shows that the sample surfaces treated by copper-phosphorus and copper electrodes are worn out much 
less than the surface of AKS5M7alloy. 

Discussion and Conclusions. The results obtained validate the efficiency of the ESD of AK5M7 alloy surface 
by various electrodes and with different energies. Surfaces obtained by this method wear out less. When using a copper 
electrode, the wear is reduced by 1.37 times, when using a copper-phosphorus electrode — by 1.57 times. 

Upon reaching the minimum value of wear after the ESD, the surfaces have the maximum values of micro- 
hardness and microstrain, as well as the lowest CSR values. At the same time, there is a significant increase in the con- 
centration of intermetallic phases formed after the ESD: Al,Cu — from 5 to 9 times with copper and copper-phosphorus 
electrodes, respectively; Cu3;P — 8 times with a copper-phosphorus electrode. The reason for the increase in wear re- 
sistance of surfaces after the ESD can be various distortions and defects formed after melting, as well as the formation 
of new intermetallic phases. 
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Ucc1eqoBanne 3ABHCHMOCTH OT 4aCcTOThI KOHCTAHT NOJIApPH3OBaHHOl MWbhe30KepaMuKu 
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*}Oxnerit cbeyepanbuerit yaupepcuter, r. Pocros-na-Jlony, Poccnticxas Dexeparna 


Introduction. It is acknowledged that electroelastic modules 
do not depend on the amplitude and frequency of oscilla- 
tions. This approach is reflected in the Russian and foreign 
standards for determining the complete set of electro-elastic 
piezoceramics modules. For example, to determine d3, pie- 
zo-module of a disc-shaped sample, it is required to take 
measurements in three frequency domains: in the first and 
second resonances, in the antiresonance region, and at fre- 
quencies much below | kHz. Accordingly, it is assumed that 
when determining d3,, the modules of the ceramic under 
study in the frequency range from 1 KHz to the second reso- 
nance are independent of frequency. The work objective is to 
study the frequency dependence of electro-elastic ceramic 
modules. In this case, a disc-shaped sample from LZT (lead 
zirconate titanate) is used. 

Materials and Methods. Techniques of setting and solving 
problems of the stationary electroelasticity and sections of 
the electrical engineering basics are applied. To implement 
the finite element method, the perturbation technique and the 
ANSYS application package are used. The experimental 
results are processed in the MATLAB environment. 
Research Results. For the LZT piezoelectric ceramics, the 
frequency dependences of various modules (piezoelectric 
d3,, dielectric ©, and elastic modules of compliance - 
S11, 512,513) were investigated. Radial oscillations of a disc- 
shaped sample with electrodes on the ends were considered. 
The sample thickness was 1 mm, the diameter was 40 mm, 
and the oscillation range was up to 700 KHz. First, the fre- 
quency dependence was studied for the elastic ceramic mod- 
ules from the determination of ten resonance frequencies. 
Then, the frequency dependence of d3,and £43 modules 
was determined from the measured values of the sample 
conductivity. For this purpose, we used the expression for the 
electrical conductivity obtained from the solution of the radi- 
al oscillations of the disc considering its thickness. 


* The research is done within the frame of the independent R&D. 


“E-mail: epohrbats@gmail.com, igorro@rambler.ru, kruglov@sfedu.ru 
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PaOota BBINOJIHeHa B paMKax HHHUMaTHBHOH HUMP. 


Beedenue. IIpuHato cautTaTb, 4YTO BIeKTpoympyrve MoyyIH 
He 3aBHCAT OT aMIVINTYIbI HM 4acTOTbI KONeOaHHH. ITO NoD- 
XO], OTpaxKeH B poccHiickuHx UH 3apyOexKHBIX CTaHapTax [yA 
ompezereHua MomHoro HaOopa 9eKTpoympyrux Mosyueli 
Tibe30KepamMuKkH. Hanpumep, Wa onpeyeseHHA Mbe3oMOTyIA 
d3, oOpa3ia B opMe WMcka HeOOxOAMMO TIpoBecTH H3Mepe- 
HHA B TpeX YaCTOTHBIX OOMacTAX: B OOAaCTH MepBoro UM BTO- 
poro pe30HaHcos, B OOacTH aHTHpe30HaHca HW Ha yacTOTax 
MHoro Hwxe | KI'y. B cBa3M c 9THM Mpesnomaraetca, 4TO 
IIpH ompeyeneHuu dz, MOAyIM uccleayemMol KepaMHkH B 
qMama30He uactoTt oT | KI'y 40 Broporo pe30HaHca He 3aBH- 
CAT OT UaCTOTHI. 

Lembro WaHHO paOoT! ABIIdeTCA UccieqOBaHve 3aBMCHMO- 
CTH OT YaCTOTHI 9ICKTPOYUpyrux Mose KepamuKu. IIpu 
3TOM Hcnomb3yeTca OOpa3ell B PopMe Wucka u3 LTC (ump- 
KOHaTa-THTaHaTa CBHHIA€). 

Mamepuaaoi u memoooi. Uctiob30BaHbl MeTOAbI MOCTAaHOB- 
KH M pelleHua 3ajay cTaljMoHapHOM 9JIeKTpoympyrocTu u 
pa3eubl TeOpeTH4eCKHX OCHOB DIeKTpOTeXHUKH. JIA peasu- 
3al[MH MeTOJla KOHCYHBIX JJICEMCHTOB IIPHMCHEHEI MCTOJ, BO3- 
MYIeCHHH MW WakeT IpHKayqHbIxX nporpamM ANSYS. Pe3yub- 
TaTbI 9KCIepHMeHTOB OOpaboTaHBI B cpeqe MATLAB. 
Pezyaomamoi uccredoeanua. Jina tbe3s0KepamuKn LTC uc- 
CJI€¢OBaHbl 3ABMCHMOCTH OT YacTOTbI pa3JIMYHbIX MOLE 
(be309eKTpHyecKHX d31, Tu3eKTpu4yeckux £43 HW yIpyrux 
Moyle THOKOCTH S14,512,5,3). PaccMoTpeHbI paqHalIbHEIe 
KojleOaHHaA OOpa3ziia B (opMe WHcKa Cc 9IeKTpOAaMu Ha TOp- 
wax. TommuHa o6pa3ia — | MM, qHuametp — 40 mM, quaria- 
30H KoneOaHHi — yo 700 KIu. Cuayana uccreqopamacb 
yacTOTHad 3aBHCHMOCTS JIA yIIpyrHx Moye KepaMHKH V3 
onpeyeueHua JeciTH pe30HaHCHbIX 4acTOT. 3aTeM M10 H3Me- 
PCHHBIM 3HayeHHAM TIPOBOAMMOcTH OOpa3ila Oba onpeye- 
jleHa 3aBHCHMOCTb OT YacTOTBI Mosyueli dz, u $3. C 9dTObI 
I[eJIbYO HCHOMb3OBaIOCh MOJyYeHHOe B padoTe BbIpaxKeHHe 
AIA SICKTpHYeCKON MpOBOAMMOCTH M3 pelieHHA paasIbHBIX 
Kojle0aHHli JMcKa C y4eTOM eFo TOJMIHMHEI. 
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Discussion and Conclusions. A technique is developed for 
determining the frequency dependence of LZT piezoelectric 
ceramic modules. The disc-shaped sample was studied in 15- 
650 KHz frequency range. It is shown that in the range up to 
650 KHz, $11,512,513 elastic modules with E superscript (it 
is omitted) or measured at dc field are practically independ- 
ent of frequency. In the specified range, ¢43, d31, kp con- 
stants have an insignificant frequency dependence for the 


considered radial oscillations. 


Keywords: piezoelectric ceramics, electroelastic modules, 
equivalent circuit, frequency dependence, disc, finite-element 
technique, perturbation method. 


For citation: V.V. Madorsky, et al. Study on frequency de- 
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O6cyacdenue u 3axmouenua. Pa3spadoTaHa MeTOaMKa TIA 
oMmpeeseHuA 3aBHCHMOCTH OT YaCTOTEI MO/Yleli Mbe303IeK- 
Tpuyeckoi KepamuKu I[TC. OOpa3er B popMe ANcKa ucce- 
WoBasica pH Wuama30ue yactot 15—650 KTu. Toxa3ano, uro 
B quana30He 20 650 KI ynpyrue Mody $44,542,513 C 
BepXHHM HHyeKcoM £ (oH onyjeH) WIM W3MepeHHble mpu 
TIOCTOAHHOM 9JIEKTPHYCCKOM Hosle WpakTHueckKH He 3aBHCAT 
oT YacToTH. B yka3aHHOM jMalla30He KOHCTAHTE! €43, d31, ky, 
II pacCMaTpHBaeMBIX paMasIbHBIX KoONeOaHHii HMeIOT He- 
3HAYHTeJIBHYIO YACTOTHYIO 3ABHCHMOCTB. 


KsroueBbie c10Ba: Mbe303IeKTpH4ecKkad KepaMHka, 9JICKT- 
poympyrve MO/lyIM, cxXeMa 3aMeII[CHHA, 3aBHCHMOCTb OT 
YacTOTBI, JHCK, M€TOJ, KOHCUHBIX 9IEMeHTOB, ANSYS, MetToy, 
BO3MYIIeHHH. 


O6pa3zeu dna yumupoeanua: Mayopcxuii, B. B. Uccnex 
OBaHHe 3aBHCHMOCTH OT 4aCTOTHI KOHCTAHT MOAPH30BaHHOH 
IIbe30KepaMHKH B CX€MAX 3aMEII[CHHA TPH caObIx WIeKTpH- 
yeckHx nomax (4acTp III)/ B.B. Manopcxui, VW. E. Poros, 
A. K. Kpyrsios // Bectuuk JJou. roc. TexH. yH-Ta. — 2019. — 


T. 19, Ne 2. — C.143-150. https://doi.org/10.23947/1992- 
5980-2019-19-2-143-150 


Introduction. A significant number of papers study and develop mathematical methods for solving problems 
of piezoelectric body oscillations. Being a component of the piezoelectric devices, piezoelectric elements serve to excite 
and register oscillations caused by external effects. The selection of piezoelectric material for measuring transducers 
and the analysis of their characteristics requires a large amount of information on the parameters of materials. Such in- 
formation includes: 

- a complete set of electroelastic modules [1], 

- losses, 

- analogous electrical circuits, or equivalent circuits of piezoelectric elements [2]. 

The reactive dynamic parameters (LZ, C) of the equivalent circuits are determined by elastic dielectric and pie- 
zoelectric modules, as well as by the density of piezoceramics [3]. It is considered that the electroelastic modules - the 
equivalent circuit parameters are constant; they do not depend on the amplitude (weak electric fields) and on the oscilla- 
tion frequency. All this is reflected in the Russian [4] and foreign [5, 6] standards for determining a complete set of 
electroelastic piezoelectric ceramics modules. For example, to determine d3; piezo-module of a disk-shaped sample, it is 
necessary to carry out measurements in three frequency domains: in the first and second resonances, in the antireso- 
nance region, and at frequencies much below | kHz. Accordingly, it is assumed that when determining d3,, the mod- 
ules of the ceramic under study in the frequency range from 1 KHz to the second resonance are independent of frequen- 
cy. 

The work objective is to study the frequency dependence of electro-elastic ceramic modules. In this case, a 
disc-shaped sample from LZT (lead zirconate titanate) is used, as precisely this type of ceramics is most well known. 
The ANSYS program [7] is used to verify the correctness of the presented methodology. 

Materials and Methods. Consider a piezoelectric disk with the thickness of 2/ and the radius a. We introduce 
a cylindrical coordinate system (7, ©, z) where z axis coincides with the direction of the polarization axis. The coordi- 
nate plane z = 0 coincides with the midplane of the disk. Based on the known linear piezoelectric ratios, the equations of 
the continuum dynamics [1] and Maxwell equations [8], the system of equations for the axisymmetric oscillations of the 
piezoelectric disk can be written as follows: 


O,T +037, + + pw2U = 0, 
1*rr 3 °rz r p 


O,T;z + 03 Tzz +7 + pw*W =0; 0D, +0,D, += = 0. (1) 


Henceforward, the following notations and definitions are introduced: U, W are mechanical displacements along r, z 
coordinate axes, respectively; @ is circular frequency; p is density; T,,,.,., are mechanical stresses; D,, D, are compo- 
nents of the electric induction vector; 0, and 03 are operators of derivatives with respect to r and z. 
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In case of axial polarization, the linear piezoelectric effect equations for weak electric fields in the cylindrical 
coordinates can be written as follows [1]: 


Try = C410,U + C12 < + C1303 W + €3103 &, 

Tog = C120,U + Cy - + €1303W + €3103 @, 

Tyz = C130,U + Cy3 - + C3303 W + €3303 @, (2) 
Dz = €31(0,U + =) + €3303 W — £3303 @; D, = €15(03U + 0,W) — £1101 . 


The following notation is introduced in the relations (2) and further: c,,, are elastic constants in the matrix des- 
ignation measured on samples with shorted electrodes or at the constant electric field E (E upper index of the ceramics 
constants is omitted here and further); e,,,, are piezoelectric constants; ¢,,, are clamped dielectric constants; © is electric 
potential where £ =—grad o [2, 8]. 

Assume that the boundary conditions are set at the electrode ends and on the lateral surfaces of the disk [1, 2]: 

atz=+h T,, =0; T,, = 0; p = 2V, (3) 
atr=a T,z, = 0; T,, = 0; D, = 0. (4) 
In (3), 2V value is electric potential difference applied to the ends [2, 8]. Let us introduce dimensionless coordinates and 

















quantities from the formulas: 
r z h p o. p cis 
=-;0=7;€=-; Q=oh [—; 0, =a fp; Cy, = C1. —-—. 
5 a’ 6 h’ a’ Cha ch? 11 Bea 4 


The solution to the boundary value problem (1-3) consists of the sums of two solutions: 








- homogeneous solution at zero boundary conditions at z = +h; 

- particular solution that satisfies only nonzero conditions at the ends (3). 

The constructed system of homogeneous solutions will enable to further satisfy the boundary conditions (4) on 
the lateral surface (as a rule, using variational methods). 

It is not difficult to construct the particular solution D2 = const and D? = 0, which satisfies automatically the 
third equation from the system (1) and the boundary conditions at the ends (3). According to the first two equations of 
the system (1), the mechanical and electrical components of the particular solution are equal: 








£33C 
To. = Tg9 = A(ea1 + = ~*)y sin(x@) + €3;K; U° =0; T2, = 0; 
33 
° = Kh¢ + Ahsin(xZ); D? = —€33K; TZ, = e33K + e33Ax cos(xG); ; 
a €33 K= V 1 (5) 
heb. (xcos(y) —k? sinQ))’ A 1 72 tan) 
é— a 
Xx 
The following notation is introduced in (5): 
X= 0. ae 64 = @33 + ae: k2 = ea 
33 33 D: 


If the vector of external forces and the electric potential are zero on the end surfaces, then the construction of 
homogeneous solutions is associated with the definition of the dispersion equation roots [9]. For symmetric oscillations, 
the dispersion equation has the following form: 

a,M, tan-1(B,) = 0, (n= 1, 2, 3). (6) 

The following notation is introduced in (6): 

An = O° Cy3kin + Cagk2n + C33k3ni Dn = KinBn + Kani 
M, = bgk33 — b3k32; Mz = b3k31 — byk33; M3 = byk32 — bk33. 
Here, Kj» are algebraic compliments to the elements of the third determinant row of the system for symmetric oscilla- 
tions (1); B,, is binary cubic root from [9], a is non-dimensional wave number of oscillations along r axis. 

To find the roots (a) at the given values of Q, it is necessary to solve the dispersion equation (6) in combina- 
tion with the binary cubic. A detailed analysis of the dispersion equation roots of symmetric lossless oscillations for the 
piezoelectric layer is given in [9], considering losses — in [10]. It is the sum of the particular and homogeneous solutions 
that will allow satisfying the boundary conditions both at the ends and on the lateral surfaces of the disk. 

With an arbitrary ratio of disk sizes, the inverse problem of its forced oscillations (1-4) is very difficult to ana- 
lyze and has a finite analytical solution only in some special cases (for example, oscillations of a thin disk or without 
considering its thickness when ¢ « 1). Therefore, when determining the modules of ceramics, it is more convenient to 
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solve the inverse boundary problem using approximate methods with allowance for thickness corrections for low- 
frequency radial oscillations of the disk. In this case, we can get an analytical solution in the form of finite formulas. In 
this paper, we seek the solution in the form of expansion with respect to ¢ small parameter: 

a? = £702 + yet 8 + ne... U = Uy + £702U2 + e*04U;, +... (7) 








Here, U,, is vector function with U(U, W, ~) components; n is the order of the constructed approximate theory of sym- 
metric disk oscillations considered in the paper; y, 7 are unknown constants which depend on the modules of piezoelec- 
tric ceramics, they are determined from (1, 2) and satisfy the homogeneous (zero) boundary conditions (3). 

Omit the cumbersome manipulations. We confine ourselves to the terms with e° to determine a” wave number, 
the terms with e* for U vector function in (7), and we give the final result of the considered boundary value problem 
with the boundary conditions (3) at the ends [11]: 








2 2 
— &13 _ 512, pn 2d34 ~, 13 S813, 
2? ~ P~ oT , <1 = , 
3¢33 Si’ €33(S11 + S12) C33, S41 + S42 
2 Me 
=yl=—t + piezo; 
eed er Sous prow 
t? d31 


piezo = ————- [ c7,(2td3, +d33)]?. 


a5ch,e3(1-K2) suits. 
Here, S,,, are flexibility modules with E constant; d3, is piezo-module; k, is planar coupling coefficient; v is planar 
Poisson ratio; £33 is free dielectric constant of the disk. 

Next, we introduce the following definitions and notation: a = e702 is approximate zero-order wave number; 
a3 = 702 + ye*04 is approximate wave quadratic number; a3 = €70,,7 + ye*0,,* + ne°Q,° is approximate quartic wave 
number; C is disk capacity. 

Table | presents the results of the exact solution of the wave number (a) of the PZT4 ceramics depending on the 
frequency from the dispersion equation (6) and an approximate calculation from (7) at ¢ = 0.033 for various approxima- 


tions of Qn. 
Table 1 


An example of calculating wave number (a) at various frequencies for the piezoceramic disk under study 




















f, kHz Q a. from (6) Oo On Q4 c piezo = 0 04 
50 0.085 0.045055 0.045049 0.045055 0.045055 0.045055 
250 0.425 0.2262 0.2252 0.22603 0.22605 0.2261 
500 0.85 0.4578 0.4505 0.4568 0.4574 0.4575 
700 1.1903 0.6522 0.6307 0.64777 0.6514 0.6518 


























The results given in Table | show that for the disk 1 mm thick and with frequencies up to 700 KHz, the cal- 
culation of the wave numbers from the dispersion equation (6) and the approximate calculation for a, almost coincide, 
and the piezo correction can be neglected. In this case, the decomposition (7) for a? depends only on the disk geometry, 
density, and moduli of flexibility — 514,542, 513, 

Omitting the relatively cumbersome manipulations, we give the expression of Y conductivity for the piezoelec- 
tric disk. It is obtained from the approximate solution (7) for various 7. 
Y000 is zero approximation of ¢ = 0, or the known equation of radial oscillations of a zero-thickness disk: 

Ji) 
XJo(x) — A - wa) 


0 


Qa 
x=; Yo00 = w6 (1- KB + KCL +) 
E 


Y040 is approximate quartic wave number: 


= 4, = _ pp2 2 IG 
x=“4; 040 =wC (1 kp + kB(1 + v) =): (8) 
Y042 is approximate quartic wave number, quadratic particular solution: 
07e3, Ji(x) 
+k5(1 +v) asc aite} 
3£33C3s Jo(x) — 1 — wa) 
Y242 is approximate quartic wave number, quadratic homogeneous and particular solution: 


1 033013 
b@ es se): 


€31 C33 


Y042 = wl (: — k2 + Cag ae 


= + k2(1 +v) 





07 e2 3 
Y242 = w6 (1-63 + cee . 
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zn = [xJo(x) — (1 —v)Jy(x)] [1 + e%at (= + =)] - e?aBy xJo(x); (9) 
E33C d 
Eo cS 
= i €33 S41 1 Sy2 


Table 2 shows the numerical analytic calculation results. The conductivities of the exact solution (5, 6) and 
approximate solutions (8, 9) with respect to e? were found using the MATLAB program [12]. For the disk made of 
PZT4 [13] with the thickness of 1 mm and the diameter of 30 mm, numerical calculations were carried out in the AN- 
SYS and ACELAN systems [14]. The latter software package was developed by the employees of the Southern Federal 
University (SFU) and focused on the calculations of piezoelectric devices. 














Table 2 
Conductance calculation for various frequencies 
Frequency, kHz | ANSYS,Cm/m | Solution (5,6),Cm/m | Y000,Cm/m | Y040,Cm/m | Y042,Cm/m | ¥242,Cm/m 
350 5.03 5.032 4.9688 5.014 5.0457 5.028 
554.4 8.74169 8.741029 8.478 8.5864 8.73237 8.735 
700 10.98 10.97 10.376 10.682 10.957 10.96 





























Table 2 shows that in the construction of approximate theories of the type (8, 9), it is required to use the de- 
composition (7) to calculate the wave number at least of the second order. 

Research Results 

1. Study of the dependence of the elastic ceramic moduli on frequency. For a sufficiently well-studied PZT 
piezoceramics, we will investigate the frequency dependence of the following modules: piezoelectric d3,, dielectric £4, 
and flexibility modules - s,,,5,2,5,3. Consider low-frequency radial oscillations of a disc-shaped sample with elec- 
trodes on the ends. The sample thickness is 1 mm, the diameter is 40 mm, and the oscillation range is up to 700 KHz. 

First, we study the frequency dependence for elastic modules. To do this, according to [15], we precheck the 
first three resonant low frequencies f,.(the major resonance and its two overtones). Elastic constants are determined 
from the solution of the frequency equation (9) of the radial oscillations of the disk considering ¢ relative thickness: 
three equations for three unknown variables. It is ¢ relative thickness that distinguishes the above frequency equation 
(9) for radial oscillations of the finite thickness disk from the known frequency equation for radial oscillations of the 
disk with “zero thickness” [15]: 

R JfR) = (1 — v) A(R). 

The introduction of ¢ thickness correction into the solution of the known equation of radial oscillations of a 
disc-shaped sample improves the accuracy and measurement informativeness of the elastic modules. 

The elastic compliances determined by [11, 15] through the technique of three resonances for the considered 
ceramics turned out to be equal: 

Sy, = 12.29e- 12, 54. =-4.05e — 12, 5,3 =—-5.28e — 12. 
Table 3 shows the first ten resonant frequencies for the PZT19 ceramics. In the “Experiment” stock, the fre- 





quencies measured on the WK 6500B conductivity measuring device at the Institute of High Technologies and Pie- 
zotechnics at the Southern Federal University are given. In the “Analysis” line there are 511,512,513, calculated by the 
formula (9) for the frequency-independent elastic modules defined by [11, 15]. Measurement errors did not exceed the 
values recommended by the standard [4]. 











Table 3 
First ten resonant frequencies for PZT19 ceramics 
Resonant frequencies f,, KHz 
Experiment 50.8 | 132.8 | 211 2884 | 364.6 | 439.4 | 513.2 | 585.4 | 654.4 | 722.8 
Analysis 50.8 | 132.8 | 210.95 | 288.1 | 364.8 | 440.7 | 516 | 590.4 | 663 735.4 









































Table 3 shows that in the frequency range up to 600 KHz, the difference between the calculated and experi- 
mental data with constant elastic moduli does not exceed 1%, therefore the following conclusion can be made: “For the 
ceramics considered, $14,512,513 elastic modules with F superscript (it is omitted) or measured at a constant electric 
field are practically independent of the frequency in the range up to 650 kHz.” 

2. Study of the dependence on the frequency of d3, and £33 modules. We use the expression (8) for Y040 
conductance from the solution of the radial disk oscillations considering its thickness, as well as the values of imaginary 
parts of the complex conductivity of the piezoceramics under study measured at room temperature. In this case, it is 
possible to investigate the dependence of d3; and €3 modules on the frequency in the range from 10 KHz to 600 KHz. 
To determine two unknown d3; and £43, measure Y at two frequencies — f\, fs. The difference between d3; and €43 is 
selected so that we can neglect the dependence of d3, and £33 modules on the frequency in f|— f; range. 
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In this paper, two conditions are introduced: 

1) f:—fi = 200 Hz; 

2) the conductivities are calculated for an ideal piezoelectric body or without considering losses, therefore the 
frequencies of resonances or in the neighborhood of resonances are not taken into account. 

As a result, we obtain an easily solvable system of two linear equations with respect to two unknown Heu3- 
BeCTHBIX ky and €33. 

Fig. 1 shows kj, d3, and €43/€, dependences on 10-650 kHz frequencies. 
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Fig. 1. Piezoelectric parameters — frequency dependences 


Here, € is the dielectric constant of the vacuum. In the range from 15 to 650 KHz, Ks coupling coefficient 
first increases with frequency growth from 0.28 to 0.34, and then decreases to 0.26. d3; piezo-module has a similar fre- 
quency dependence: at 15 KHz, from —1.50e '° to maximum of 1.661e —'°, and then decreases to 1.31e '°. The relative 
dielectric constant decreases monotonically with increasing frequency from 1766 to 1455. 

Fig. 2 shows frequency dependences in the range of 1-10 KHz measured and calculated from the formula (9) 
of the ceramics capacity with the constant modules defined according to OST [4, 15]. 
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Fig. 2. Piezoceramics capacity — frequency dependence 


There is sharp difference between two PZT ceramics capacities for low frequencies. This probably explains the 
dependence of the ceramics parameters on the frequency below 15 KHz (see Fig. 1). This implies a higher degree of 
dispersion — first of all, of d3, piezo-module d31. This issue will be considered in greater detail in the next paper 
(€33T, d31, ky constants are investigated in the low-frequency range). 

Discussion and Conclusions. A technique has been developed for determining frequency dependence of the 
PZT piezoelectric ceramic modules. The disk-shaped sample was studied within the frequency range of 15—700 KHz. It 
is shown that in the range up to 650 KHz, $41,542,513 elastic modules with E superscript (it is omitted) or measured at 
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a constant electric field are practically independent of frequency. In the indicated range, £33, d31, k, constants for the 
radial oscillations considered have an insignificant frequency dependence. 

Only the frequency dependence of real part of the ceramic modules was studied, the losses were not consid- 
ered. Therefore, the experimentally obtained frequency dependences of the imaginary part of the ceramics conductivi- 
ties were measured far from resonances, where the effect of losses was absolutely null. Even so, in the future, both the 
real and imaginary parts of the modules must be considered. This means that losses are included. This problem is sup- 
posed to be studied considering the frequency dependence of the complex modules of the ceramics under study. 

In the present paper, the radial oscillations of the sample are considered to obtain more complete information 
when measuring a larger set of constants (except for the piezoelectric and dielectric modules, s,1, 512, 513 elastic mod- 
ules are measured). This technique can be extended to other forms of samples (rods, plates, etc.), since in these one- 
dimensional problems, there is a simple replacement of Bessel functions by trigonometric functions. 
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Introduction. Polymeric materials (PM) are increasingly used 
in various industries and agriculture. Under the action of UV 
light, PM are destroyed. UV stabilizers are used to protect PM 
from photodegradation. Their action is based on the absorption 
of the photoactive sunlight component or on the deactivation 
of excited molecules that have already absorbed a light quan- 
tum, as well as on the inhibition of dark light-induced reac- 
tions. The work objective is to provide a comparative analysis 
of compounds of the sterically hindered phenols (SHP) series 
and nitroxide radicals (NR) as the PM photodegradation inhib- 
itors. 

Materials and Methods. Reagents of the “purum’” grade, 
LDPE premium grade 15803-020 polyethylene film were used 
for the investigations. The stabilizer was applied to the film by 
dipping. Carbonyl groups in polyethylene were determined by 


IR spectroscopy. IR spectra were recorded on Varian-640 


instrument. 
Research Results. 2,4,6-tri-tert-amylphenol (1), 2-methyl-4,6- 
di-tert-butylphenol (2), 4-acetylamino-2,2,6,6- 


3-carboxamido-2,2,5,5- 


tetramethylpyrrolin-1-oxyl (4) were tested as stabilizers for 


tetramethylpiperidin-l-oxyl (3), 


photo-oxidative degradation of polyethylene. It is known that 
the accumulation of carbonyl and hydroxyl groups is recorded 
in PM samples under irradiation in the process of photodegra- 
dation through the IR spectroscopy. The absorption band of 
the carbonyl group at 1720 cm | appears in the IR spectra of 
oxidized polyethylene. The IR spectra analysis shows that the 
content of carbonyl groups in the check samples is significant- 


ly higher than in the samples treated by stabilizer solutions. 





“ The research is done within the frame of the independent R&D. 


Beedenue. lommmepuste Matepuast (IIM) sce mupe ucrosb- 
B3YIOTCA B pa3JIMYHBIX OTPaCJAX MPOMBILVICHHOCTH H CeJIbCKO- 
ro xo3aiictBa. log nelicrsuem Y@-csera IIM pa3pymatrotca. 
7Ina 3anmutsr TIM or dotogectpykyuu upuMensotca YO- 
ctaOvmm3atoppl. Mx yelicraue ocHoBaHO Ha MOrionjeHHu 
(OTOXMMMYeCKH AKTHBHOM KOMIOHCHTEI COJIHeYHOTO cBeTa 
WIM Ha JeakTHBalHH BO30yxKeCHHBIX MOJICKyII, yxKe MOT0- 
TMBIUMX KBaHT CBeTa, a TakoKe Ha TOPMO%KeHHH TeCMHOBBIX 
peak, HHYUMpoBaHHEIx cBeToM. Llemb paboTsr — cpaB- 
HUTCIbHBIM aHasIM3 COeqMHeHHM u3 pAyja MpocTpaHcTBeHHO- 
3aTPyHeHHBIX (enous (113M) u HuTpOKCHIBHBIX payMKasIoB 
(HP) B kayectBe uHTMOuTOpoB doTogxectpyKyHH ITM. 
Mamepuaavi u memoovi. Jia uccneyqoBaHva VCTIOb30BaIIM 
peaKTHBbI KBaIM@UKayMH «4, TMOMMITHICHOBYyIO IIeHKy 
TISB], B/e 15803-020. CraOunm3atop HaHocHM Ha MIeHKH 
MeTOJOM OKyHaHHa. KapOonuibuble rpynibl B NosMITHIeHe 
onpeyzemamm Metoyom MK-cnexrpockonun. WK-cnexrpsi pe- 
THcTpupoBaM Ha mpHoOope Varian-640. 

Pesyismame ucciedoeanua. B KayectBe cTaOnmM3aTOpoB 
(OTOOKMCIUTeIbHOM eCTpyKUHH MOWMSTHIeHa MCTIbITAHBI 
2,4,6-Tpu-TpeT-amunpenona (1), 2-mMeTuil-4,6-M-Tpet- 
OyTumdeHona (2), 4-aljeTHaMuHO-2,2,6,6- 
TeTpaMeTwmunepHanH-l-oxcuna (3), 3-KapOoKcamuyo- 
2,2,5,5-TeTpaMeTusimuppommu- | -oKxcusia (4). 

VispecTHo, 4To npu oOsyyeHun B OOpa3yjax IIM B mpomecce 
dotoyectpyKkuHu Metoyom MK-cnextpockonuu dbukcupyetca 
HakoleHHe KapOOHWIBHBIX MW THAPOKCHIIBHBIX rpyum. B UK- 
CHeKTpaxX OKHCICHHOrO HOMMITHIeHa MOABIAeTCA Tomoca 
MOrIOWjeHuA KapOOHMIBHON rpynnol upu 1720 em |. Anamu3 
VK-cnekTpoB moka3bIBaeT, YTO cOfepxKaHve KapOOHWIbHBIX 
Tpyll B KOHTPOJIBHBIX OOpa3zllax 3HAaYHTeIbHO BbIWe, YeM B 


oOpa3rjax, oOOpaboTaHHBIX pacTBopaMH cTaOMIM3aTOPOB. 
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Discussion and Conclusions. The experiments show that ni- 
troxide radicals of 3-carboxamido-2,2,5,5-tetramethylpyrrolin- 
l-oxyl series and 4-acetylamino-2,2,6,6-tetramethylpiperidin- 
1-oxyl are the best photostabilizers of polyethylene. Moreover, 
there is no significant difference between the radicals of the 
2,2,6,6-tetramethylpiperidine and 2,2,5,5-tetramethylpyrroline 
series. Sterically hindered phenols, under photodegradation, 
have a far smaller stabilizing effect, falling short of nitroxide 


radicals. 


Keywords: polyethylene, photostabilizers, phenols, nitroxide 
radicals. 


For citation: 1.Yu. Zhukova, et al. Comparative characteristics 
of sterically hindered phenols and nitroxide radicals as stabi- 
lizers of polyethylene photodegradation. Vestnik of DSTU, 


O6cyorcdenue u 3aki04uenuA. OKCIepHMeHTEI TOKa3aIu, 4TO 
TYWIMMH otToctabumm3aTopaMuH MONMSTHIeHa sABIIAIOTCA 
HUTPOKCHJIBHBIC 


pagukasl 3-KapOoKcaMuyo-2,2,5,5- 


TeTpaMeTwymuppommH-l-oxcu =u  4-alleTHIaMHHO-2,2,6,6- 
TeTpaMeTwmunepHyHH-1-oxcns. Upwaem Het cymjecrBeHHoi 
MOK]Y 2,2,6,6- 


TeTpaMeTwinunepHywHa u  2,2,5,5-rerpameTHImupposmua. 


pa3HHIbl paqukasamu paqa 
IIpoctpaHcTBeHHO-3aTpyHeHHbIe (PeHOIILI B YCIOBUAX (POTO- 
TeCTPyKUMM OKa3bIBalOT 3HAYMTeIbHO MeHbilee cTaOwIH3H- 


pyromree TevicTBue, yceTyliad HATPOKCHJIBHBIM pajiMKasiaM. 


Ksrouesble c10Ba: TOuMaTMIeH, (oTOcTabuM3aTOpsl, PeHo- 
JIbI, HATPOKCHJIBHbIC payuKaliBl. 


O6opaszey OA YUumuposanua: CpaBHUTellbHad XapakKTepucTu- 
Ka IIpOCTpaHCTBeHHO-3aTpYAHCHHBIX (PCHOOB HM HATPOKCHIJIb- 
HBIX payjMKasIOB Kak CTaOMNIM3aTOpOB *oTOAecTpyKUHH NosM- 
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Introduction. Polymeric materials (PM) are increasingly used in various industries and agriculture. With a 
sufficiently high resistance to climatic factors, the PM performance decreases when affected by humidity and tempera- 
ture, as well as the ultraviolet (UV) component of the solar spectrum: in this case, radical particles are generated in the 
polymer, and transformations that cause material destruction occur [1, 2]. If in a PM there are products of its oxidation, 
for example, ketones, then they are photoinitiators of the decomposition process of the polymer. To protect the PM, it is 
required to ensure the interaction of radicals and highly efficient light stabilizers [3]. For this purpose, UV stabilizers 
are used which can act as follows: 

- they absorb photoactive components of the sunlight; 

- they deactivate excited molecules that have already absorbed a quantum of light; 

- they inhibit dark light-induced reactions; 

- they destroy or deactivate photoactive impurities and photoreaction products [4, 5]. 

Sterically hindered phenols (SHP), sterically hindered amines (SHA), and nitroxide radicals (NR) are widely 
used as UV stabilizers for polymers [6, 7]. The SHP act as chain terminators. Since the addition of mobile hydrogen to 
the primary radicals of the polymer photolysis is an oxidative reaction, the stabilizers themselves are easily oxidized [3]. 
Under photodegradation, the SHA form stable NR interacting with the alkyls of polymers which causes the photodegra- 
dation chain breaking [2]. 

The work objective is to provide a comparative analysis of compounds from the SHP and NR series as the PM 
photodegradation inhibitors. 

Materials and Methods. Reagents of the “purum” grade of the Aldrich company, and high-pressure polyeth- 
ylene film (LDPE 15803-020, GOST 16.337-77) 0.2 mm thick were used for the investigations. Photo-oxidative degra- 
dation stabilizers of the compound: 

- from the NR series — 3-carboxamido-2,2,5,5-tetramethylpyrrolin-l-oxyl (1), 4-acetylamino-2,2,6,6- 
tetramethylpiperidine-1-oxyl (2); 

- from the SHP series — 2-methyl-4,6-di-tert-butylphenol (3) and 2,4,6-tri-tert-amyl-phenol (4). 

The surface treatment of the film for applying stabilizer solutions was carried out as follows: samples of 5x5 cm 
were fixed in a cassette without damaging the working part of the sample, washed in warm water, blotted with cotton 
cloth, degreased in acetone and dried for three minutes in air at 20 —25°C. 

To prepare stabilizer solutions, 50 ml of toluene were poured in each of twelve 100 ml capacity measuring 
flasks, and the stabilizers (1-4) were added in a volume of 0.01; 0.005; 0.0025 mol. Then, the contents of the flasks 
were made up to the mark with a solvent. The resulting solutions were used to process film samples. 

The prepared film was dipped into a beaker with the solution of a UV stabilizer [8, 9] of the desired concentra- 
tion for 30-40 min. Then, the sample was removed from the solution and kept over the beaker for 3—5 minutes so that 
the excess solution drained from the treated surface. The sample was completely dried being suspended. 

For each stabilizer solution concentration, 5 film samples were tested. 120 samples were tested in total: 60 
pieces — under the natural conditions and as many — under the simulated conditions. 10 check samples (CS) were not 
treated with stabilizers. 
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After testing, the IR spectra of the check and experimental samples were taken. IR spectroscopy determined 

the content of carbonyl groups in the samples [10, 11]. IR spectra were recorded on Varian 640 instrument. 
Research Results. Field tests are the most reliable method for determining the PM resistance to UV exposure. 
To reduce the experiment time, quick tests in a lab environment are used. The samples were subjected to photodegrada- 
tion with a UV stabilizer ((1—4) compounds) and without it. In this case, the samples were affected by: 
- the sun, a natural source of UV rays — within 365 days (field tests according to GOST 9.708-83); 
- a xenon lamp, an artificial source of UV rays (flux density of UV radiation of 68 W/m” in the range of 
wavelengths of 300-400 nm, GOST 9.708-83) — for 504 hours. 

It is known that in the process of photodegradation through IR-spectroscopy, the carbonyl and hydroxyl groups 
are accumulated in the PM samples [1, 12]. But the yield of hydroxyl groups is a few dozen times less than carbonyl 
groups[12], therefore, the rate and degree of the polymer degradation is determined by the concentration of carbonyl 
groups according to the respective absorption bands in the IR spectrum [10]. 

The analytical IR-spectroscopy application to the determination of the film oxidation rate is based on the dif- 
ference in the spectra of the processed and not processed (check) samples. After UV irradiation, the absorption band of 
the carbonyl group at 1720 cm ' (stretch vibrations of the carbonyl group are C = O) appears in the IR spectra of poly- 
ethylene (PE). This band is analytical for the quantitative analysis. At the maximum of the absorption band at 1720 cm™ 
' the optical density value was determined by the baseline method. Calculations of the content of carbonyl groups in the 
samples were carried out by the well-known methods [10]. The research results are shown in Tables | and 2. 

Table 1 
Carbonyl groups content in polyethylene samples * treated with stabilizer solutions, under irradiation of xenon lamp, 
artificial ultraviolet source 





Stabilizer solution 
No concentration 
mol/l 0.1M 0.05 M 0.025 M 
Stabilizer 


name 





Content of C = O groups, % 





1 3-carboxamido-2,2,5,5-tetramethyl-pyrrolin-|-oxyl 




































































(1) 0.02 
2 4-acetylamino-2,2,6,6-tetramethyl-piperidine-1- = = 0.025 
oxyl (2) 
3 2-methyl-4,6-di-tert-butylphenol (3) 0.10 0.25 0.40 
4 2,4,6-tri-tert-amylphenol (4) 0.10 0.30 0.45 
5 Check sample” 2.50 2.50 2.50 
* No carbonyl groups were found in the sample before testing. 
** Not treated with stabilizer solution. 
Table 2 
Carbonyl groups content in polyethylene samples* treated with stabilizer solutions, under irradiation of the 
sun, natural ultraviolet source (exposure period was 365 days) 
Stabilizer solution 
No concentration 
mol/l 0.1M 0.05 M 0.025 M 
Stabilizer 
name Content of C = O groups, % 
1 3-carboxamido-2,2,5,5-tetramethyl-pyrrolin-1-oxyl = = 0.016 
(1) 
2 4-acetylamino-2,2,6,6-tetramethyl-piperidine-1- = _ 0.021 
oxyl (2) 
3 2-methyl-4,6-di-tert-butylphenol (3) 0.07 0.095 0.16 
4 2,4,6-tri-tert-amylphenol (4) 0.06 0.09 0.17 
5 Check sample” 0.35 0.35 0.35 





* No carbonyl groups were found in the sample before testing. 


** Not treated with stabilizer solution. 
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The IR spectra analysis showed that the content of carbonyl groups in the check samples was significantly 
higher than in the treated samples. Fig. 1 shows the IR spectra of CS and the samples treated with 0.025 M solutions of 
stabilizers (2) and (4) after irradiation with an artificial ultraviolet source. 
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c) 
Fig. 1. IR spectra of PE-film samples: treated with 0.025 M solution of 4-acetylamino-2,2,6,6-tetramethylpiperidine-1-oxyl (2) (a); 
treated with 0.025 M solution of 2,4,6-tri-tert-amylphenol (4) (4); check sample (c). Absorption level is shown on the right vertical 
scale. 
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The onset of the degradation process in the sample treated with SHP (4) is characterized by the appearance of a 
low-intensity absorption band at 1720 cm ' in the IR spectrum (Fig. 1, b). Such a band is practically absent in the sam- 
ple treated with the stabilizer (2) from the NR series (Fig. 1, a). Consequently, the contribution of NR to the decelera- 
tion of photodecomposition processes is much higher than of the SHP. 

UV radiation destroyed the CS. This is evidenced by the effect of accumulation of carbonyl groups in the sam- 
ple manifested in an increase in the intensity of the absorption band at 1720 cm ' (Fig. 1, c). 

Using IR spectroscopy, it was shown that UV-irradiation of PE films causes the accumulation of carbonyl 
groups associated with photodecomposition of PE carbon chains. Fig. 2 shows the kinetic curves of the accumulation of 
carbonyl groups during the simulated film photooxydation depending on the exposure time. 
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Fig. 2. Dependence of accumulation rate of carbonyl groups on time of simulated UV irradiation of samples treated with: 

1 — 0.025 M solution of 3-carboxamido-2,2,5,5-tetramethylpyrrolin-1-oxyl; 2 — 0.025 M solution of 4-acetylamino-2,2,6,6- 

tetramethylpiperidine-1-oxyl; 3 — 0.025 M solution of 2-methyl-4,6-di-tert-butylphenol; 4 — 0.025 M solution of 2,4,6-tri- 
tert-amylphenol; 5 — check sample was not processed 


The analysis of the data obtained shows that NR (1, 2) and SHP (3, 4) essentially slow down the photoinitiated 
oxidation destruction of PE and reduce the number of breaks of macromolecules. This is evidenced by the appearance 
of an induction period with an increase in the exposure time (see Fig. 2). 

At the initial stage of the experiment, changes in CS are observed. Gradually, the concentration of carbonyl 
groups in the film not treated with a stabilizer increases. This is associated with an increase in the rate of oxidation pro- 
cesses causing the structure destruction of the PE carbon chains and the formation of low molecular weight active com- 
ponents (radical particles). In the PE samples treated with stabilizers (1-4), no changes are observed. 

After about 150 hours, carbonyl groups are recorded in the samples treated with SHP (3, 4) solutions; after 220 
hours — in the samples treated with NR solutions (1, 2). The accumulation rate of carbonyl groups increases linearly 
within 150-200 hours. Further, the dependence is maintained during the whole process of photo-oxidation (up to 504 
hours). 

The induction period in the presence of NR and SHP can be explained by the ability of stabilizers (1-4) to cap- 
ture radicals quantitatively from the first moment; as a result, the process of destruction is inhibited at this stage. In this 
case, NR quickly interact with alkyl radicals, and, involving the polymer, they are converted into hydroxylamines capa- 
ble of restoring radical particles [6, 7, 12]. SHP can turn into quinones or mono- and diesters of hydroquinones [6], 
which are also traps for radical particles. 

The samples treated and not treated with stabilizers were also irradiated by a natural ultraviolet source accord- 
ing to GOST 9.708-83 for 365 days (see Table 2). The indicator of natural photo-oxidation of PE films was monitored 
after 1, 3, 6, 9, and 12 months (Fig. 3). 
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Fig. 3. Monitoring of indicator of concentration of carbonyl groups in PE film samples treated with 0.025 M stabilizer 
solution (1—4) and check sample (5) 


Comparison of samples treated with NR (1, 2) and SHP (3, 4) shows the following: in the first case, the con- 
centration of carbonyl groups increases insignificantly, in the second case, the content of carbonyl groups is sequence 
higher. In CS, a significant increase in the concentration of carbonyl groups was recorded after 1 month. 

It should be also noted that in the course of testing of CS, PE films became brittle, and at the end of the exper- 
iment, they broke down, in contrast to the samples treated with stabilizer solutions. 

Discussion and Conclusions. It is established that the best photo stabilizers of polyethylene are nitroxide radi- 
cals. Significant difference between the radicals of the 2,2,6,6-tetramethylpiperidine and 2,2,5,5-tetramethylpyrroline 
was not detected. Under the conditions of photodegradation, the SHP have a less stabilizing effect, yielding to nitroxide 
radicals. 

The inhibitory effect of NR and SHP is associated with their ability to interact with radical particles formed as 
a result of the polymer photodegradation. 

The tested compounds from NR class effectively slow down the process of photo-oxidation degradation and 
can be recommended for application as polymer stabilizers. 
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Study of the working body mechanism in forging-and-stamping equipment’ 


K. O. Kobzev! 


' Don State Technical University, Rostov-on-Don, Russian Federation 


VWsyuenne Mexanu3ma pagouero oprana B Ky3He4HO-IITaMMOBOuHOM OOopys0BanHun 


K. O. Koéses! 


' Jlouckoii rocynapcTBeHHEIi TexHuYecKHi yHuBepcute;r, r. Pocros-Ha-Jlouy, Poccuiickaa Denepauna 


Beedenue. CtaTba MocBaAlleHa UCCiIeqOBaHHIO CHHTe3Hpo- 
BaHHOM MpMHUMMMaIbHOW CXeMbI ()pHKUMOHHOTO KOHTaKTa 
TBEPAbIX Tel B KY3HCUHO-LITAMNOBOUHBIX MaluMHax. Ycta- 
HOBJICHa BO3MO2KHOCTb HOJIyYYeHHA MaKCHMyMa Harpy304HOl 
XapaKTepHCTHKH (PPHKUHOHHOTO KOHTaKkTa BHYTPH HHTepBa- 
Ila W3MeHeHHA Koopduunenta Tpenua. BelaBseHbI WBe cye- 
TYy!OWHe BO3MOXKHOCTH CHI TpeHHA (PpHKUMOHHOTO KOHTAaK- 
Ta: Ha rpaHwMwax yKa3aHHOro MHTepBasia OHM OyAYT paBHEI 
IIpH HasIM4H MaKCHMyMa paBeHCTBa; IPH JaHHbIX yCJIOBHAX 
OHH JOCTHraIoT HaHOobUlelt CTaOWIbHOCTH. 

Mamepuaavi u memoovi. Wipu v3MeHeHHu BesIM4YHHBI yrsia 
MeHACTCA MOOKCHHE TOUKH MaKCHMyMa. DITO NpHBOAMT K 
HapylIeHHIO paBeHCTBa CHJI TpeHHA Ha rpaHulax HHTepBasa 
W3MeHeHHA KooppuuNeHtTa TpeHus. B Takom cytyyae Ko93p- 
(PHUHeHT TOUHOCTH OMKeH OMpeseIATECA OTHOLICHHeEM 
MaKcHMyMa (YHKUMH K HaHMeHbUIeMy TpaHHyHomMy 3Haye- 
Huto. Ja sTOro HeOOXOAMMO ycTaHOBHTb TeHeHUHM H3Me- 
HeHHA TpaHWyHbIx 3HayeHH (yHKUMH, CBA3aHHbIe C BapbHi- 
poBaHHeM BeuIMUMHEI yruia. Jina pelienua sTO 3aqaqH HO- 
By!IO BevIM4HHy TaHTeHca yrla JaBNeHHA TpeycTaBHIH B 
Be MpousBeyqeHuA KOIpPpuyHeHTa BAapbupoBaHHsA Ha 0a30- 
Boe 3Ha4eHHe TaHreHca yrla. 

Pesyismamei ucciedoeanun. Homy4enubie pe3yibTaTbl 110- 
Ka3bIBaIOT BbICOKYIO CTaOHJIbHOCTb CHJIbI TpeHHA Ip Tmpo- 
CKaJIb3bIBaHHM Tel PPHKUMOHHOTO KoHTakta (OK). OqHako 
TIpH OoOMLWIMX BeIMYHHAX yryia WaBIIeHHA YyBCTBUTeNIbHBIX 
QICMCHTOB JaTavKa-lpeoOpa3z0BaTeaA MaKCHMalIbHad CHa 
TpeHHA KpaTKOBpeMeHHO MOxeT OBITh MponopuMoHuabya 
TeKyLIeMy 3HayeHHIO KOIPPUuLMeHTa TpeHHA. 

O6cyarcdenue u 3akniouenua. MoxepHu3upoBanHasd TIpHHyy- 
lHalibHad cxemMa MK no3BollaeT TeopeTH4eCcKH NOT 
OUCH BEICOKYIO CTAOHJIBHOCTh CHJIbI TpeHua. DK He JoIKeH 
oOpaljaTbca B HOJIb B MHTepBasle M3MeHeHHA KoopPuuneHTa 
TpeHHA BBIXONHOTO MapaMeTpa OCHOBHOM dpvKIMOHHOM 
rpyuner (OM!) u npu Hamm4uHu MakcuMyMa (pyHKIMM Harpy- 
304HOH criocobHoctu MK. HeoOxogzuMsIM ycioBuemM 3TOTO 
ABJLACTCA Mepeyaua YYBCTBUTCJIbHBIMH 3JIEMCHTaMH JOMOJI- 
HUTeIbHOH (pvkyHOHHOM rpynmer (JT) ee nonnor 
Harpy3KH. 


“ The research is done within the frame of the independent R&D. 
* B-mail: 5976765@mail.ru 
““ Padota BBITIOJHeHa B paMKax MHHWMaTuBHOH HUMP. 


Introduction. The synthesized basic diagram of the frictional 
contact of solids in forging-and-stamping machines is consid- 
ered. The possibility of obtaining the maximum load charac- 
teristics of the frictional contact within the variation interval of 
the friction factor is determined. The following two possibili- 
ties of frictional contact forces are indicated: they will be 
equal at the boundaries of the specified interval if there is 
maximum balance; they achieve the greatest stability under 
these conditions. 

Materials and Methods. When the angle value changes, the 
position of the maximum point changes. This causes violation 
of the friction forces balance at the boundaries of the variation 
interval of the friction factor. In this case, the accuracy coeffi- 
cient should be determined by the ratio of the maximum of 
function to the least boundary value. Doing this requires estab- 
lishing trends of changing the boundary function values asso- 
ciated with the angle variation. To solve this problem, a new 
value of the pressure angle tangent was presented as a product 
of the coefficient of variation by the base value of the angle 
tangent. 

Research Results. The results show high stability of the fric- 
tion force under slipping of the frictional contact (FC) bodies. 
However, at large values of the pressure angle of sensing ele- 
ments of the transducer, the maximum friction force can be for 
a short moment proportional to the current value of the friction 
factor. 

Discussion and Conclusions. The upgraded FC basic diagram 
enables to theoretically obtain very high stability of the fric- 
tion force. The FC should not vanish within the variation in- 
terval of the friction factor of the output parameter of the basic 
friction group (BFG) and at the maximum of function of the 
FC load capacity. A necessary condition for that is the transfer 
of full capacity of the additional friction group (AFG) by the 


sensing elements. 
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Introduction. The synthesized schematic diagram of frictional contact (FC) of solids in forging-and-stamping 
machines is considered. The analysis shows the possibility of obtaining maximum load FC characteristics within the 
variation interval of the friction factor. The following two possibilities of frictional contact forces are identified: they 
will be equal at the boundaries of the specified interval if there is maximum balance; they achieve the greatest stability 
under these conditions. The results obtained show high stability of the friction force during slipping of FC bodies. How- 
ever, at large values of the pressure angle of the sensitive elements of the transducer sensor, the maximum friction force 
can be briefly proportional to the current value of the friction factor. 

Materials and Methods. The FC scheme shown in Fig. | is free from this disadvantage. 


10 MF2 MFI 








Fig. 1. Basic diagram of frictional contact 


The master friction group (MF1) includes bodies 1, 2 and 4; the additional friction group (MF2) contains bodies 
3, 5 and 10. Between | and 2, 2 and 3 bodies, sensitive elements in the form of rolling bodies 7 and 8 are placed in the 
profiled sockets. 

The difference of the considered scheme is the division of body 2 into two parts, due to which the connection 
between bodies 2 and 5 is carried out using body 10 interacting with body 2 through the springs 11. This allows trans- 
ferring half of the total load of MF2 from body 2 directly to the rolling elements 8 making them a leading element in the 
MF2. The second half of the total load of MF2 is transmitted through the springs 11 from body 2 to body 10, and further 
through friction — to body 5. In addition, the place of F, moving force is transferred from body | to the rolling elements 


7. These differences allow excluding the impact of the pressure rolling 7 on the load distribution between the friction 
surfaces in MF1 and MF2 when the friction factor value changes [1]. 

To restrict to the linear nature of the friction between the bodies | and the supporting surface 9, as well as be- 
tween the bodies 2 and 10, the rollers 6 and 12 are installed. 

The formula for determining the value of the FC friction force depending on the disturbance input is as follows 
(at the same values of the feedback gain factor of MF1 and MF2]: 


4F, f 
F,=—,, 
uh (I+ figa:) 


where +, is total ultimate friction force between bodies 1, 2, 4 and 3, 5, 10; ; F,, is initial clamp force of friction 


(1) 


couples; f is current value of the friction factor between the mentioned bodies; a is pressure angle between rolling 


elements 7, 8 and the socket. 
It was also established that the function (1) within fiiin---fmax range of values does not have a maximum since 


the friction force of MF1 vanishes at the value of the friction factor of fg =1/tga. (Here, fin and finax are, respec- 


tively, the smallest and largest values that the friction factor can take in actual FC operating conditions for the adopted 
combination of materials of friction couples.) The same value corresponds to the imaginary maximum of the function 
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(1) since in f,..-fimax Interval, the FC load is transferred through MF2friction group, the friction force of which in- 
creases with f growth. This does not ensure the FC friction force stability despite the fact that it is somewhat higher 
than that of the first-generation FC (at values finin = 9-1, fimax = 9.8 and tga = 1.125, the accuracy coefficient is equal 
to K,=2.5 and K,= 3.67). 

The real maximum of the FC load characteristics within the variation interval of the disturbance input is possible 
through modifying MF2 according to the scheme (see Fig. 1). This modification is to reduce the number of friction cou- 
ples. To this end, it is necessary to exclude elements 10, 11 from the scheme and to rest body 5 directly upon body 2 


through the rollers 12. 
In accordance with this, we find: 


Fro = (Fo) 


where F7, is friction force of MF2; F,> is control action of MF2 (expansion force on the rolling elements 8). 


But 
Fo = Fy tga 7 
hence, 
Lad. 
Lg = . . (2) 
l+figa 


Friction force between MF1 pairs is equal to: 
F. 
Fy = Fr, [Ft Jtea| 1 , 

Considering the relation (2), we find 

=2F, ——_ (3) 
(I+ftga.) 

The relation (3) does not contain the difference in the numerator; therefore, this function does not vanish on 

any conditions. It only asymptotically approaches zero under the following conditions: 


- theoretically unlimited increase in friction factor; 
- a maximum at the point corresponding to f, =1/tga value [2-5]. 


Tl 


Summing up the friction forces in the relations (2) and (3), we obtain 
3+ figa 
DF =F, fee. (4) 
(I+ ftga) 
Differentiating the function (4) by / argument and equating the derivative to zero, we find the value of the 


friction factor corresponding to a maximum of the function: 


tee 3 


ie (5) 


on 


Substituting in sequence fin, /max Values in the expression (4), and equating the obtained relations to each 


— ati | © 


Here, m is relative width of the friction factor interval: m= finax / nin - 


other, we find: 





At the found tga value, the function (6), having a maximum at the point corresponding to the value (5), takes 


the same values at the boundaries of the variation interval of the friction factor. 
In this case, the coefficient of accuracy is determined by the ratio of a maximum of the function (4) to its any 


boundary value (for finin OT finax Values). Based on this, we get: 
9(+fmax 20)” ” 
8 fmax tga(3 a Fmaxtga) 


Set the value of parameter (tga ) at which K, value is minimal. When tga value changes, the position of the 





T 


maximum point ( /", ) changes as well. Only the relation (6) establishes the equality of the friction forces on the bound- 


aries of the variation interval, therefore, changing the position of the maximum point of the function (4) violates the 
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above equality. 

In this case, the accuracy coefficient should be determined by the ratio of a maximum of the function (4) to the 
smallest boundary value. To do this, it is required to establish trends in the boundary values of the function (4) associat- 
ed with tga value variation [6 —10]. 

To solve this problem, let us present a new value of the pressure tangent as a product of K variation coeffi- 
cient and the base value of the tangent angle in accordance with the expression (6), i.e.: 

tga; = Ktga. (8) 

Based on this, we can write considering (4): 

Gm + Kfmaxt8)fmax . B+K/maxt8)fmax 
(m+Kfrnax t@o) (1+ Kf max t@o.)” 
Here, the left side corresponds to the FC friction force at f = fiyi,, and the right side — to the friction force 





at f = fmax- The solution to the latter inequality with respect to the coefficient of variation is: 


2 
+1)- +1)" +12 
K 003 eee) is UC; +0). 
2f max tga. 





Since tga > 0, the fraction in parentheses of the left side of the solution is negative, therefore, for K > 1, the re- 
sulting inequality is undoubtedly fulfilled, i.e. when f, maximum point is shifted to smaller values, to calculate the 
accuracy coefficient, it is necessary to use the value of the FC friction force corresponding to fiy,, Value and converse- 
ly, when 1> K > 0, it is necessary to take the friction force value corresponding to /,,;, value [11-13]. 

Using this conclusion, we make the inequality of K,, > K., form where K,,, is the accuracy coefficient calculat- 
ed with consideration of the equality (8). We have: 

(14+ Kfimax 80)” (1+Fimax 1804)? 
KQB+Kfinax tga) 3+ finax tga 
The inequality is set up given that K > 1. 
Transformation of the resulting inequality to the form 


3+ (14K) finax tego. + Kf max tea > 0 
shows the validity of the assumption that K,, > K,. 





We investigate the ratio of the values of the coefficients of accuracy when the maximum point of the function (4) 
is shifted to the region of large values, that is, at K <1. Then, K,» > K, or 


2 2 
(m+ Kfnax tga) > (+finax t8) 
K(m+Kfrpaxt80.) 3+ frnaxt@O | 
Transformation of the resulting inequality to the form 


3m +(m+ K) frnax tga. — Kf max ta. > 0 





and the solution of the latter give 
m(3+finax t8t) 
(fmax tga. — Dfinax tga 
The fraction on the right side of the obtained solution is equal to one at the value 


(m+1)+4{(m+1)?+12m 
tga = 
2f max 
which fully corresponds to the solution (6). Therefore, the inequality K,. > K,is performed if K <1. 





K< 





’ 


Research Results. Let us take the final judgment on the stability of the FC friction force for various forms of 
load characteristics. For this purpose, we determine the FC accuracy coefficient when the function (4) has a maximum 
at f = fmax- In this case, the function monotonically increases in the range of the friction factor variation. At that, 


JSimax =3/tga equality is true. We will get: 


DF fin.) _ (m3)? 
Df.) 8(m + 1) 





Ky 
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Considering K,3 > K,,, we will find: 


3m(3 — mf max t Zo.) +m(3—m) f max tea. +91 = finax t@ox) <0. 
For all friction materials used as FC friction couples, m> 3. According to (6), finax tga >1, therefore, it is ob- 


vious that the differences in the brackets of the obtained inequality are negative and the above assumption is true. 

Thus, the analysis performed shows that the greatest stability of the FC output parameter is in the case when 
the function (4) has a maximum within the interval of the friction factor variation and takes the same values at its 
boundaries. 

When m= 8 and fray = 0.8, we get tga +14. Then, /’, 0.214. In this case, maximum of the function (3) will 


be at f,, 0.071 which is almost identical to the lower boundary of the interval f = fyi, = 0.1. The function (3) de- 
creases within the interval of the friction factor variation. With such initial parameters, K,~1.04 and K,3 ~1.68. 


Discussion and Conclusions. As is obvious, the modified FC block diagram enables theoretically to obtain 
very high stability of the friction force. However, due to the relatively large value of tga parameter, F,, force is used 


inefficiently. The FC should not vanish in the variation interval of the friction factor of the output parameter of the mas- 
ter friction group (MFG) and in the presence of maximum of the function of the FC load capacity. 

A necessary condition for this is the transfer of full load of the MFG by the sensitive elements of the additional 
friction group (FGD). An additional condition for the existence of a maximum can be formulated as follows: with an 
equal number of friction couples of both friction groups, the MFG sensitive elements transfer part of its full load; with a 
smaller number of friction couples than in the AFG, the sensitive elements transfer the full load of the MFG. 
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AnasM3 3PPeCKTHBHOCTH CHCTeMbI BHELIHerO OOOrpeBa MepHOH EMKOCTH HACOCHOM YCTaHOBKH OT BbIXJIONHOL 
CHCTeMBI ABHTaTesA BHYTpeHHero cropanusa 


C. O. Kupees', B. H. Crenanos’, M. B. Kopuaruna’, A. B. Epumos*™* 


1.2.3.4 Tlonckoli rocyapcTBeHHEIit TexHMyeckuii yHuBepcuter, Tr. Poctos-na-Jlony, Poccuiickas Dexepanus 


Introduction. The heating efficiency or regulation of the 
operating temperature of the fluid located in the gauge tank 
of a mobile pumping unit (PU) using the exhaust system of 
the internal combustion engine (ICE) is studied. The paper 
objective is to improve the design process and to reduce 
costs of field experiments. 

Materials and Methods. To solve this problem, numerical 
analysis methods were used when calculating the heating 
systems of the gauge tank with the external location of the 
pipeline from the ICE exhaust system by the example of the 
widely used installations (TsA-320, UNB, AChF, etc.). Sie- 
mens STAR-CCM+, a non-stationary non-linear solver of 
gas-dynamic processes, was used, which evaluates the cor- 
rectness 

of the problem statement and reduces significantly the costs 
of full-scale field tests. 

Research Results. The study was conducted for the operating 
conditions of the cementing unit on the Kamaz-43118 chas- 
sis with the SIN-32 pump and a drive from the power take- 
off attachment on the gearbox of the chassis engine. In the 
calculation model, a convection-type heat transfer was ap- 
plied between the body of the gauge tank and the surround- 
ing air; between the body of the gauge tank and the liquid; 
between the chassis exhaust system duct and ambient air; 
between the exhaust system duct and exhaust gases. The 
following study results were obtained: characteristic curves 
of the temperature variation of the gauge tank liquid at the 
watch points; the distribution of the temperature field of the 
liquid in the gauge tank; the distribution of the ambient air 
temperature field; lines of flows and ambient air speed field 
and of the ICE exhaust gases. 


Beedenue. CraTba MocBalleHa UccieqOBaHHIo 3:)peKTHBHO- 
cTH oOorpeBa HIM NoAepxKaHHA padouel TeMrepaTypbl 2%KUL- 
KOCTH, Haxoyjauelica B MepHOH eMKOCTH Tepe BwKHOM 
HacOCcHOM YCTaHOBKH, OT BbIXJIOMHOM CHCTeMbI J{BuraTesA 
BHyTpeHHero cropanna. Lem paborbl — coBepliieHCTBOBaHHe 
Tipowecca MpoeKTHpOBaHHA HW CHWKeHHe 3aTpaT Ha HaTypHBIe 
9KCIICpHMeCHTBI. 

Mamepuaaoi u memoooi. Jia petieHua TocTaBeHHoH 3aqauH 
ObLIM MCIOb30BaHbI MeTOJbI YACICHHOTO aHasIv3a TIpH pac- 
yeTe CHCTeM OOorpeBa MepHOl CMKOCTH C BHeINHMM paciio- 
JIOHKCHHEM TpyOoONpoBosa OT BBIXJIONHOM CHCTeMBI JIBHTaTesa 
BHyTpeHHero cropanua (JIBC) Ha puMepe mIMpoKo pacmpo- 
CTpaHeHHbIx ycraHoBoK I[A-320, YHb, AU® u gp. Wpumenen 
HeCTaliMOHapHbIi HeIMHeHHBIM pelllaresb ra30qHHaMHyecKHx 
TIporeccos (Siemens STAR-CCM+), no3BosA1ONMi OL[CHHTb 

IIpaBHJIBHOCTbh MOCTaHOBKH 3afauH MW 3HaYMTebHO CoKpalia- 
FOINM 3aTpaTbl Ha MOJHOMACcHITAaOHbIe HaTypHble MCIbITAHHA. 

Pe3yaomamoti uccredoeanua. UccneqoBpaHve NpoBoaHs0cb AIA 
ycCNOBHH paOoTbI IeMeHTHpOBOyHOrO arperaTa Ha IlaccH 
«Kama3-43118» c Hacocom «CMH-32» u npHBoyoM oT Kopo6b- 
KH oTOOpa MOIIHOCTH Ha KOpoOKe MepeKHOUeHHA Mepezad 
Tpuratesa Waccu. 

B pacueTHoi MojesIM IpHMeHeH KOHBeEKIMOHHBIM THI TerJIO- 
Tlepeqa4u Mex ]y KOpIlycOM MepHOH EMKOCTH HM OKpy2KaltouHM 
BO3YXOM; Me%K LY KOPIyCOM MepHOH CMKOCTH VM X2KHKOCTBIO; 
MOXY TPyOOMpOBOOM BBIXJIOMNHOM CHCTeMBI IaccH HU BO3- 
TyXOM OKpyKarolllel cpeybl; MexKZy TpyOonpoBoOM BBI- 
XJIOMHOM CHCTeMBI H BEIXJIOMHbIMH ra3aMu. 

B pe3ybTaTe HccieqOBaHHA TMOUy4eHbI rpapuyeckue 3aBH- 
CHMOCTH H3MeCHeCHHA TeMMepaTypbl %KHAKOCTH MepHOM eMKO- 
CTH B KOHTPOJIBHBIX TOUKaX; paciipeyemeHve MosA Temmepa- 
TYP 2KHKOCTH B MepHOM Oake; paciipezeseHve MosId Temre- 
paTyp OkpyKaloulero BO3yxXa; JIMHHM MOTOKOB H TOJIA CKOpO- 


cTe OKpykarollero BO3YyxXa HW BbIXJIOMHBIX ra30B JIBC. 
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Discussion and Conclusions. The study has revealed poor 
effectiveness of the considered structure. The data analysis 
allows us to offer an improved design of the gauge tank heat- 
ing. The results obtained in the work can be used in the cal- 
culations of such devices used in the hydraulic fracturing 


units and grout machines. 


Keywords: mobile pump unit, gauge tank, water feeding 
pump, high pressure pump, internal combustion engine 
(ICE). 


For citation: S.O. Kireev, et al. Effectiveness analysis of 
external heating system of gauge tank of pumping unit using 
exhaust piping of the IC engine. Vestnik of DSTU, 2019, 


O6cyacdenue u 3akuouenus. B xoye UccieqOBaHHA BbIACHH- 
Jlacb HefocTaTouHad 3(pPeKTHBHOCTh paccMaTpHBaeMOl KOH- 
cTpykuun. AHaiIM3 2aHHbIX MO3BOJIAeT MIpeIO%xKUTb ycoBep- 
II€HCTBOBaHHY!O KOHCTpyKuMio OGorpeBa MepHOii eMKOCTH. 
Pe3yibTaTbl, MOyYeHHbIe B padoTe, MOryT OBITb HCIIOIb30Ba- 
HbI TPH pacueTax MOMOOHBIX ycCTpOiicTB, MpHMeHAeMBIX B 
ycTaHOBKax Tupopa3pbiBa MlacTa MW IeMeHTHpOBaHHA CKBa- 
2KVH. 


KsoueBble c10Ba: TepeqByKHad HacocHasd yCcTaHOBKa, Mep- 
Had CMKOCTb, BOJOMOAaIOMM Hacoc, Hacoc BbICOKOrO jaBJIe- 
HHA, JBUraTeb BHYTpeHHero cropanua (JIBC). 


O6pazeu ona yYumupoeanua: Ananuz 39pdpeKTHBHOCTH cH- 
CTeMbI BHeLIHero OOorpeBa MepHOM eMKOCTH HacocHol ycTa- 
HOBKH OT BBIXJIONHOM CHCTeMbI JBHTaTe/IA BHYTPeHHero cro- 


vol. 19, no. 2, pp. 164-169. https://doi.org/10.23947/1992- panna / C. O. Kupees [nu gp.] // Becrumx JJon. roc. TexH. yH- 
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Introduction. The deposits of hard-to-recover hydrocarbons in the northern latitudes of Russia are confined to 
low-permeable, low-drained, heterogeneous and segregated reservoirs. Whereas it is necessary to develop such oil and 
gas areas, need arises for a guaranteed the above-zero temperature of the liquid that is in the gauge tank of the installa- 
tion. This problem remains relevant in the construction and workover of wells. 

The quality of hydraulic fracturing (HF) fluid, drilling fluid and cement slurry is affected by the following fac- 
tors: 

- types of rocks that interact with these systems; 

- well depth (pressure); 

- bottomhole temperature. 

We will accept two statements of the TNK-BP standard “Quality and Compliance under HF”: 

1) liquid chemicals specific for mixing at the well pad should be stored or transported above 15°C; 

2) storage, transportation, testing and injection of all liquid reagents should be carried out above 15°C, regardless 
of weather conditions [1]. 

The temperature of the system should not be lower than 8—10° C for all types of work on the preparation of the 
tempering fluid, pressurizing systems and systems of drilling and flushing solutions [2]. 

The 1979 oil and gas equipment literature mentions the use of a chassis or deck engine (TsA-320A cementing 
unit) for heating the exhaust gas system of the engine [3]. There are references to the vent dehumidifying system of the 
booster pulser hydraulic unit in two operating manuals: “Pumping units UNB-1000, UNB-800, UNB-630” and “Ce- 
menting hydraulic fracturing unit ACF-1050. UPETROM”. In the operating manuals of “Two-pump mobile unit UNP- 
320 x 40” and “Pumping unit UNB-125-50 SO”, the items on operating the drying and heating systems are highlighted. 

In the present paper, the effectiveness of external heating system of gauge tank of pumping unit using exhaust pip- 
ing of the IC engine is analyzed. 

Materials and Methods. Methods of numerical simulation were used. As examples, the widespread installations 
of the TsA-320, UNB, AChF, and others were considered. 

The application of numerical methods permits to assess the correctness of the problem statement (with an error 
of an idealized model selection), reduces significantly the cost of full-scale field tests or fully proves the design- 
engineering inefficiency of the problem. In the absence of statistical data on the problem under study, the results of nu- 
merical methods narrow the spectrum of the field problems to be solved. The non-stationary nonlinear solver of gas- 
dynamic processes (Siemens STAR-CCM +) [4] was applied as a numerical method tool. 

For the numerical study, the problem of heating a two- compartment gauge tank was selected. For this purpose, 
the cementing unit on the Kamaz-43118 chassis with the SIN-32 pump and a drive from the power take-off attachment 
on the gearbox of the chassis engine were used. In the rear zone of the chassis truck, capacity of 2 x 2 m? was installed. 
It is an all-welded sheet-metal construction which has two compartments with a sloping bottom and bottom valves. In 
the scheme proposed for calculation and analysis, the exhaust pipe (heating pipe) passes under the bottom of the gauge 
tank and rises upwards along the wall of the left tank (Fig. 1). 
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Watch points 1 Liquid Flow 
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P=1100 kg/m? T=5°C | v=lm/s T=5°C 
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bottom valve 







Exhaust gas flow 






Fig. 1. Idealized design model of cementing unit on “Kamaz-43118” chassis with “SIN-32” pump 


Research Results. The study was conducted for the following operating conditions of the unit at the field: 

- pump unit is fixed; 

- ambient temperature (air) is 5°C above zero; 

- wind speed is 1.0 m/s; 

- liquid density in both sections of the gauge tank is 1100 kg / m3; 

- liquid temperature in gauge tank is 5°C above zero. 

In the calculation model, a convection-type heat transfer was applied between the following media was used: 

- body of the gauge tank and the surrounding air; 

- body of the gauge tank and liquid; 

- chassis exhaust system duct and ambient air; 

- exhaust system duct and exhaust gases. 

To simplify the task, air was used as exhaust gases to eliminate the application of a multiphase solver. 

Air parameters were as follows: density was 1.184 kg m*; thermal conductivity was 0.026 W/mK; specific heat 
capacity was 1003.6 J/kg-K. 

The parameters of the liquid in the gauge tank were as follows: density was 1100 kg/m*; thermal conductivity 
was 0.569 W/ m:K; specific heat capacity was 4217 J/kg-K. 

The material parameters of the steel body of the gauge tank and the exhaust system duct were as follows: density 
was 7832 kg/m’; thermal conductivity was 63.9 W/m-K; specific heat capacity was 434.0 J/kg-K; thickness of the body 
sheet of the gauge tank was 4 mm; exhaust pipe thickness was 3 mm [5]. 

The exhaust gas outlet temperature of the ICE manifold was 450° C [5, 6]; the exhaust outlet flow of the engine 
manifold was 35 m/s [5]. 

The calculation (see Fig. 1) was made from the ICE working condition for 1800 seconds (30 minutes). The test 
points of the liquid temperature were obtained for the left and right compartments of the gauge tank at the bottom 
valves at a distance of 200 mm from the bottom, at a distance of 1000 mm from the bottom (T1), and at a distance of 
1000 mm from the bottom and 200 mm from the back wall (T2). The liquid temperature lines during the ICE operation 
for 30 minutes are shown in Fig. 2. The temperature fields and current lines are shown in Fig. 3-6. 
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Fig. 2. Liquid temperature lines of gauge tank at watch points: 1 is temperature at top point of left gauge-tank compartment; 
2 is temperature in left compartment (point 1); 3 is temperature in left compartment (point 2); 4 is temperature in left 
gauge-tank compartment at bottom valve; 5 is temperature in right compartment (point 1); 6 is temperature in right 
compartment (point 2); 7 is temperature in right gauge-tank compartment at bottom valve. 
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Fig. 3. Distribution of liquid temperature field in gauge tank (left view) 
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Fig. 4. Distribution of liquid temperature field in gauge tank (right view) 





Fig. 5. Distribution of ambient air and gauge-tank body temperature field 
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Fig. 6. Flow line and speed field of ambient air and ICE exhaust gases 


Discussion and Conclusions. The calculation results obviously demonstrate that for 30 minutes of the ICE op- 
eration and convection heat exchange between the hot exhaust pipe, air and liquid in the gauge tank at the test points at 
the bottom valves, the liquid is heated up to 0.1° C in the left compartment and to 0.1 ° C in right compartment. 

This brings us to the following conclusion: at the considered variant of the constructive solution, most of the 
thermal energy is directed to the operator’s cabin heating. This is a positive factor, but not a solution to the problem. 
The studies and graphical dependences show that it is required to place the exhaust tube directly in the gauge tank for a 
more efficient distribution of heat flows inside the gauge tank and in the environment. Given the design features, it is 
advisable to place the pipe inlet at the bottom of the right side of the tank and, passing through the left side, direct it 
along the tank up to the relief passage. 

In order to determine the possibility of heating the gauge tank, it is necessary to carry out calculations for a new 


schematic construction. 
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a ARK 
K Bompocy odecneueHua TepMeTH4HOCTH COCAHHCHHH Ha OCHOBE TexHOJIOIMUVeCKOrO HHAYUMpoOBaHHnn 


T. A. Homma’, I. P. Meipukos’, E. A. Wamduuos’, B. B. Kanycrun’ 


'23.4E DsSHCKMIi FOCy,apCTBEHHBIli TEXHH4eCKHii yHuBepcnter, Fr. BpaucK, Poccuiickas Dexepana 


Introduction. Some theoretical and engineering aspects of 
sealing joints through magnetostriction, as well as the polari- 
zation of the sealed medium under the external induction are 
considered. Control of surface roughness of joined parts to 
increase the joint density when induced by an external mag- 
netic field is studied. The creation of electromagnetic barriers 
for the migration of molecules of the sealed medium through 
a sealer is considered. The work objective is to validate the 
technological conditions for sealing movable joints in the 
cases described above. 

Materials and Methods. The conditions for ensuring the joint 
density are shown as a result of the contact problem solution 
and as a factor determined by the molecular-mechanical 
friction theory. Geometric, operational and tribological con- 
ditions of joint tightness are accepted. Damping properties of 
the fixed friction contact are determined by the molecular 
component. The theoretical and calculated analysis of the 
factors affecting the joint density is presented. Decrease in 
the smoothing depth, reduction of the ratio of transverse and 
longitudinal roughness steps, and increase in the contact area 
are indicated as the target results of the process preparation 
of the surfaces of the joint parts. Loss of tightness is defined 
as a specific transfer of molecules. They are transferred to 
the area of the joined surfaces or migrate freely through the 
sealer at the stages of sorption, diffusion and desorption. The 
predominance of any stage occurs when the entropy changes, 
and it is due to temperature and pressure. The schemes of 
sealing joints in the controlled magnetic field and of the de- 
pendence of magnetostriction and magnetostrictive stresses 
on the magnetic field strength are visualized. 

Research Results. The stability of sealers in highly volatile 
and gaseous media during their polarization and magnetiza- 
tion in an external field is experimentally investigated. In the 
former case, the magnetic induction vector was first oriented 
perpendicular to the longitudinal axis of the joint. A drop in 
the magnitude of the magnetic flux was observed when the 





* The research is done within the frame of Contract No. 9.10677.2018/11.12. 


Beedenue. CTaTba MOcBALleHa HCKOTOPbIM TeOpeTH4YeCCKHM HU 
TCXHOJIOFHYCCKHM aciiekTaM oOGecrieyeHHA TepMeTH4HOCTH 
cCoeqHHeEHHM Ha OCHOBe MarHHTOCTpHKIMEH, a TakKe NOuApH- 
3aljM TepMeTH3MpyeMOl CpeAbI B COCTOAHHM BHeIMHeTO HH- 
yyupopanua. ViccneyqopaHo yipaBlleHHe UlepOxOBaTOCTBIO 
TIOBEPXHOCTeH CTBIKYeMBIX eTaei TIA MOBbILMeHHA MIOTHO- 
CTH CTbIKa IIpH UX MHAYUMpOBaHHH BHeLIHHM MarHuTHbIM 
Trouiem. PaccMoTpeHo co3qaHve 3JIeKTPOMaTHUTHBIX OapbepoB 
WIA TepeMeljeHHA MosIeKyJI TepMeTH3HpyeMOl cpebI Yepe3 
repMetu3atop. Len padorrr — oOocHoBaHHe TexHouorH4ye- 
CKHX YyCNOBHH oOecleyeHHaA TepMeTHYHOCTH MOZBYKHBIX 
COeCAHHeHHH B ONMCaHHBIX BBIIMe CIydaix. 

Mamepuaavi u memoovi. Ycnosua oOecneyeHua TWIOTHOCTH 
CTbIKa TOKa3aHbl Kak pe3yJIbTaT PeWICHHA KOHTAKTHOM 3aqa4qH 
WM kak (akTop, OMpeeuAeMbIM NOMO%KCHHAMH MOJICKyIIAPHO- 
MexaHWyeckoli Teopuu TpeHua. II[punaTs! reomMeTpuueckue, 
OKCIIyaTalHOHHble HM TpHOoNOorM4ecKHe yCOBUA TepMeTH4- 
HOCTH coequHeHnH. Jlemndupyroue cBolcTBa Heo BMDKHO- 
TO (pHKUMOHHOTO KOHTaKTa OMpeyeyeHbI MOJIeKyIApHOH 
coctaBuaromei. IIpeycrapieHo TeoperuuecKkoe WH pacueTHOe 
oOocHOBaHHe (PakTOPOB, BUIMAIOWIMX Ha MIOTHOCTE CTEIKa. B 
Ka4eCTBe IeJIeBbIX Pe3YIbTATOB TEXHOMOTHYCCKON MOATOTOB- 
KH MOBepxHOcTeli JetTasel CoeqMHeHHA yka3aHbI yMeHBIIeHHe 
TIYOMHBI CriaxkKMBaHHA, YMCHbUIeCHHe OTHOMCHHA WlaroB T0- 
TepeyHoli H MpoOAOMbHOM WlepoxoBaTOcTH UM yBerMueHHe 
TIowjagqu KoHTakta. [lorepa1 repMeTHuHOCTH oMmpeyeeHa Kak 
chelyH@uueckHi TpaHcdep Moyekyn. Onn Tepenocatca B 
30HY CTBIKYeMBIX IIOBEPXHOCTel HI CBOOOTHO MepemMelyaroT- 
ca Yepe3 TepMeTH3aTOpP Ha CTaquAx COpOuMOHHOrO Norsowye- 
Hua, Dupdy3uu u yecopOunu. I[peo6naganue kakoli-1160 
cTaqMH IpOHCxOZHT Ip W3MeHeHHM IHTPONHU u OOycOBsIe- 
HO TeMMepaTypoi u aBieHuem. BusyamM3upoBaHbl CXeMbI 
TepMeTH3allHH COCAHHeEHHM B yIpaBJI4eMOM MarHHTHOM Mosle, 
3aBHCHMOCTH MarHHTOCTpHKUHH WM MarHHTOCTpHKUMOHHBIX 
Hallp@KeHUMM OT Hallps2»KeHHOCTH MarHHTHOPO MOA. 
Pesylemambel ucciedosanuA. IKCTIePHMCHTaJIbHO HCCIIeqOBa- 
Ha YCTOMYMBOCTL repMeTH3aTOPOB B JIerKOeTYYHX HW Ta30BbIx 
cpeflax Ip HX NOJApH3allHH MW HaMarHH4YMBaHH BO BHeLWHeM 
tome. B nepBomM cylyyae cHayasla BeEKTOp MarHHTHOM HMHyK- 
IMM OpHeHTHpoOBasIcA TepleHIHKyApHO TpoOONbHO ocn 
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compound was under the on-load operation for 268 hours. 
The total operating time of the joint was 1070 hours. If the 
magnetic induction vector was oriented longitudinally to the 
shaft axis, the operating time to the correction of the field 
strength was 87 hours. In the gas environment, the operating 
time of the connection to the adjustment of the tension was 
187 hours with a total operating time of 935 hours. 
Discussion and Conclusions. The penetrating ability of pres- 
surized media decreases in the “gas — vapor — liquid” series. 
It depends on the temperature at the joint contact. Depressur- 
ization can be traced through changes in the magnetic flux 
determined by the intrinsic magnetic permeability of the 
molecules of the sealed medium as they penetrate the inter- 
face surface. 

To increase tightness, it is required to suppress the activity of 
molecules. For this purpose, ionization and induction in the 
constant and alternating magnetic field with the intensity of 
<60 kA/m are used. 


Keywords: tightness, joint, machine parts, induction, densi- 
ty,magnetic field, crystallographic ordering, contact rigidity, 
anisotropic effects. 
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coequHenua. YcTaHoBsIeHO MaleHve BeJIMYMHbI MarHHTHOTO 
TOTOKa Mp padoTe coeyqHHeHHA MOA Harpy3KoH B TeyeHHe 
268 4acop. OOmad HapaboTKa coeqHHeHuaA cocTaBHa 
1070 uacos. Eciu xe BeKTOp MarHHTHOM HHYKIMM OpveHTH- 
POBaJICA TIPOAONbHO OCH Balla, HapaOOoTKa JO MOMeHTa Kop- 
PeKUHM BeIMYMHBI HallpsKeHHOCTH Moa cocTtaBusta 87 waco. 
B ra30Boli cpeyje NpOAOIMKUTeMBHOCTS paooTE! CoeqHHeHHA TO 
MOMeHTa KOppeKTHPOBKH HalipsxKeHHOCTH cocTaBua 187 
yacoB lip obujel HapaOoTKe 935 yacos. 

O6cyarcdenue u 3akmouenua. I[ponuKaroulai crocoOHOCTB 
TePMeTH3HPYeMBIX Cpe], YMeHbIMaeTcA B pay «ras — nap — 
2KUKOCTb»>. OHA 3aBHCHT OT TeMIepaTypbl Ha KOHTaKTe CO- 
equHenua. PasrepMeTH3al[HIO MO%KHO OTCJICXHTb 10 H3MeHe- 
HHAM MarHHTHOrO MOTOKa, olpeyessemoro coOcTBeHHO! 
MarHHTHOH MpOHMaeMOcTbIO MOJIeKy TepMeTH3Mpyemoli 
CpefbI IPH UX MpOHMKHOBeHHH Ha MOBepXHOCTb CTEIKa. UTo- 
Obl HOBBICHTh TePMETHYHOCTL, HEOOXOAMMO MOaBUTb aKTHB- 
HOCTb Mosekys. J[JIt STOTO IpHMeHAIOTCA HOHM3allHA W HHy- 
I{MpOBaHHe B HOCTOAHHOM HM TlepeMeHHOM MarHHTHOM Molle 
HallpsKeHHOCTBIO <60 KA/M. 


KsoueBbie c10Ba: TepMeTH4HOCTb, CoeqMHeHHe, WeTau 
MallHH, HHAyWMpoBahHe, MWIOTHOCTb, MarHHTHOe Tose, KpH- 
cTamorpapuyeckoe yiopsyOueHHe, KOHTaKTHad %KeCTKOCTb, 
aHH30TPOMHBIe 9(pPexTEI. 


O6pazeu 01a Yumupoeanua: K Bompocy obecreyeHHaA rep- 
MeCTHYHOCTH CoeqMHeHHM Ha OCHOBe TeXHOJIOrM4eCKOrO HH- 
ayyupopanua /T. A. umoumua [u gp.] // Becraux Jou. roc. 
TexH. yH-ta. — 2019. — T.19, Ne2. — C.170-178. 


https://doi.org/10.23947/1992-5980-2019-19-2-170-178 


Introduction. The theoretical basis of hermetology is presented in the works of G.A. Golubev, A.V. Chichi- 
nadze, V.P. Tikhomirov, L.A. Kondakova, G.V. Makarova, H.H. Wachter, E.Mayer, R.H.Waring, etc. It is known that 
tightness is mostly due to two factors: 


- stability of the processing density of abutting rough surfaces of the sealer parts; 


- stability of the sealer material properties and the sealed medium excluding the penetration and migration of 
molecules through its sealer beyond the interface. 


The conditions for ensuring the joint consistency are justified through solving the contact problem and are 
largely determined by the principles of the molecular-mechanical friction theory [1, 2, 3]. The mentioned approach re- 
mains one of the fundamental in hermiticity engineering. 

Penetration of sealed media (liquids, gases) through the sealer is described by the physical model of S. Slichter, 
I. Kozeni, T. Carman and D.K. Kollerov (for porous polymers). From its analysis, it follows that the volume flow of the 
medium through the seal is determined by the length of the free path of the molecules [4, 5, 6] and depends on the fol- 
lowing: 

-  microroughness height of solid surface - Rina, 
- joint density - P, 

- relative approximation - ¢, 

- geometrical dimensions of seal. 

Specifically, the fretting operation is initiated by boundary structures which are different polymer seals, liquids 
and films on the interfaces at the amplitudes of contact vibrations less than 0.05 mm. These media can act as catalysts; 
that causes a decrease in the technologically secured design value of the specific pressure on the seal and increases 
leakage. 

The amplitude and frequency of the relative displacements of the joint parts can be reduced through controlling 
the joint density, for example, through increasing the nominal pressure on the seal. The studies performed on polymers 
in the friction units [7] show that the joint density and contact rigidity depend largely on the compressive force of sur- 
faces and on their intermolecular activation: 
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Fp =aS+OP,, 


where F,,,, is friction force; a is average intensity of the molecular component of the friction force; Sis true contact area; 
b is coefficient of the mechanical component of the friction force; P,, is surface compression force. 

Consequently, the molecular component aS is proportional to: 

- the true contact area, 

- the intensity of mutual adhesion of contacted materials. 

Both of these indicators can be increased through plastic saturation along the interface plane, as well as during the 
alignment of the joint parts. Similar conditions can be created on the basis of magnetostrictive effects, which are pro- 
nounced to a greater degree for textured materials. Such conditions are applicable also in cases when reloading of joint 
parts is hindered or not provided for by the engineering performance standards. 

A decrease in the penetrating power of sealed media is predicted when a mobile spherulitic structure of the sealer is 
formed or modified by sealing with low-modulus materials [8]. However, the conditions that develop at drastically neg- 
ative temperatures and under cracking of the polymer sealer require special consideration. In these cases, differences in 
the temperature linear expansion coefficients with metals are observed; therefore, to control the permeability of the 
sealed medium molecules through creating electromagnetic barriers for them is of interest. 

Now then, this work objective is to study the technological conditions for ensuring the joint tightness based on 
the induction effects created in the materials of the joint parts and in a sealed environment. 

Materials and Methods. The conditions for joint tightness follow from the conclusions of the Lame problem, 
in which the boundary factors are: 

- geometric (for example, / - length of the contact of surfaces and d - diameter of the mating parts); 

- operational (R,.. - longitudinal axial force or M,,, - torque tending to move one part relative to another); 

- tribological (f- friction factor) determining the moment of pressing or turning. 

Friction locking is provided while minimizing the amplitude of the relative displacements of the surfaces, at 
contact stresses not exceeding the region of the damping capacity of materials [7]. 

Damping properties of the fixed friction contact are determined, in particular, by the molecular component. 
Moreover, the interface with the textured (anisotropic) structure in the direction of ensuring high damping capacity al- 
lows for a higher value of the critical vibration loads in a wide frequency range [8]. 

In addition, when an external induction changes R, roughness parameter due to the magnetostrictive effect, 
they can affect the joint density. This is explained by the comparability of magnetostriction values and the sizes of 
blocks of structural components (about 10 *-10~ m). Moreover, the technological texturing in crystallographic direc- 
tions with pronounced maxima of the elastic or plastic properties of materials causes a maximum magnetostriction ef- 
fect at the unidirectional position of the magnetic induction vector. This causes an increase of the nominal pressure val- 
ue in the preloaded demountable joints since part of the stresses spending on the parts' compliance is compensated by 
magnetostrictive stresses within the elastic properties. In permanent joints, on the contrary, part of the stresses can be 
spent on the plastic saturation of the contact. For the same reason, the actual preload is greater than the calculated one. 

P is joint density coefficient with magnetostrictive effect: 


P=(C marty +Vomcnedy Kandy ! Vobuy + (1) 


Here, Vi.nxn 18 the material volume of the rough layer with A&/ texture; V2.4 1s the material volume of an in- 
compressible seal extruded into the gaps; & is the magnetostriction coefficient; Vg, is the total volume of the rough lay- 
er. 

It is known that the presence of waves causes a decrease in the bearing contact area by a factor of 5—10 [9]. 
Therefore, the solution to the contact problem in the aspect of ensuring the joint tightness requires consideration of the 
waviness of the surfaces. 

The height and step parameters of roughness are interrelated; therefore the technological preparation of the sur- 
faces of the joint parts should be aimed at the following results: 

- reducing the depth of smoothing (R,); 

- reducing the ratio of steps of transverse and longitudinal roughness; 

- increasing the contact area and R,,/R, ratio. 

As atule, there is no direct relationship between the tolerance value and the height of irregularities; however, it 
should be noted that the depressions on the surface profilogram turn out to be reservoirs in which by-products accumu- 
late, which initiate surface weakening. 

While ensuring the joint density due to the preload, Sj compliance of the compound is represented as: 


Sy =x(1) + (Eqnry(Bk )/ N’ (2) 
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where x (f) is normal contact displacement in the function of time; &(,,, (B) is the magnetostriction tensor at a regulated 
crystallographic ordering considering plastic deformations in the coating and changing with magnetic induction; k* is 
coefficient of crystallographic ordering; N*is dynamic load. 

When the compounds operate under the hydrostatic pressure of the sealed medium (liquid or gas), loss of tight- 
ness is expressed in a certain form of transfer of molecules. They are transferred to the area of the joined surfaces or 
migrate freely through the sealer at the stages of sorption, diffusion and desorption. The predominance of any stage oc- 
curs when the entropy changes, and it is due to temperature and pressure. 

The simplest case of phase transfer in seals is the flow of a viscous incompressible fluid in a porous medium. The 
determination of the leakage value follows from the Darcy law. However, in this case, it is not considered that the seal- 
ing joint, in contrast to the porous solid body, is formed as a result of the contact of two surfaces. Therefore, under the 
load increase and the alignment of the surfaces, individual contact spots can form isolated volumes (effective channels). 
Their number decreases with increasing load until the formation of a hydrodynamic film along the seal perimeter [10- 
11]. Thus, considering the conductive capacity of most sealed media, it is advisable to investigate tightness from the 
standpoint of polarization and external induction. 

According to [12], when a medium passes through an electric current source or ionizer, molecules obtain q 
charge. If a charged particle with an initial velocity v moves in a uniform magnetic field along the magnetic induction 
lines B of the external field, then the angle a between the vectors of v and B is equal to zero. The magnetic field does 
not affect the particle (the Lorentz force is also zero), and it moves uniformly and rectilinearly. 

If a charged particle moves with velocity v in a magnetic field perpendicular to vector B, then the Lorentz 
force F’ = q[vB] is constant in magnitude and perpendicular to the particle trajectory, which prevents its displacement 


(Fig. 1, a). 
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Fig. 1. Schemes for sealing joints in controlled magnetic field: in uniform field of toroidal inductor (a); in inhomogeneous 


field solenoid (4); 1 is shaft; 2 is sleeve; 3 is magnet; 4 is charged medium particle; 5 is source of electric current (or ionizer); 6 is 
automatic solenoid control unit with voltage generator); v is orbit of positively charged medium particle 
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Assuming the direction of medium motion, knowing its electric potential, viscosity, density, penetrating ability, 
and changing the field strength and the position of the magnetic induction vector, it is possible to provide the Lorentz 
force in quantities sufficient to completely suppress the motion. Thus, the effect reactive to the hydrostatic pressure is 
realized. 

Similarly, in an inhomogeneous magnetic field, the component of magnetic induction B’ creates force F ', 
which pushes the particle into the region of a weak field (see Fig. 1, 5). Initially, the particle moves along the radius of 
constant magnitude with velocity v. However, its own magnetic field counteracts the externally induced one, and this 
explains the extreme instability of its state. As a result, the particle is pushed into the field region with a reduced volt- 
age. 

The sealing effect is based on creating — by the magnetic field — the reactive repulsion of pre-charged medium 
molecules penetrating the interface. Whereby, the degree of sealing depends on the magnetic inductance, the value of 
which is set through the strength of the current considering the relative magnetic permeability and form factor of the 
inductor. 

It can be assumed that the conditions of joint vacuumization are created when the surfaces of the mating parts 
come closer to the appearance of molecular (adhesive) interaction. In this case, the magnetic permeability of the materi- 
als of the joint components is represented by the ratio of the magnetic fluxes in the material ® and in vacuum 0 (p= 
@/@®,). Under the assumption of ® = Dg, the magnetic permeability of vacuum can be taken as a unit (u = 1). 

In the medium with constant magnetic permeability, the magnetic field induction is proportional to its density. 
Consequently, the magnetic flux from the external field, initially defined for vacuum, depends on the magnetic permea- 
bility of the medium. Its presence at the interface, due to its own magnetic permeability, either increases the magnetic 
flux (in paramagnetic media), or decreases it (in diamagnetic media). 

Thus, it is possible to evaluate and correct the performance of the sealer in the induction state through changing 
the magnetic flux and, accordingly, magnetic induction. 

Research Results. The crystallographic areas described below can be taken as regulated directions under the 
process part texturing. 

1. For detachable joints without a polymer sealer: [110] (for structures with body-centered cubic (bcc) lattice), 
[111] (with face-centered cubic (fcc) lattice), [1000] (with hexagonal close-packed (hcp) lattice). 

2. For permanent joints and connections with a sealer: [100] - [111] (with bcc lattice), [100] - [110] (with fcc 
lattice), [0001] - [2110] (with hcp lattice). 

When the magnetic induction vector is established in the indicated directions, the magnetostriction A in the 
crystals turns out to be field oriented. This is because its magnitude is non-linear to the field strength and is determined 
mainly by the position of the easiest magnetic axis [13, 14]. 

In the absence of preliminary texturing, a mean change in roughness for structures with cubic lattices is repre- 


sented as: 
= 9 F 
05 AR <( Zit 2h] . (3) 
5 5 
where Ajo, 4111 are magnetostriction constants along the [100] and [111] crystallographic directions; d* is an average 
size of the dispersed structural components. 


In the case of texturing in the directions [100] or [111], when the field is oriented along the normal to the surface 


(Fig. 2), the changes in roughness (R,) will be: 
OS AR. <Awd/k, (4) 


05 AR <2hnd/k, (5) 


where k is demagnetization loss factor. 
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Fig. 2. Magnetostriction A — magnetic field intensity H dependence 


Valves with a seat made of X12M steel with a coating in the direction [111] of CoFe,O, cobalt ferrite coating 
with thickness of 20x10 m (Rinax 2-4 wm, 7) = 52x10° and 12 = 50x 10° m) and a brass valve by the amount of leak- 
age flow of Q, uPa * m°/s of helium (7 = 77 K) are tested. It is established that, during the compound magnetization in 
the field of 60 kA/m, this value is 7830-10240. The number of shutter cycles is 200-5000, respectively. Under similar 
test conditions for a saddle-valve pair in the absence of a magnetic field and prior texturing of saddles, the O value was 
960-15340 p1Pa « m’/s, respectively. 

With an increase in the magnetic field strength, no major changes in the connection health are established. Ap- 
parently, this is due to the fact that effective magnetostriction stresses are formed in fields <40 kA/m (Fig. 3). This con- 


clusion is consistent with the results of the magnetostriction study. 


Steel P6MS 
Steel U1OA 
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0 5 10 15 20 25 30 35 H, kA/m 
Fig. 3. Magnitude of magnetostriction stresses o — magnetic field intensity H dependence 
Changes in the properties of the sealed medium, the materials of the sealer and the temperature affect the pene- 
trating ability and tightness in general. The permeability of media decreases in the “gas - vapor - liquid” series. With 
that in mind, the performance of sealers in volatile and gaseous media during their polarization and magnetization in an 


external field was investigated. 


Machine building and machine science 


175 


http://vestnik.donstu.ru 


176 


Vestnik of Don State Technical University. 2019. Vol. 19, no. 2, pp. 170-178. ISSN 1992-5980 eISSN 1992-6006 





At the end of a rigid detachable shaft — copper hub connection, a source of electrical current was attached, one 
of the poles of which was fixed to the ends of the hub. At that, the second pole had a minimum clearance with the shaft 
generatrix. 

The joint was placed between the poles of a toroidal electromagnet connected to the automatic control unit 
with a voltage generator, a voltmeter and a webermeter. The magnetic induction vector was oriented perpendicular to 
the longitudinal axis of the joint. The joint was installed in a chamber filled with diamagnetic benzene. Pressure of 30- 
35 kPa was developed in it. The sleeve was rigidly fixed, after which the source of electric current and the inductor 
were turned on, and the magnetic field strength was set at 60 kA/m. The shaft was mechanically loaded in the radial and 
axial direction relative to the sleeve in a symmetrical alternating cycle with the load of 1.5 N with the frequency of 50 
Hz. The tests were performed until the reduction of the magnetic flux recorded upon the webermeter indication. This 
moment characterized the violation of tightness expressed in the penetration of benzene to the joint surface. The induc- 
tor automatic control unit corrected the voltage in the toroid, due to what the magnetic field strength was set equal to the 
initially predetermined one. 

A drop in the magnitude of the magnetic flow was observed when the joint was under the on-load operation for 
268 hours. The total operating time of the joint was 1070 hours. The criterion for the completion of the tests was the 
reduction of the interval of operating time until the moment of correction of the magnitude of the magnetic force. It 
should exceed 60 kA/m. This is due to the need to increase the operating time to the value corresponding to the previous 
cycles. 

With the orientation of the magnetic induction vector longwise the axis of the shaft under similar test condi- 
tions, the operating time until the correction of the field strength was 87 hours. 

The simulation of the operating conditions of the joints in the gas environment was carried out in a similar 
way. An ionizer was mounted on the shaft which was attached to the ends of the sleeve. The joint was placed in a sole- 
noidal inductor, the density of the turns of which increased from the end of the sleeve. As a result, a non-uniform mag- 
netic field was generated, the intensity of which increased from the end of the solenoid. 

The joint was placed in a chamber under the pressure of 30-35 kPa with paramagnetic nitrogen, after which the 
ionizer and inductor were switched on. The shaft was loaded in the modes described above with periodic recording of 
the magnetic flow. The moment of changing its initial value characterized the violation of tightness and the penetration 
of nitrogen to the interface. In this case, the inductor automatic control unit provided a reduction in voltage in the sole- 
noid winding to the value corresponding to the magnetic field strength of 60 kA/m. 

The operating period of the joint until the adjustment of H value of was 187 hours with a total operating time 
of 935 hours. The long duration of the tests caused the correction of the magnetic field strength with the time interval of 
less than 187 hours. 

Discussion and Conclusions. Reduction of leakage through the seals is provided by the joint density (specific 
pressure p,) management while reducing the amplitude and frequency of the relative displacements of the joint parts. 
The latter is obtained through increasing the true contact area and the intensity of mutual adhesion of the contacted ma- 
terials, in particular due to magnetostriction effects. In this case, the level of magnetostriction stresses should not cause 
plastic deformations of the sealer (ie, p, ~ 1.1 6,,) and can be provided in the fields up to 40—50 kA/m. As the complex 
roughness parameter A changes due to magnetostriction, an increase in the contact bearing area is observed, mostly due 
to the material whose elastic modulus is lower and the magnetostriction is higher. 

The impact of the bearing contact surface area, which increases with magnetostriction, grows with decreasing 
the smoothing depth R, and the ratio of of transverse and longitudinal roughness steps. As a result, the nominal pressure 
in the joints grows due to the compensatory effect from the side of magnetostrictive voltages. Magnetostriction is stimu- 
lated by crystallographic magnetic anisotropy, in particular under technological texturing. In the “gas - vapor - liquid” 
series, the penetrating power of sealed media decreases. It is dependent on the temperature at the joint contact. Depres- 
surization is detected through changes in the magnetic flux under the influence of the intrinsic magnetic permeability of 
the molecules of the sealed medium as they penetrate the interface surface. You can increase the tightness if you sup- 
press the activity of molecules. For this purpose, ionization and induction in the constant and alternating magnetic field 
with the intensity of <60 kA/m are used. 
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Introduction. The problem of adequate modeling of physical 
processes and objects is an urgent task of modern science. 
Automation of work in this area increases the accuracy of 
modeling and saves money on full-scale modeling of objects 
under study. The research objective was to create a computer 
program that would automatically form dimensionless simi- 
larity criteria based on the initial set of parameters. 

Materials and Methods. The calculation algorithm is based 
on m-theorem of similarity theory and the basic provisions of 
the dimensional theory. 

Research Results. The computer program is developed. It 
enables to get all the theoretically possible similarity criteria 
for the simulation of the physical process or object of inter- 
est. 


Discussion and Conclusions. The results obtained can be 
used in the industrial and scientific modeling of physical 
objects of research, calculation of new similarity criteria, 
solving problems of describing complex processes, etc. 


Keywords: similarity criterion, theory of similarity, dimen- 
sional theory, theorem, dimensionless complex, Kirpichev- 
Gukhman theorem. 


For citation: 1.M. Tsibirova. Formation of similarity criteria 
for physical objects and processes based on NonDimCrit- 
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Beedenue. CoppeMeHHad HayKa IIpv3HaeT akTyaJIbHOM mpo- 
Onemy ajleKBaTHOrO MOjerIMpoBaHHA (pu3H4ecKHX MpoleccoB 
wu OObekToB. ABTomMaTH3alMa paOor B 9TOH ciepe MOBbIMIAaeT 
TOUHOCTb MOJeMpOBaHHA H IKOHOMHT cpecTBa Ha HaTypHOM 
MoyeuMpoBaHHH OObeKTOB H3yyeHHaA. Llerbro WaHHOrO Uccie- 
OBaHHA ObIIO CO3LaHHe KOMIbIOTepHO! MporpamMMBbl, KOTO- 
pasa 10 HcxoqHOMy HaOopy HapaMeTpos aBToMaTuueckH (pop- 
MupyeT Oe3pa3MepHble KpHtTepuu no706u4. 

Mamepuaanei u memoovi. B ocHoBe pacueTHoro ajropuTma 
JI@KUT T-TeopeMa TeopHv MosOONA MH OCHOBHEIe MOOKeHHA 
TeOpHU pasMepHocTeii. 

Pe3ytbmamet uccredoganua. Co3faHa KOMIbIOTepHad TIpo- 
rpaMMa, KOTOpad HO3BOJIAeT MOJIYYHTh BCe TeopeTH4ecKH 
BO3MOXKHBIe KPHTepHu MoOA00uA AIA MOJeMpoBaHHA w3y4ae- 
Moro (bu3HyecKOro Mporecca HIM OObeKTA. 

O6cyoarcdenue u 3axmouenue. HomyaeHHble pe3yIbTaTbI MOryT 
ObITh HMCIOb30BaHbI B MPOMBILMJICHHOM HM Hay4HOM MOJfeJIH- 
poBaHHH (Pu3H4eCKHX OObEKTOB HCCeAOBaHHA, pacueTe HO- 
BbIX KPHTepveB MO WOOuA, pelieHHH 3a,ad ONMCAaHHA CIOKHBIX 
HIpoueccos UT. J. 


Ksrouesbie copa: KpuTepHii mozobus, Teopua moyzobus, 
TeopHA pa3sMepHoctTeli, 1-Teopema, Oe3pa3MepHbIi KOMIUIeKc, 
Teopema Kupmmuepa — Tyxmana. 


O6pazeu Ona yYUumuposeanua: LuOuposa, U. M. Dopmuposa- 
Hue KpHTepHes NosoOusa cbu3H4ecKHX OOLEKTOB MH MpOLeccosB 
Ha OcHOBe KOMIIbIOTepHOH TIporpaMMBI 
NonDimCritFormer 1.0/ W.M. Un6upopa // Bectauk Jon. 
roc. TexH. yH-Ta. — 2019. — T.19, Ne2. — C. 179-184. 
https://doi.org/10.23947/1992-5980-2019-19-2-179-184 


Introduction. Similarity criteria are dependencies between dimensionless combinations of parameters describ- 
ing an object or physical process. On the basis of similarity criteria, it is possible to create mathematical models of the 


objects of interest [1-4]. 


Similarity criteria enable to establish the correspondence of the model and the object under study. They be- 


come irreplaceable when the mathematical description of processes according to experiments or observations is not 
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formed. In fact, similarity criteria are dimensionless power complexes of sets of physical parameters describing objects 
and processes. They are identical in form and numerically equal for similar objects or processes. 

The algorithm is based on the second theorem of similarity theory, the so-called 2-theorem. In the source [5], it 
is formulated as follows: any complete equation of the physical process, written in a certain system of units, can be rep- 
resented by a functional dependence between the similarity criteria obtained from the parameters involved in the pro- 
cess. From this, it follows that analyzing and combining with each other the physical parameters describing the object 
of study, it is possible to make the similarity criteria of the object or physical process [5, 6]. 

When modeling objects of research on the basis of similarity theory, the analysis of the dimensions of the 
physical parameters describing the object at its primary points [6, 7] is most frequently used. 

Researchers select the determining physical parameters that describe the object of interest, in particular, pres- 
sure, viscosity, density, speed, temperature, resistance, current, and, working with their dimensions, obtain a theoretical- 
ly possible number of similarity criteria [8]. These criteria can then be used as the basis for mathematical or full-scale 
modeling. 

In the case when the initial object under study and the model have similarity criteria that are numerically equal 
in the critically important zones, then, according to the third similarity theorem (Kirpichev -Gukhman theorem) [9], it 
implies their physical similarity [6]. Therefore, to create, for example, a model of flow along the contour of pipes, it is 
necessary to find the appearance and numerical value of the similarity criteria describing the real object, and then to 
build a model or laboratory installation so that the values of the similarity criteria in it are numerically equal to the val- 
ues of the same criteria on the real object. 

The NonDimCritFormer 1.0 program presented in this paper enables a person who studies a physical object or 
process to select significant physical parameters (see Fig. 1), which determine its basic parameters and calculate auto- 
matically all theoretically possible similarity criteria and dimensionless complexes (see Fig. 2) formed on their basis. 
The basis of this algorithm is the methods of analysis of dimensions of the selected physical parameters. 

2. Description of the program algorithm 

The basis of the algorithm implemented in the program of formation of similarity criteria based on the dimen- 
sions of the selected physical parameters is a technique that includes five steps. 

Step 1. Determination of physical parameters that affect significantly the physical properties of the object un- 
der study. The accuracy of the results depends on the correct choice of parameters, neglecting nonessential ones. A 
large number of selected parameters causes an increase in the number of obtained dimensionless complexes, but does 
not increase the accuracy of the model formed on their basis. 

Step 2. Compiling a matrix of dimensions for the selected physical parameters. In it, the dimensions of the se- 
lected parameters in the corresponding degrees are arranged in rows. The dimension of the matrix corresponds to the 
number of parameters and their units of measurement, that is, in the matrix there are as many rows as the selected pa- 
rameters, and there are as many columns as the units of measurement. For example, a string corresponding to the physi- 

2 


kg 


cal parameter “electrical resistance, R”, which has a dimension of four units =o 
A 





, B MaTpulle pa3sMepHocteit OyzeT 
s 


uMeTb Buy (2 | -3 -2). 

Step 3. Determination of the number of independent parameters. The number of independent parameters k cor- 
responds to the rank of the dimension matrix, i.e. the largest order of the non-zero determinant which is based on the 
rows of the matrix under study. 

Step 4. Formation of groups of independent parameters. At this step, combinations of k pieces of m parameters 
that make up a dimensionless complex are formed. That is, the units that are included in the dimension of the selected 
parameters are mutually reduced, which leads to the fact that this combination of the selected physical parameters has 
no dimension. Generally speaking, the complex consists of such parameters, the partial matrix of dimensions of which 
has at least one determinant of k-th order unequal to zero. 

According to the combinatorial placement formula, the theoretically possible number of such complexes can be [1] 

ck _ m! 
kt (m—k)! 
where m is the total number of parameters, k is the rank of the complete matrix of dimensions. 
Not all obtained complexes of the selected parameters make up a complex with a matrix of k-th order, there- 


fore, very often the number of groups after the test is less than Ci ; 


Step 5. Determination of expressions for similarity criteria based on calculated groups of independent parame- 
ters. Expressions are formed by the following algorithm. 
1. Take another group of parameters. For example, set #1- P,...P,. 


Tsibirova I. M. Formation of similarity criteria for physical objects and processes based on NonDimCritFormer 1.0 computer program 





2. We begin the calculation of similarity criteria with dividing the parameters of the next group by the parame- 
ters of the current group. In particular, the following criteria can be drawn up for the set #1 


Fa _ En 
Ty = is ve) Wink : 
Ma, Pv. Perk a Perks 
PMP Fatt. PI 





Ordinary division of a parameter by a set does not produce a dimensionless expression. To do this, you need to 
calculate the values of powers x; }...x;;, for each of the criterion expression. 


3. Considering that the dimensions of the parameters in the numerators and denominators of the criterion com- 
plexes should be equal, we form and solve a system of linear equations. Obtaining such a system is due to the equaliza- 
tion of degrees at the same units of measurement in the numerators and denominators of the criteria, which are the 
products of combinations of dimensions selected at the beginning of the calculation of physical parameters. As a result, 
we obtain a system of g equations with k = q unknowns. 

We obtain sets of m-k criteria for each group of parameters. The maximum possible number of dimensionless 
complexes, which can be calculated on the basis of m parameters, is 


n, =C* .(m—k). 
This value can be quite large. For example, for m = 7 and k = 3, the total number is 1, =140. But, due to the 


fact that the list of calculated criteria has a lot of the similar and reciprocal ones, the total number of selected criteria is 
less than 2, . The proportion of criteria which failed verification depends on the initial set of physical parameters select- 


ed for study. 

3. Program interface and functionality 

The program NonDimCritFormer 1.0 is written in the language VisualBasic.Net. It has an intuitive Windows 
interface and is easy to operate. 

The user chooses physical parameters, which are essential, in his opinion, for the simulated object or process 
(see Fig. 1), then, by clicking the mouse on the button “Generate all criteria” and “Remove repetitions of criteria’, all 
similarity criteria (dimensionless complexes), which can be made on the basis of the selected physical parameters (see 
Fig. 2) are displayed. 


ws. CritMaker 1.0 -|5) x} 
atin Mpaska Tectb! CTaHaapTHble KpHTepHn 
BeiGepure wcxogHble napamerpbi 
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Hanpsikenne, U Di ckopocte punbTpauun, wFilt 
HaTA@DXKeHHE NOBEpKHOCT., Fpoy [M ckopocts, w 

ofmem, ¥ OC conpotwen. snextpuy., R 
MAOTHOCTe TOKA, rol D conpotnenenne ya., RUd 
MAOTHOCTE, 10 Co temneparypa, T 

nnowanp, S Ci tennoemn, maccosaa ya., OCmUd 
NoTeHuWan TepmaguHam., potTD C1 tennoemn. monapHaa ya., QCMolyrUd 
noTeHuHan xuMHYeCKH, potHim Otennoem«. offoemnaa ya., QCVUd 
notok TennoBoi, potocO C1 tennoemkocts cucremer, QC 
NpoBogumocth ya., GUd OC tonwuna, th 

nposoaumocts 3neKrpyu., Gel Clyrnosaa cxopoct, ugw 
NpoBoguMoct, G C1 yrnopoe yckopenue, uga 
npowssogurenbHocts ya. PrevUd =] yron naockuii, ugPl 

patora,A Clyron tenechi, ug 

pacxoa Maccosbili, rashm [4 yckopeHve cBo6.nageHua, g 
pacxog ofpemubiii, rash DO yckopenue, a 

cuna caeta, Fsv QO yacrora, ch 

cua ToKa, | O sksueanenr 3 nextpoxum., 
cuna, F UO snextponposoae. ya., GUd2 
CKOPOCT SBYKa 6 cpene, Ww2v Osneprua, E 

cKOpocTb Maccosaa, wm QO serponua cueremei, Entrp 
CkKOpocTb Macconepea., wmPer 





Jasnenue, p 
Javamerp. d 


emkocTb 3nekTpHy., C 

SaPAN SNEKTPHYECKH, ge 
WHQYKTHEHOCTE, L 
KOnMYeCTBO BewectBa, kolv 
Konw4ecteo Tennoter, kol@ 
KOHWeHTpayHA Maccosaa, Cm, 
KOHWUeHTPaUHA Ofpemnaa, CV 
Koa @.Temneparyponpos, kT 
Koshpuu. auedysuu, kD 
Kos@onu. maccooraayn, kmot 
Kospuu. Macconepegayn, kmper 
Kos epuu. TennonposaaH., kk. 
Kos epuu.tennoorgaun, kh 
Macca monapxaa, mMolyr 
macea, m 





OOOOOOOODOOOOOOO00090OU000 


Cpopmuposarb Bce 








U6pare noeropey | Hairy e | Dobaeurs Hoses Upanur Copopmuposare Fre 
KpHTEpHH Kpt 5 napamerp: BbIENeHHbIe sem ucKogHyio 6asy 
Fig. 1. Selection of initial physical parameters 
Sn athe 
@aiin Mpsexa Tecte: CraxaaprHeie prrepin 
BeiGepure HexoaHele napamerpes 


Kpurepuu 6e3 nopTopexnk 





Fig. 2. Calculated similarity criteria 
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Using the “Find standard criteria” button, you can detect the existing (named) similarity criteria in the list. 
NonDimCritFormer 1.0 also enables to print the entire list of criteria or part of it, to save the list as a separate file, to 
edit the base of standard criteria. 

4. Examples of using the program 

Calculation example 1. Under constant voltage conditions, we find similarity criteria describing the current 
flow in a resistor, choke, and capacitor connected in series [8]. 

Step 1. Select the physical parameters essential for the simulated object. They include: current, I, A ; induct- 

2 4 42 2 2 
ance, L, ee capacitance, C, oe resistance, R,; =; voltage U, a2 
s°-A kg-m s-A s--A 





and time ¢, s. Total parameters 


are m = 6 and base unit are g = 4. 
Step 2. The matrix of dimensions of the specified parameters has the form: 


kg s 
t}O O 1 0 
Cil-2 -1 4 2 
L|2 1 2 -2 
U/}2 1 -3 -l 
Z}0 0 O 1 
R|2 1 -3 -2 








Step 3. The rank of this matrix is k = 3. It follows that the number of independent parameters in the group will 
be three, and the number of criteria for each independent group will be m—k=3. 

Step 4. The total number of groups of independent parameters can be: a = res = 20. Combinations of pa- 
rameters are made automatically with the help of the described program. The algorithm checks the rank of the partial 


matrix for each composed combination. 
For example, for an ILU combination, a partial matrix would look like this: 


0 0 0 1 
IU=\;2 1 -2 -2 
2 1-3 -l 


The rank of this matrix is 3, therefore, this combination forms a group of independent parameters. Not all of 
the possible 20 combinations are dimensionless similarity criteria. For example, CL combination is not a criterion be- 
cause the rank of the matrix is not equal to x. In the context of the initial set of physical parameters under consideration, 
only 15 can be considered as similarity criteria. 

Step 5. Drawing up similarity criteria passes through the solution of systems of linear equations, which are 
formed on the basis of criterion complexes. For example, for a combination of /LU: 

fru. C™ R* 


ON ooh Fi peak an 8 Gee ye He? 
[2.73.4 [022 . [23 .U 4 [32 - £33 -U34 





after solving the system, which will result in finding the values of powers of dimensions x, the following crite- 
ria will be obtained: 
- t a Cc 7 R 
poe 9° Pa 








Ty 7 
| mee se 


For each of 15 possible combinations of initial physical parameters, 3 criteria will be obtained. Among these 
45 criteria, there are often repetitive and reciprocal ones. For this example, NonDimCritFormer 1.0 gives 11 unique 


criteria: 
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Calculation example 2. Consider the process of fluid motion in a horizontal pipe. To do this, we form similari- 
ty criteria that describe the pressure drop in this process [10]. The most important parameters that describe the pressure 
drop in the pipe during the passage of the liquid medium, the following values can be called: 


Ap = f (,p, HL), 


where Ap is the pressure, 





k, : ., ms, : ki P 
S, ; @ is the flow rate of the fluid, —; p is fluid density, = ; pis the coeffi- 
m:s S m 
. er . kg 
cient of dynamic viscosity, ———. 
m:s 
Total parameters are m = 5, and the base units are g = 3. The complete matrix of dimensions of the selected pa- 
rameters has the form: 


m kgs 
uwf-l 1. -l 
Ap j-l1 1 -2 
Lil 0 0 
p}3 1 
o/l O -l 








The rank of the reduced matrix is k = 3. The number of independent parameters in group is three and the num- 
ber of criteria for each independent group is m — k = 2. Total groups can be ce = C? =10. 


According to the calculations of the program, in this example, out of ten possible combinations, only nine are 
independent groups. It follows that a total of 18 dimensionless complexes can be obtained. After removing the repetitive 
and reciprocal ones, the program gives 6 unique combinations. 








0,5 
L- Ap Ap Laka’? <p" 
Ty = > IQg= >? 13 = 3 
H:® p-@ H 
po? L:p-@, LT? -Ap-p 
MH et oS 8 Se 
Ap”? HW u 


It is important to note that the complexes 2 and 75 are expressions of the Euler and Reynolds criteria, respec- 
tively [3]. 

4. Conclusion 

The computer program NonDimCritFormer 1.0 is proposed. With its help, on the basis of selected physical pa- 
rameters that determine the characteristics of the object under study or the physical process in its characteristic spatial 
and temporal points, it is possible to form dimensionless complexes — similarity criteria describing the object of study. 
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CpaBauTesbHbiii anasu3 MOJMMUUMPOBAaAHHON NOCTKBAaHTOBOH KpHnTorpapuyeckol cucrembr1 NTRUENcrypt 


set 


MW oOmjenpHAATON Kpunrocucrembl RSA 


II. B. Pasymos’, MH. A. Cmupnos’, M1. A. Inaunenko’, A. B. Cenépa’, JI. B. Uepkecona” 


1,2,3,4,5 


Introduction. The NTRUEncrypt cryptographic system, the 
calculation of the algorithmic complexity of the development 
of the NTRUEncrypt cryptosystem and its modifications are 
considered. The study objectives are to develop NTRUEn- 
crypt, an efficient post-quantum cryptographic algorithm, 
which has high cryptographic resistance to quantum comput- 
er attacks, to work out a modification of the proposed algo- 
rithm, to analyze and experimentally validate its advantages. 
Materials and Methods. A description of the NTRUEncrypt 
encryption system is proposed. The modification of the con- 
sidered algorithm is studied; the block diagram of the im- 
plementation of the software based on it is presented. An 
example of the software operation and its characteristic is 
given. The reliability of the results is proved using the 
Mann-Whitney U test. During the experiment, the third-party 
software implementation of the RSA cryptosystem was used. 
A Stopwatch class element was introduced in the source code 
of all three programs of NTRUEncrypt, RSA, and NTRUEn- 
crypt modifications. This class provides a set of methods and 
properties that can be used for the precise measurement of 
the execution time. Thus, it became possible to record the 
results of the time spent on all three basic stages: key crea- 
tion, encryption and decryption of the message. 

Research Results. The advantages of the developed cryp- 
tosystems in terms of the performance characteristics are 
proved. An experimental comparison of the implemented 
NTRUEncrypt algorithm and its modification is performed. 
All advantages of the latter are indicated. 

Discussion and Conclusions. The advantage of using the 
NTRUEncrypt algorithm modification is experimentally 
validated. The new application is 25% faster to perform gen- 


eral work on key generation, encryption and decryption. In 


“The research is done within the frame of the independent R&D. 


Jlouckoii rocyyqapcTBeHHbI TeXHuYeCKHH yHuBepcutert, r. Poctos-Ha-Jlony, Poccniickaa Degepauna 


Beedenue. CtaTba NocBaljeHa UCCHeqOBaHHIO KpunTorpadu- 
yeckolH cucrembt NTRUEncrypt, pacuery auroputmuyeckoit 
CiO%KHOCTH pa3paboTKH KpuntTocuctempr NTRUEncrypt u ee 
MoguduKkayun. Lenmm uccneqoBpanua: paspadotKa 3pdekTuB- 
HOrO MOCTKBaHTOBOrO KpHMTorpauyeckoro ayIropHTMa 
NTRUEncrypt, o6mayqatoujero BbICOKOM KPHITOCTOHMKOCTBIO K 
aTaKaM C KBaHTOBOTO KOMMBbioTepa, a TakxKe pa3paOoTKa MO- 
WupukalHu WpeoxeHHOrO asIrOpHTMa, aHasiM3 MU IKCHepH- 
MCHTAJIBHOC JOKA3ATCJIBCTBO CFO IIPCHMyLIeCcTB. 

Mamepuaavi u memoovi. Upennoxeno onucanHe CHCTeMBI 
wumm@posanua NTRUEncrypt. Wsyyena moguduKanna pac- 
CMaTpHBaeMOTO asIrOpHTMa, peycTaBleHa OOK-cxema pea- 
JIM3al1MM OCHOBaHHOYLO Ha HEM TiporpaMMuHoro cpeyctsa. IIpu- 
BeleH TIpHMep paOoTs! NporpaMMHoro cpeyCTBa HW JaHa ero 
xapaktepucTuka. JIOcTOBepHOCTh pe3yIbTaTOB OOOCHOBaHAa C 
moMomb10 U-Kputepua Manna — Yuruu. IIpu npopeyzenuu 
9KCHepHMeHTa HCHOJIb30BaHa CTOPOHHAA NpOrpaMMHas peasin- 
3allda KpuntTorpapuyeckort cuctempl RSA. B ucxoqHerit Koy 
Bcex Tpex mporpamMmM NTRUEncrypt, RSA, mogzudbuKkanuu 
NTRUEncrypt Obi BHeypeH 9IeMeHT KIacca Stopwatch. J[an- 
HbIM Kacc peqoOcTaBaeT HAOOp MeTOOB H CBOMCTB, KOTOPBIe 
MO2KHO HCHOJIb30BaTb Jd TOUHOTO H3MepeHHaA BpeMeHH, 3aTpa- 
YeHHOTO Ha BbIMOHeHHe. TakuM 00pa30M, HOABHJIACb BO3MO%K- 
HOCTb (PUKCHPOBATh Pe3yIIbTATHI 3ATPaeHHOTO BPeMeHH Ha BCeX 
TpeX OCHOBHBIX 9Tallax: co3qaHve KOveH, WMppoBaHve u 
pacimM@popanne coobmjeHns. 

Pesyismamei uccaedoeanua. JioKa3aHbl TpeHMyllecTBa pa3- 
paOoTaHHBIxX KpHITOCHCTeM Ilo xapakKTepHCTHKaM IIpovu3BO]H- 
TCJIbHOCTH. BmOIHeHO 9KCIICPHMCHTaIbHOe CpaBHeHHe pea- 
JIM30BaHHOTO alroputTMa NTRUEncrypt u ero MoquduKayun. 
IIpu sToM o6o03Ha4eHBI BCe IpeHMylecTBa MocueqHeH. 
O6cyorcdenue u 3akiloveHUuA. DKCTICPHMCHTAIbHO OKa3aHoO 
TIPCHMyLIecTBO HCHOJb30BAaHHA MOMPUKAaYMN asropuTma 
NTRUEncrypt. Hopoe npunoxenue Ha 25 % ObicTpee BBITIOI- 
HAeT OOLIyIo paOoTy NO reHepalluu KrIO4eH, WIMdppoBaHHtIo u 
pacmimposannio. Ilomumo 3Toro ONTHMH3HpyeTcaA MCIONb- 
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addition, the internal memory usage is optimized through 
reducing the weight of the source program file and the size of 
the secret key. When attempting to crack a ciphertext, cryp- 
tographic robustness and complexity of using quantum algo- 


rithms are shown. 


Keywords: cryptographic system, post-quantum crypto- 
graphic algorithm, cryptographic strength, Mann-Whitney U 
test, encryption. 
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XOJJHOTO cbaiina TporpaMMbI HW pa3Mepa CekpeTHOrO KJIrOUAa. 
IIpu TIOMbITKe B3JIOMa mmMdpotekcta TIPOABJIACTCA KDPHIITO- 
rpaduyeckas CTOHKOCTb HW CJIOXKHOCTb MCIOUb30BaHHA KBaH- 


TOBBIX aJIFOPHTMOB. 


KsoueBble c10Ba: KpHITorpapuyeckas CHCTeMa, MOCTKBaH- 
TOBbIM KpHnTorpapwyeckul aIroputM, KpunTocTolkocts, U- 
kputepHii Manna — Yurun, uimdposanne. 


O6paszey Ona YUmUposanua: CpaBHUTeIIbHbIM aHasiM3 MOJH- 
(PHUMpOBaHHOM MOCTKBaHTOBOM KpuMTorpaduyeckow cucTe- 
mbt NTRUENcrypt u oOmjenpuHuatoi Kpuntocuctemp: RSA / 
II. B. Pasymos [ gp.] // Bectnux Jlon. roc. TexH. yH-Ta. — 
2019. — T.19, Ne2. — C.185-194. https://doi.org/ 
10.23947/1992-5980-2019-19-2-185-194 


Introduction. The paper [1] gave an impulse to the development of a new cryptographic system. It shows that 


quantum computers potentially threaten to hack into all widely used cryptographic algorithms. 


The software presented in this paper is cryptosecure versus possible quantum attacks and surpasses its counter- 


parts (for example, RSA cryptosystem) in the algorithm speed characteristics and in the quantity of spendable resources 


[2]. This explains the urgency of the research. 


The object of the study is the NTRUEncrypt cryptosystem. 


The subject of the research is the algorithmic complexity of the development of the NTRUEncrypt cryptosys- 


tem and its modifications. 


The work objectives are to develop an efficient post-quantum cryptographic NTRUEncrypt algorithm, which 


has high cryptographic resistance to quantum computer attacks, to work out modifications of the proposed algorithm, to 


analyze and experimentally validate its advantages. 


To achieve the objective, the following tasks were defined. 


1. To investigate the algorithm of the NTRUEncrypt cryptosystem. 


2. To develop an algorithm for modifying NTRUEncrypt. 


3. To implement the NTRUEncrypt cryptosystem software and its modifications. 


4. To analyze and compare two programs with each other and with their counterpart - RSA cryptosystem. 


Materials and Methods. Consider the description of the encryption NTRUEncrypt system. The cryptographic 
system with the public NTRUEncrypt key uses operations over the Z [X] / (X% — 1) ring of polynomials of degree not 


exceeding N— 1 [3]: 


aA=ao ta, *X1 +a) *X* +--+ ay_y *¥ XN}, 


where dg, Qj, Az... Ay_4 are integers. 


The operations of addition and multiplication are performed as usual, except that X" is replaced by 1, X“*? is 


N+2 
,X 


replaced by X? , 80 on, 


The cryptosystem is determined by a number of parameters, the key parameters are: N, p and q. To maintain 


the algorithm strength, it is necessary for p and g parameters to be coprime. 
To provide high resistance of the algorithm to various attacks, it is recommended to use the following parame- 


ters (Fig. 1): 






































Indication N q p df dg dr Guaranteed resistance 
NTRU167:3 167 128 3 61 20 18 Moderate level of resistance 
NTRU251:3 251 128 3 50 24 16 Standard level of resistance 
NTRUS503:3 503 256 3 216 72 55 Highest level of resistance 
NTRU167:2 167 127 2 45 35 18 Moderate level of resistance 
NTRU251:2 251 127 2 35 35 22 Standard level of resistance 
NTRUS503:2 503 253 2 155 100 65 Highest level of resistance 

















Fig.1. Recommended parameters 
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Research Results 

Key generation. Bob wants to send a message to Alice. For this, he needs public and private keys. Therefore, 
he chooses randomly two small polynomials f and g from the ring of truncated polynomials R. The smallness of poly- 
nomials means that the small polynomial will be much less than g with respect to the arbitrary polynomial modulo gq, in 
which the coefficients are uniformly distributed [4]. To determine the smallness of polynomials, the numbers df and dg 
are used which Bob chooses independently. 

The polynomial f will have df coefficients equal to one, (df— 1) coefficients equal to minus one, and the rest 
coefficients equal to zero. 

The polynomial g will have dg coefficients equal to one, as many coefficients equal to minus one, and the rest 
coefficients equal to zero. 

Bob should keep the selected polynomials in secret since anyone he learns them will be able to decrypt the 
message. 

Further, Bob calculates the inverse polynomials f,, and f, modulo p and q, respectively, such that: 

f X fp =1(mod p) u f x fy = 1(mod q). 
If by chance these inverse polynomials do not exist, then Bob goes back and re-selects the polynomial f 
The secret key is the pair (f, f,), and the public key A is calculated using the formula: 
h=pxX fy x g (mod q). 

Encryption. Alice wants to send a message to Bob using the public key h. To do this, Alice needs to present her 
message as a polynomial m with the coefficients modulo p selected from the range (—p/2, p/2]. Then, Alice needs to 
choose another small polynomial 7 which is called “blinding”, and calculate the ciphertext: 

e = (rxh+m)(mod q). 

Decryption. Bob receives an encrypted message e from Alice and wants to decrypt it. First, using his secret 

key, Bob calculates: 
a=f Xe(mod q). 
Since Bob calculates the value a modulo q number, he should choose its coefficients from the range (—q/2, q/2] 
and then calculate: 
b=a(mod p). 
Finally, Bob, using the second part of the secret key, receives the original message from Alice: 
c = fp, X b (mod p). 

Modification of the NTRUEncrypt algorithm. As can be seen from the description of the algorithm, the poly- 
nomial f shall comply with the following requirements: 

- the polynomial fis invertible modulo p, 

- the polynomial fis invertible modulo q, 

- the polynomial fis a small polynomial. 

In the algorithm itself, the invertibility modulo p and g was guaranteed as follows. If the polynomial f being not 
invertible in one of the moduli was generated, then it was discarded and the next one was generated — and so on as 
long as the required polynomial was found. 

The proposed modification is to replace the polynomial fby a polynomial of the form: 

f=1+pF, (1) 
where F is a small polynomial. 


This approach has the following advantages. 
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1. From the expression (1), it is clear that the polynomial fis always invertible modulo p. This fact accelerates 
the key generation since it is not necessary to additionally calculate /,. 

2. Since f-' = 1 mod p, then decoding does not require additional multiplication by f~', which speeds up the de- 
cryption process itself. In this case, the private key will not be the pair (f, f,), but (f). 

Key generation. As in the original algorithm, Bob first selects the encryption parameters N, p, q and the num- 
bers df, dg. Then, he selects randomly two small polynomials F and g from the ring of truncated polynomials R. 

He calculates the modified polynomial fusing the formula (1). 

Next, Bob calculates the inverse polynomial f, modulo q: 

f X fy = 1(mod q). 
If by chance an inverse polynomial is not found, Bob goes back and re-selects the polynomial /- 
The secret key is the polynomial f, and the public key / is calculated as follows: 
h=pxX fy X g (mod q). 

Encryption. Encryption remains unchanged; all is quite as in the original NTRUEncrypt algorithm. 

Alice wants to send a message to Bob using the public key 4. To do this, Alice needs to present her message as 
the polynomial m with coefficients modulo p selected from the range (—p/2, p/2]. Then, Alice needs to choose another 
small polynomial r, which is called “blinding”, and calculate the ciphertext: 

e = (rxh+m)(mod q). 

Decryption. Bob receives the encrypted message e from Alice and wants to decrypt it. In the first place, using 

his secret key, Bob calculates: 
a=f Xe(mod q). 

Since Bob calculates the value a modulo qg number, he should choose its coefficients from the range (—q/2, 

q/2], and then calculate: 
b=a(mod p). 

That is all, calculation is finished at that; we have received the original message from Alice: b = m [5]. 

Proof of the modified algorithm. To prove the algorithm, consider the decryption process itself. 

Alice’s encrypted message looks like: 

e = (rxh+m)(mod q). 
Bob uses his secret key — polynomial f: 
a=fxe(modq) =(f X (r xX h + m))(modq) =(f X (r X pfy X g + +m))(mod q). 
As a result: 
a=(pr X g +m X f)(mod q). 
Thereafter, Bob obtains the polynomial b through decreasing the coefficients of the polynomial a modulo p: 
b = a(modp) = m x f(modp) = (m + m X p X F)(mod p) = m (mod p). 

Thus, we have checked and proved that the polynomial 5 is the original message m indeed. 

Implementation of the algorithm. The programming language used is the object-oriented programming lan- 
guage C#, which belongs to the family of languages with a C-like syntax. Development framework was Microsoft Vis- 
ual Studio 2015 Enterprise. The primary advantage of this software product is the application with a graphical interface 
that allows the user to make short work with the device and the operation scheme of this software product. 


Fig. 2 shows a generalized block diagram of the software. 
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Fig. 2. Flowchart of software implementation 


Example of the software. NTRUEncrypt uses three constant parameters: N, p, g. The user enters them in the 
System Settings panel. 

Next, there is the Key Generation panel. The parameters entered by the user in the previous step are input. 
Then, two polynomials fand g are randomly formed, so that: 

- the number of coefficients equal to one and zero to be equal to the number that was previously determined in 
the program; 189 
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- the degree of polynomials to be corresponded to the input parameter N. 

Then, the Euclidean theorem on computing the greatest common divisor for polynomials and its backward 
course is used to calculate two polynomials inverse to fmodulo p and gq, respectively. The public key h is computed. 

The user enters the binary message type m in the Encryption panel in the Message window. It should be men- 
tioned that the original message should be divided into blocks of N bits, each of which will be processed separately and 


converted into a polynomial with the coefficients {—1, 


0, 1} (in this program, instead of the value “—1”, “2” is used). 


For an encryption operation, the program needs to perform preliminary preparation — to generate one blinding 
polynomial. It is formed in the same way as f and g polynomials. 

The resulting public key and the polynomial r make it possible to encrypt message m using the appropriate 
formula. Then, the result is checked on the membership to the ring of truncated polynomials of a degree not exceeding 


N-—1, and is displayed on the screen. 


The next block implements the decryption mechanism. For this, actions are consistently performed considering 


the ring ZLX)/(X% — 1). 
The operation result is shown in Fig. 3. 


= NTRUEncrypt 





Key generation 


Cryptosystem . 
parameters Key generation Encryption 
4 
. | Common message: 
Ne iW | dia 4,dg=3 Be: | 10011101001 
~ “7 
-“.- ' Generate £ = (xO) 4x 1ox"2ox"3+{x"6) +(0° 10) 
qe {7 Encrypt 
f 4x" © (37D) © 63°) © 6305) 07 
Q =x Dela” 3) (73) ox “Sof °S) ox 10 ; 
Decryption 
Clear ; 
fq = 333-3-3230233 Ciphertext 210-1-31-211-31 
Output h= 1102-3312-202 


Public key: (h); 
Private key; (f), 


Fig. 3. Appearance of NTRUEncrypt program 


Fig. 4 shows the software results in the form 


@: text.tt —Broxnoy fi shifr_text.t¢ 
®aiin Mipaska Dopmart| Pain Mpasxa 
Texct 3To: nucbmeH|’Ab p Pu EAT ,¢ 


BbICKa3bIBaHMA un ci°AB&p MuB EAB p 
peanu3aumn; Boipaxen| €A.. t MAB] A..,d 0 
3Hakos (He3zasucumo [M@l Ab ~ uB EAL 
peanu3aumn) yyecte) IIAT"*¥ Mu EP... 
T.4. autepatypxoro [°Ab&Y Ev@l: ,,(d 
KommMyHuKaumM, o6na, EAP d YaB €P....p 
(yenoctHocTEw) mo} MBIEPb"dYIMu@ Ab 
Texct Bo BTopom 3H\Du@l A. ~YIEWBIy 
B TpeTbem 3Ha4eHun @u@ yA..(t¥ 

oTNpaBHow TouKon B/BIEP..$p AB pr" 
ycTaHosnexve Ha ocil Af"> 0) wBIEP...,~ 
oTHOweHHii Mexay er) $Y LuBIEA] DY Ea 
Hepenko nog texctoi€Pb$d 6 Ul ,,"d¥ 
OT OTQZeNbHOM mopdel P....p (Av@lEA > YIN 
ynotpe6nexue paamni| W@IEA,,.~ Ev@leA.. 
uccnegosatensckux | (br.d¥ En@le *t 









eS Texct 8 tTpetbem 3H\€P..$d¥ w@l Ar. 
= «peyesoe nponssene| W@lyYAb p H BIEPb 
77) «Teopuw uWenoro Tex)Qa IAb qd 

g ceoux acted, un (\MBIIP,, p¥ J 
= o6nagaet ussectHon A... » YIE@BIVA..."py¥ 
ws cneuguyecko oprajl] AT" L BIIPh$d 
™= ‘opraHu3zaynei atux |Ew@: T"t¥l uw OIA, 
s BHYyTpeHHei cTpyKTy|/P [Pa@ A,(d a 
o 

S Fi 
= 
f=y 

— 

a= 


of source text, encrypted and decrypted. 


te) text.tt — Bnoxnot = o x 
®opm @Maiin Mpaska Dopmat Bug Cnpaska 


TekcT STO: MMCbMeHHaA MAM NeYaTHAaA oukcauna a 
peyesoro BbICKa3bIBaHMA MNM COObUeEHHA B 
NPOTMBONONOXHOCTb yCTHOK peanu3aunn ; BoipaxeHHaa “ 


BakpenneHHaA NOCPEACTBOM ASBbIKOBbIX 3HAKOB 


} (He3aBMCMMO OT NMCbMeHHOM MAM YCTHOM dopmbl Hx 


peanv3zauMv) YyBCTBEHHO BOCNPMHMMaemMaA CTOpOHa 
peyesoro, B T.4. AMTepaTypHOrO, npowsBemeHHA; 
MunwManbHaA eQMHuua peyesod KOMMyHMKaLIMM, 
o6nafawujan OTHOCUTeNbHEIM eQMHCTBOM 
(ueNOCTHOCTBH) MW OTHOCMTeENbHOM aBTOHOMMeN. 

TekcT BO BTOpOM 3Ha4eHMH, ABAAWWYMUCA OAHWM M3 
acnekToB TekcTa B TpeTbeM 3HaYeHMM (a MMEHHO, 
ero «nmaHoM BbIPaxeHMAY), CAYXMT OTNPAaBHOM TOYKOK 
Bcex Ownonormyeckux Npouegyp. Ero ycTaHoBneHve 
Ha OCHOBe KPMTMYeCKOrO M3YYeHMA MCTOPMM TeKCTa M 
OTHOWeHW Mexay ero BapwaHTamu CoCcTaBnAeT B3aga4y 
Tekctonoruu. Hepeako nog TeKCTOM NOHMMaeTCA 
Takxe nw60H peyesow oTpe30K — OT OTZeNbHOK 
Mopdems! M0 AUTepaTypbi B WeNOM, OfHAKO Takoe 
ynotpe6neHne pasmbiBaeT CneuMduKy NOHATHA MW BeneT 
K HeYeTKOCTH uccnegopatenbcKhx pesynbTaTos. 
TekcT 8 TpeTbem 3Ha4eHMM («COObWeHKe», «YenbIi 
TekcT», «peyesoe npoussegeHve») — npegmeT 
oco6oro pa3fena ounonormu, «Teopwuw Yenoro 
TekcTa», WAM AMHFBMCTMKM TeKCcTa. B oTau4ne OT 


4. Software result 


Experiments. The experimental research objective is to show and prove the advantages in the performance 


characteristics of the developed cryptosystems in comparison with the analogue. (In this paper, an RSA public key cryp- 
tosystem was chosen as an alternative [6].) In addition, it was necessary to compare the implemented NTRUEncrypt 


algorithm and its modification in order to note the advantages of the latter and provide experimental evidence. 
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To verify the results, the statistical Mann — Whitney U-test was used, which is intended for comparing two in- 
dependent samples by the level of any feature measured quantitatively. The criterion enables to determine the degree of 
difference between samples and is more powerful than the Rosenbaum criterion. 

For the experiment, a third-party software implementation of the RSA cryptographic system was used. An el- 
ement of the Stopwatch class was introduced in the source code of all three programs — NTRUEncrypt, RSA, and 
NTRUEncrypt modifications. This class provides a set of methods and properties that can be used to accurately measure 
the time spent on the execution. Thus, it became possible to record the results of the elapsed time at all three main stag- 
es: key creation, encryption and decryption of the message. The experiment results are shown in Table 1. 

Table 1 
Experiment results 








































































































Ne NTRUEncrypt operate time, s cea ne raas RSA operate time, s 
Key generation 
1 0.0239676 0.0190014 0.0964144 
2 0.0149743 0.0109994 0.1023557 
3 0.0129915 0.0099936 0.1372658 
4 0.0170503 0.0099772 0.0491137 
5 0.0139885 0.0099773 0.0869555 
6 0.0139880 0.0099782 0.0587503 
7 0.0129761 0.0099936 0.0986608 
8 0.0139931 0.0109775 0.0707417 
9 0.0139918 0.0099773 0.0595015 
10 0.0169748 0.0099777 0.0517874 
Encryption 
1 0.0069961 0.0069770 0.0453772 
2 0.0049961 0.0059796 0.0598658 
3 0.0059954 0.0049965 0.0265461 
4 0.0049966 0.0069966 0.0402407 
5 0.0069765 0.0060854 0.0783415 
6 0.0049961 0.0069989 0.0097463 
7 0.0059795 0.0059954 0.0247876 
8 0.0069961 0.0059954 0.0272320 
9 0.0089956 0.0049798 0.0494907 
10 0.0059954 0.0049744 0.0100916 
Decryption 

1 0.0009843 0.0009843 1.0624476 
2 0.0010212 0.0010026 0.3692022 
3 0.0010193 0.0009974 0.765 1423 
4 0.0009989 0.0009975 0.9922103 
5 0.0030012 0.0009988 0.3056423 
6 0.0010017 0.0009975 0.3757183 
7 0.0009989 0.0009988 1.2844272 
8 0.0010021 0.0009989 0.5902735 
9 0.0010026 0.0009984 0.9932694 
10 0.0009844 0.0009970 0.8740392 
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The Mann — Whitney U-test is applicable to the results given in Table 1. Compare the key generation rate of 
the NTRUEncrypt algorithm and its modifications. 

Following the algorithm of the Mann-Whitney U-test, we criterion write the computation step by step. 

1. Create a common ranked list of both samples assigning a lower value to a lower rank. 

2. Divide the total ranked list into two consisting of the elements of the first and second samples. 

3. Calculate the sum of ranks for the first and second samples separately, as shown in Fig. 5. 


a a 


0.0129915 0.0099936 
0.0170503 0.0099772 











0.0139885 0.0099773 


 _—.— iS aa 
7 a 
0 ne9r08 a rs 2 J 


a 


Fig. 5. The second step of calculating the Mann-Whitney U-test 











4. Calculate the value of the Mann-Whitney U-test: U= 9. The critical value of the criterion for the data n1 and 
n2 should be determined from to the table of statistical significance level (Fig. 6). 





Fig. 6. The forth step of calculating Mann-Whitney U-test 


Since the values of m1 and n2 are the same for all experiments, this table will be used in each calculation. 

Hence it follows that the obtained empirical value U = 9 is in the zone of significance (Fig. 7). Consequently, 
there is a significant difference between the speed of the NTRUEncrypt program and its modification. If this fact is ex- 
pressed percentagewise, it turns out that the NTRUEncrypt modification is 28% faster than the program itself. 


Axis of significance: 









U, 
a Zone of 
insignificance 





19 27 


Fig. 7. Axis of significance 


Similarly, in shorthand, we present the results of the Mann-Whitney U-test and the percentage superiority for 
the remaining cases. 

Key generation 

1. NTRUEncrypt modification is 28% faster than NTRUEncrypt. The obtained empirical value U = 9 is in the 
significance zone. 

2. NTRUEncrypt program is faster than RSA by 80%. The obtained empirical value U = 0 is in the area of sig- 
nificance. 

3. NTRUEncrypt modification is faster than RSA by 86%. The obtained empirical value U = 0 is in the area of 
significance. 

Message encryption 

1. The obtained empirical value U = 47 is in the zone of insignificance. 
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2. NTRUEncrypt program is f80%aster than RSA. The obtained empirical value U = 0 is in the area of signifi- 
cance. 
3. NTRUEncrypt modification is faster than RSA by 84%. The obtained empirical value U = 0 is in the area of 
significance. 
Decryption 
1. NTRUEncrypt modification is 17% faster than NTRUEncrypt. The obtained empirical value U = 24 is in the 
zone of uncertainty. 
2. NTRUEncrypt program is faster than RSA by 99%. The obtained empirical value U = 0 is in the area of sig- 
nificance. 
3. NTRUEncrypt is99% faster than RSA. The obtained empirical value U = 0 is in the area of significance. 
Discussion and Conclusions. In the framework of this research, the following is developed: 








- software that implements the NTRUEncrypt cryptosystem operation; 

- software that implements this cryptosystem modification. 

One of the NTRUEncrypt advantages over its counterpart - the RSA cryptosystem, a higher speed of operation 
can be specified. Performing encryption and decryption operations requires O(n’) operations, unlike O(n*) in the same 
RSA. As for the experimental data, the NTRUEncrypt program wins significantly on the speed of the algorithm com- 
pared to RSA. In addition, there is a slight increase in durability with the same key length. The disadvantage of the sys- 
tem is the necessity for using the recommended parameters. 

Regarding the NTRUEncrypt resistance, after the creation of quantum computers, the problems of fast factori- 
zation and discrete logarithmation will be solved [7]. In this case, RSA, DSA, and similar algorithms will become use- 
less. The relevance of NTRUEncrypt will remain: it will be fully applicable in the “post-quantum” era since there is no 
algorithm that solves the problem of the shortest lattice vector. 

The advantage of using the NTRUEncrypt algorithm modification is experimentally proven. The developed 
application performs the general work on key generation, encryption and decryption 25% faster. In addition, it optimiz- 
es the use of internal memory through reducing the weight of the original program file and the size of the private key. 
When attempting to crack a ciphertext, cryptographic robustness and the complexity of using quantum algorithms are 
manifested. 
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> Voxesckuii TOproBo-2KOHOMH4eCKHii TeXHUKyM, Pr. VoKesck, Poccuiickaa Penepanus 


Introduction. The solution to the problem of optimal control of 
the biogas process under its conversion in two digesters is con- 
sidered. The work objectives are to propose a mathematical 
model of this process and to develop an optimal control algo- 
rithm. 

Materials and Methods. The developed mathematical model 
describes the biomethanation from animal waste through the 
downstream processing of the substrate in two digesters. Cases 
of the same and different temperature media (mesophilic and 
thermophilic) are considered. An optimal control problem is 
defined as a Lagrange problem for this model. Its modifiers are 
the rates of substrate entry into the digesters. The algorithm for 
solving this problem is proposed; it is based on the numerical 
implementation of the Pontryagin maximum principle. When 
optimizing, a hybrid genetic algorithm was used with an addi- 
tional search in the neighborhood of the best solution through 
the conjugate gradient method. 

Research Results. A new mathematical model is developed. It 
describes the biomethanation during the downstream processing 
of the substrate in two digesters. A numerical algorithm for 
solving an optimal control problem is proposed and software- 
implemented. The numerical studies have shown that the biogas 
production rate is nearly twice as high for a thermophilic medi- 
um as for a mesophilic one. It is established that the down- 
stream processing of the substrate in two digesters with the 
same temperature medium allows the biogas production rate to 
be doubled. If the temperature media in the digesters are differ- 
ent, then in the first of them, the mesophilic medium should be 
used, and in the second - the thermophilic medium. At this, the 
biogas formation rate is somewhat lower compared to the case 
when there is a mesophilic medium in each of the digesters; 
however, the degree of the substrate processing is by 10-15% 
higher. 


“ The research is done within the frame of the independent R&D. 
* B-mail: stkj@mail.ru, MaykovD@yandex.ru 
*“* Pa6ota BBITIOJHeHa B pamkax uHnMaTuBHon HUMP. 


Beedenue. CraTba MocBaljeHa pelleHHio 3aa4un ONTHMAaJIbHO- 
TO ylipaBlleHHua TIpoleccoM mowyyenua Ovora3a pH Herpe- 
PbIBHOM pexkHMe ero MepepaOoTKH B AByx MeTaHTeHKax. Lem 
paOorTsl: npeycTaBHTb MaTeMaTHYeCKylo MOJeIb JaHHOrO TIpo- 
yecca HM pa3paOoTaTb aIropuTM BbIOopa ONTHMasIbHOrO yupaB- 
JIeHHA. 

Mamepuaavi u memooovi. Co3qaHHad MaTeMaTH4ecKkad MOjeJIb 
OMMcbIBaeT TOyyeHHe Ouora3a H3 OTXOOB %KHBOTHOBOJCTBA 
IIpH ocdeqOBaTeNbHON MepepadboTKe cyOcTpaTa B WBYX Me- 
TaHTeHKax. PaccMaTpuBaloTca Cilyya OAMHAKOBBIX VM pa3yIH4y- 
HBIX TeMMepaTypHBIX cpey (Me30PuJIbHOH uM TepModuIbHO!). 
J\na qaHHol Moye ctopMysMpoBaHa 3afa4a ONTHMasIbHOTO 
ynipaBieHHa B Bue 3aqauH Jlarpamxa. Ee ynupaBiaroljuMu 
TlapaMeTpaMH ABIAIOTCA CKOPOCTH HOcTyleHuA cyOcTpata B 
MeTaHTeHKH. IIpeqnoxeH anropuTM pellieHua JaHHol 3aqa4H, 
OCHOBaHHBIM Ha 4HC 

JIeCHHOM peasM3alHu UpHHyuna MaxcuMyma Iloutparuua. IIpu 
ONTHMH3alHH MpHMeHsAsicd THOpHAHBIM TeHeTHYeCcKHH ayro- 
PHTM C JOMOJHHTeIbHEIM MOMCKOM B OKPeCTHOCTH JIy4Lero 
pelieHHaA MeTOJOM COMpsmKeHHBIX rpay{HeHTos. 

Pesyismamei uccaedosanua. PaspadoTaHa HOBad MaTeMaTH4e- 
cKad MOJesIb, ONMCbIBaIOMad Mpolecc Nomyyenua Onora3a pu 
TIOCHIeOBATeMbHOH MepepadorKe cyOcTpata B ABYX MeTaHTeH- 
kax. IIpequioxkeH HM MporpaMMHO peasM30BaH YMCICHHBIM all- 
TOpHTM pelleHHa 3afa4H ONTHMaIbHOTO = ylipaBue- 
Husa.4HcHIeHHbIe MCCeqOBaHHA MoKa3aiu, 4TO WIA TepMo- 
(UIBHOK CpesbI CKOpOcTb OOpa30BaHua Ovora3za WpakTH4eckH 
BBOe BbIle, YeEM Wd Me30uIbHOM. YcTaHoOBIeHO, 4TO T0- 
cileqOBaTebHaad Tepepadotka cyOcTpata B JByX MeTaHTeHKax 
C OJMHaKOBLIMH TeMIepaTypHBIMU CpeaMu MO3BOsAeT BBO 
YBeJIMYHTb CKOpOcTb oOpa30BaHHa Ovora3a. Ecuu Temmepa- 
TYPHbIe CpeybI B MeTAHTeHKaX pa3JIM4HbI, TO B TepBOM 3 HHX 
CileyeT UCMOUb30BaTb Me30@HJIbHyI0 Cpezy, a BO BTOPOM — 
TepModunbyyro. [pu stom cKopoctb oOpa30BaHHa Ouora3a 
HeCKOJIBKO HWKe 10 CpaBHeHHIO CO CylyyaeM, Kora B Kaxk OM 
M3 Me€TaHTCHKOB Me30QHJIbHad Cpefa, OHAaKO CTereHb Mepe- 
paOorku cy6ctpata Bbrme Ha 10-15 %. 
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Discussion and Conclusions. The results obtained can be used 
for the calculation and design of biogas plants, as well as in the 


development of appropriate software. 
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O6cyarcdenue u 3axmouenua. TlomyueHHble pe3ysIbTaTbI MOTyT 
OBITb HCHONb3OBaHEI DJIA pacueTa H KOHCTpyupoBaHHua Onora- 
30BbIX YCTaHOBOK, a TakoxKe Ip paspaOoTKe COOTBeTCTBYIOLIe- 


To IporpaMMHoro oOecneueHHa. 


KsoueBble c10Ba: MeTaHoreHe3, Ov0ra3, MeTaHTeHK, Mepepa- 
OoTKa OTXOJOB 2%KHBOTHOBOJCTBa, MaTeMaTH4eCKad MOJ[eJIb, 
cucTtema yuddepenuHabHbIX ypaBHeHHi, 4MCeHHOe pele- 
HHe, MeTOJbI ONTHMH3ALHH, ONTHMAIbHOe yiipaBsleHHe, TIpHH- 
yu Makcumyma Ilontpsruna. 


Oopazeu Ona yUumupoeanua: Kopones,C. A. Onrumu3anna 
TBYXCTaMMHOTO pexkHMa MeTaHoreHe3a Ha OCHOBe MIpHHIHIa MaK- 
cumyMa Ilontpsruna/ C. A. Koposes, J]. B. Matixos // Becrauk 
How. roc. tex. yu-ta. — 2019. — T.19, No2. — 
C.195-203. https://do1.org/10.23947/1992-5980-2019-19-2-195- 


Introduction. When livestock enterprises, in particular, poultry farms, pig farms, cattle farms, are operated, a 
large amount of waste is generated. As a result of their fermentation in special tanks (digesters), fuel gas (biogas) and 
valuable organic fertilizers can be obtained [1]. This process is called methanogenesis. It can occur under the periodic 
and continuous delivery of the substrate. In the first case, the digester is filled once and emptied completely upon com- 
pletion of the fermentation. In the second case, two processes occur simultaneously and continuously: the supply of the 
substrate and the removal of its processed portion. 

Usually, mesophilic medium (at a temperature of 25—38° C) or thermophilic (45—60° C) is used for methano- 
genesis. The optimum temperature for the mesophilic medium is 37° C, for the thermophilic medium, it is 56° C. The 
fermentation time for these media is 25 and 12 days, respectively. 

The economic efficiency of biogas production depends on various factors: type and quantity of raw materials, 
climatic conditions [2], etc. In addition, the rate of substrate input into the digester affects significantly the biomethana- 
tion. This parameter value depends on the volume of the digester and the type of raw materials. To find the optimal val- 
ue of the specified magnitude, it is required to solve an optimal control problem. 

Under the continuous mode of fermentation, the substrate does not have time to go through full processing. To 
increase the biogas production, it is necessary to use two digesters so that the substrate is sequentially processed in each 
of them. Various aspects of this process are studied in [3, 4], and its technical implementations are reflected in the pa- 
tents [5—7]. 

In the papers on the mathematical simulation of methanogenesis, for example [8, 9], the search for optimal 
control is not presented or its asymptotic value is found [10—12]. An analytical solution to the optimal control problem 
is obtained in [13]. However, the mathematical model described there differs significantly from that presented in this 
paper. The numerical method proposed here for solving an optimal control problem is applicable to a wide class of 
models. 

Materials and Methods. The scheme of downstream processing of the substrate is shown in Fig. 1. 


Digester 1 Digester 2 


q\ 
X,, L, X,, L, 


Fig. 1. Scheme of two-stage mode of methanogenesis 





Ly, 
fo Fo, q> 











Assume Lo is the concentration of nutrients in the substrate fed to the first digester (kg/m*); L, and Ly are the 
concentration of nutrients in the substrate for the first and second digesters, respectively (kg/m*); X. , and X> are the 
concentration of methane-producing bacteria in the first and second digesters (kg/m*); Q, is the volume of the substrate 
in the first digester (m’), Q,- in the second one; qo is the rate of the substrate input in the first digester (m*/day), 
dQ dQy 

dt 


q, = — and qo = are the rate of substrate delivery from the first and second digesters, respectively. 


Losses of the substrate do not occur, therefore 
90 = = 4942 =4- 
The value of qq is determined by the volume of waste, the quantity and capacity of the biogas plants. 
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Relative rates of substrate input for the digesters are: 
1 d 
= AA and 


i= 


= uy = A 
QO dt Q Q, dt Q 


Substrate volumes in the digesters are equal to 


q q _mQ 
QO, =— and Q, ==. 
uy uy uy 


_ 1 dQ, _ 4 


To describe methanogenesis during downstream processing of the substrate, a mathematical model based on 
the mathematical model of the population dynamics of methane-producing bacteria for a one-stage methanogenesis 
mode is used [11, 14]: 





baa Mindi a} x, 











dt ay + Ly by + a 
dL Bibmei x 
Sram (Lo ii) 1Mmg141“*1 
t aq + a (1) 
IX. L b 
dXy = hain X, + Hmg242 mars Bie, 
dt ay + Ly by + L, 
dL, Bobmg2Lo Xp 
ae -(L,-Ly) 
t ay + Ly 


Here, the subscripts of the variables ( X;, L; ) and parameters (Um j>bmd i> 4» 5;, B;, uj) correspond to the 
number of the digester? € {1, 2}. The model parameters arep1,,, ; and [,,q ;, the maximum possible relative rates of 


growth and dieoff of bacteria, respectively (day '); B;is dimensionless coefficient of the substrate absorption; a; and 
b; are empirical coefficients (m*/kg); 4. is parameter equal to zero if the temperature media in the digesters are different, 
and equal to one if these media are the same. 

Values of the model parameters are determined in accordance with the selected temperature regime of the 
methanogenesis (mesophilic or thermophilic). If temperature media in the digesters are the same, 
then Unel =Hmg2 =Emg > Ud = Uimd2 = Hind > ete. 

The initial conditions are: 

X (0) = Xp (0)= Xp, L,(0) =Ly (0) =Lp, 
where XQ is a natural concentration of methane-producing bacteria in the feedstock; Lp is equal to the concentration of 
nutrients in the unprocessed substrate. 
The model (1) is built on the assumption that the digesters maintain an optimal and constant process temperature. 
If the temperature media are different, then, when the substrate is fed from one digester to another, heating or cooling to 
the temperature in the second digester should occur. 
The rate of biogas production (m*/day) in the i-th digester is equal to 
Yi Ung i LX, i 
w; = ———— 
a; + L; 
where y; is the coefficient characterizing the rate of conversion of substrate nutrients into biogas (m? « m’/kg). 
To obtain the optimal control problem, it is necessary to supplement the system of equations (1) with the criterion 
functional 





(a Y2Mng2L2X2 
V= } + 


Ja >max, (2) 
aA + Lh, ad + L, 


0 
that determine the total biogas yield from 1 m? of the substrate in the first and second digesters over time T . The opti- 
mized parameters of the problem are relative rates of the substrate delivery into the digesters uw, and win the system of 
equations (1). 

The systems (1) - (2) are Lagrange problems. In general case of the optimal control problem, there is a system of 
differential equations of the form: 


X =t(x(1), u(t), 1), re [0, T],xeR",ueRr*, (3) 


its initial conditions are: 
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x (0) =Xq- 
It is required to find the optimal control u (t) delivering a maximum to the criterion functional 
T 
J =| Fo(x(¢), u(t), t) > max. (4) 
0 
For the numerical solution to the problem, a difference grid with nodesfg=0, 4, bh... tj, bays oo 


t, =T with the constant step / =t;,, —¢; is introduced on the interval [0, 7] [15]. 


The numerical solution to the system of differential equations (3) is carried out through the fourth-order 
Runge-Kutta method: 


Xi4] =x; + 2-(ky +2ko +2k3 +k4),1=0,q-1, 


ky = f(x;,, Uu;, ti), 





h h 
ko= f(x 5k, u;, n+), (5) 
h 
2 


h h h 





To solve the problem, it is convenient to introduce gk -dimensional full control vector U = (u,), i=1,q. In 


this case, the difference approximation of the criterion functional (4) is the expression: 
q 
J =J(U)= > Fo(x;, uj.t;):2 > max. (6) 
i=0 


The problem is solved in the sequence given below. 

1) Set the full control vector U . 

2) The initial system of differential equations (3) is solved numerically using the relations (5), and the value of 
the criterion functional (4) is calculated using the difference approximation (6). 

3) The system of adjoint equations is numerically integrated (in the direction “from right to left”) according to 


the relations: 





OF : 
Pi = Pit th +G} “Pst, i=l, q-1, (7) 
i 


Py =0. 


Of (x;,U,.t; 


Here, p = p(t) are dual variables of the Pontryagin maximum principle; G; -( 3 is Jacobi matrix com- 


i 
posed for the system (3). 

4) The optimization process is performed by vector U . In this paper, the genetic algorithm with real coding 
and an additional search in the neighborhood of the best solution through the conjugate gradient method was used as the 


optimization method. 
For the considered problem of methanogenesis optimization, the vector of phase variables is equal 


to x = colon(.X;, L;, Xz, Ly), the control vector is equal to u = colon (1m ,u> ), and the Jacobi matrix has the form: 




















Hmgtlt bait i Mmg *1 4 Pied *1 0 0 
ath b+ (qe) (oeely 
Bibel Bia mei X1 r . 
_ PI mg go 
G ay +1 (a +L, a 
» 0 Hmg242 mary AM mg2X 2 bobtndrX2 
uy 2 5) + 7 
ay +Ly by + Ly (a, +L,) (b, +L,) 
' , Bobing2L2 , B249Hmg2X2 
2 et 2 
ay +Ly (ay +L, : 
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Research Results. The theoretical values of the variables found as a result of the numerical solution to the sys- 
tem of differential equations of the mathematical model of methanogenesis are recorded. The deviations of the specified 
values from the experimental values with respect to the vector of parameters [10, 16] (Table 1) are considered. The pa- 
rameter values of the methanogenesis model are estimated through minimizing the sum of squares of these deviations. 


cide. 


























Table 1 
Estimates of methanogenesis model parameters for mesophilic / thermophilic media . 
Source of raw materials 
Coefficient 
Poultry factories Pig farms Cattle farms 

Ling 0.482 / 0.821 0.346 / 0.783 0.297 / 0.563 
Und 0.353 / 0.528 0.291 / 0.423 0.254/0.351 

a 34.781 / 43.875 7.242 / 21.653 5.013 / 8.733 

b 116.457 / 14.674 37.347 / 9.278 18.722 / 5.455 

B 2.344 / 3.189 1.495 / 2.084 1.413 / 1.983 

Y 1.463 / 1.963 1.373 /1.907 1.299 / 1.813 




















The following methods for processing the substrate are considered. 
I. A single digester is used. 


II. Two digesters are used with the downstream processing of the substrate, in which temperature media coin- 


At this, two options of the medium are possible for each digester: mesophilic or thermophilic (Table 2). 
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Table 2 



































Gineciies Source of raw materials 
Poultry factories Pig farms Cattle farms 
I. One digester with mesophilic / thermophilic medium 
u,day | 0.149 / 0.402 0.128 /0.341 0.112 / 0.268 
w ,m?/ day 7.25 / 24.51 3.58 / 10.76 1.81 /5.42 
L'/ 1p, % 56/43 51/38 52/41 
II. Two digesters with mesophilic / thermophilic medium in each 
us day | 0.168 / 0.432 0.144 / 0.362 0.117 / 0.287 
U5 , day” 0.491 / 2.448 0.412 / 1.858 0.367 / 1.424 
w ,m?/ day 15.19 / 50.21 7.79 / 22.19 3.85 / 10.92 
Dll, % 45/35 38 / 33 38 / 33 














Il. Two digesters, in the fir 


st - mesophilic, in the second - thermophilic medium / in the first - thermophilic, 
in the second - mesophilic medium 


























uy , day | 0.149 / 0.402 0.128 / 0.341 0.112 / 0.268 
u>, day! 0.312/0.051 0.222 / 0.059 0.183 / 0.051 
w", m?/ day 15.35 /25.12 6.61 / 11.42 3.26 / 5.42 
Lol ig% 30/34 29/35 26 / 32 








Here, u are the optimal values of the relative rate of the substrate delivery; w is the corresponding daily aver- 


* 


age biogas yield; a is the completeness of the processing of substrate nutrients. 
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Graphs of the biogas production rate are shown in Fig. 2. 
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Biogas production rate: 
in scheme with one digester 
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Fig. 2. Biogas production rate when using waste of pig farms for two digesters with mesophilic (a) 
and thermophilic (6) medium in each 
IH. Two digesters are used with downstream processing of the substrate, in which the temperature media are 


different. 
Here, two cases are considered. 
1) In the first digester, there is mesophilic medium, and in the second - thermophilic medium. 
2) In the first digester, there is thermophilic medium, and in the second - mesophilic one. 
See the optimal parameter values in Table 2. 


Fig. 3 shows a graph of the biogas production rate. 
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Fig. 3. Biogas production rate when using pig farm waste for two digesters: in the first digester there is mesophilic medium, 
in the second digester - thermophilic medium (a); in the first digester, there is thermophilic medium, in the second digester — 
mesophilic medium (b) 


The use of a single digester in the continuous mode of fermentation enables to process the substrate nutrients 
by 40-45% for the mesophilic medium and by 50-55% for the thermophilic one. 

When using two digesters with the same temperature media, the share of substrate processing increases by 10% 
for the mesophilic and by 5% for the thermophilic media. At this, the average daily production of biogas is almost dou- 
bled. The biogas production rate for the thermophilic environment is almost twice as high as for the mesophilic one. 

Two digesters operating in different temperature media show different results. So, in the first case, the share of 
substrate processing is about 25% higher than when using a single digester with a mesophilic medium. In the second 
case, the share of substrate processing is 5—10% higher than when using a single digester with a thermophilic medium. 
Besides, the organization of the process in the first case provides a nearly twofold increase in the rate of biogas produc- 
tion (when compared to a single digester with a mesophilic medium). In the second case, a significant increase in speed 
is not observed (if compared to a thermophilic medium). It follows from the above that the second case is not advisable. 

The first option loses somewhat in the rate of biogas production in case of using two digesters with a meso- 
philic medium in each, although it provides a higher (by 10—15%) degree of the substrate processing. This is due to the 
fact that with the use of two digesters with a mesophilic medium, the optimal values of the relative (and therefore abso- 
lute) rates of substrate input are higher, that is, more of the substrate is processed per unit of time. The rate of biogas 
generation for the second option is twice as low as for the case of using two digesters with a thermophilic medium. 

Discussion and Conclusions. Mathematical models are developed that describe the methanogenesis during the 
downstream substrate processing in two digesters when the temperature media in them coincide and differ. These mod- 
els correspond to the problem of optimal control of the methanogenesis process on the basis of the Pontryagin maxi- 
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mum principle. The paper describes an algorithm for its solution. The control parameters are relative rates of the sub- 
strate delivery into the digesters. 

The numerical study results show that the sequential application of two digesters with the same temperature 
media allows for an increase in the degree of substrate processing by 5—10%. At the same time, the rate of biogas pro- 
duction is doubled. The degree of processing of the substrate for the mesophilic mode is also S—10% higher than for the 
thermophilic one. At the same time, the biogas production rate in a thermophilic medium is almost twice as high as in a 
mesophilic one. This is due to the higher intensity of the process (higher than the optimal value of the relative rate of the 
substrate delivery). 

The results of the operation of digesters with different temperature media are also shown. In the first case, the 
substrate is first processed in a mesophilic medium, and then enters the digester with a thermophilic medium. In the 
second case, on the contrary, the substrate from the digester with a thermophilic medium enters the digester with a mes- 
ophilic medium. In the first case, more substrate is processed. When compared to the use of a single digester with a 
mesophilic medium, the advantage is about 25%. When compared to the use of a single digester with a thermophilic 
medium, it is 5—10%. In the first case, the rate of biogas production is almost twice as high as in the case of a single 
digester with a mesophilic medium. In the second case, the rate of biogas production is almost the same as when using a 
single digester with a thermophilic medium. Thus, the second case is impractical. 

In the first case, the biogas formation rate is somewhat lower than in the case of using two digesters with a 
mesophilic medium in each, but at the same time, the degree of substrate processing is 10-15% higher. In the second 
case the biogas production rate is twice as low as when using two digesters with a thermophilic medium in each, which 
again shows the inefficiency of the second case. 
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